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qj • The energy eigenvalue problem of non-relativistic positronium is considered 

i—i, within the framework of Relativistic Schrodinger Theory (RST), and the results 

are compared to those of the conventional quantum theory. For the range of princi- 

^ ! 

pal quantum numbers n = 2, 3, . . . , 30, the RST predictions for the non-relativistic 

positronium energies deviate now from the corresponding predictions of the conven- 

— . tional quantum theory at an average of (roughly) 3%. These results suggest that 

I the deviations will be further diminished in the higher orders of approximation. 

■ The emphasis aims at the role played by the assumption of spherical symmetry 

^ • 

of the gauge potential. A new approximation procedure is established in order to 
regard also the anisotropic character of the interaction potential. This reduces the 
deviations of the RST predictions from the corresponding conventional results by 50 
percent, whereas the spherically symmetric corrections of higher order do amount 
only to (roughly) 0,2%. Therefore further improvements of the RST predictions may 
be expected by considering more rigorously the anisotropy of both the interaction 
potential and of the wave functions. Such an analysis (however complicated it may 
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be) seems inevitable in order to decide whether perhaps the exact RST predictions 
do practically coincide with those of the conventional quantum theory. 



PACS Numbers: 03.65.Pm - Relativistic Wave Equations; 03.65.Ge - 
Solutions of Wave Equations: Bound States; 03.65.Sq - Semiclassical 
Theories and Applications; 03.75.B - Matter Waves 
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I Introduction and Survey of Results 



The present paper studies the numerical influence of the anisotropy of the electron- 
positron interaction potential on the energy spectrum of non-relativistic positronium. In 
some preceding papers it has been found that (even up to very highly excited states) the 
corresponding RST predictions come close to their conventional counterparts up to some 
10 percent, already by very rough approximation methods being based on the complete 
neglection of the anisotropy of the interaction force. And therefore one wishes to know now 
whether perhaps this unexpected agreement with the conventional predictions becomes 
still better through including the anisotropy effect. Naturally, the physical relevance of 
Relativistic Schrddinger Theory is to be tested via the experimental verification of its 
numerical predictions in the various fields of possible applications (e.g. the presently 
considered positronium system). Concerning this specific case of demonstration, one is 
tempted to suppose that the exact RST predictions could perhaps agree (nearly or exactly) 
with the corresponding predictions of the conventional quantum theory. If such a result 
(i.e. the numerical equivalence of RST and conventional predictions) could be shown to 
hold for the whole field of atomic and molecular physics, this would then imply a partial 
dethronement of the conventional quantum mechanics. But since the latter is rather of 
probabilistic nature whereas RST is essentially a fluid-dynamic theory, it seems adequate 
and desirable to first clarify the philosophical difference of both quantum approaches, 
where Bohr's complementarity principle plays a dominant role. 



1. Particle- Wave Duality and Complementarity 

In the literature, one frequently encounters the viewpoint that the conventional quan- 
tum theory represents a logically perfect system which cannot be further improved by 
"small" modifications, neither concerning the formalism nor its predictive power. Such 
arrogant judgment arose early during the development of quantum theory when Born and 
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Heisenberg maintained 80 years ago that " quantum mechanics is a complete theory; its ba- 
sic physical and mathematical hypotheses are not further susceptible of modifications" [1,2]. 

It seems that this viewpoint represents nowadays the general conviction of the com- 
munity; however it seems to us that some criticism is well-suited here. Firstly, even if 
presently nobody is able to imagine how further improvements of the conventional theory 
could look like, one cannot be sure that future workers in this field will not invent more ef- 
ficient logical systems which on the one hand do include their Born- Heisenberg precursor, 
but on the other hand do also considerably exceed it. And secondly, the present-day form 
of the conventional quantum mechanics as the offspring of the Born- Heisenberg construc- 
tion may admittedly represent a closed logical system but it nevertheless could perhaps 
miss certain features of quantum matter which on principle cannot be grasped by any 
logical system of the probabilistic type! If this latter supposition should be true, these 
additional features (inaccessible to the purely probabilistic world view) would then be 
well-hidden behind the statistical predictions of the conventional theory. 

Concerning now a possible signal for the presence of such an additional structure be- 
yond the purely probabilistic approach, one could perhaps think of the wave-like behavior 
of quantum matter; and this would suggest that for the description of certain (but surely 
not all) quantum phenomena a fluid-dynamic approach would be better suited than the 
conventional probabilistic viewpoint. This particle-wave duality of microscopic matter has 
ever been thought to be the origin of those strange quantum effects which formerly did not 
only bother the fathers of quantum theory [3] but even nowadays do appear as a "mystery" 
for most of the workers in this field [4, 5]. Indeed, despite the widely celebrated concept of 
quantum logic, one nevertheless likes to discuss many quantum interference effects in the 
classical terms of optical imagery [4]. Naturally, such a dichotomic (or sometimes even 
contradictious) view on the quantum world must have provoked attempts which try to 
erect some unifying framework embracing simultaneously both the fluid-dynamic and the 
probabilistic aspects of quantum mechanics. However, it seems that nobody was (or is) 



7 



able to construct such a unified quantum formalism; and therefore most physicists seem 
to accept now Bohr's original idea of complementarity [2, 6]. Or, summarizing this in the 
words of Omnes[7]: the wave logic and the particle logic could not be united to a larger 
consistent logic which contains both of them as sublogics and consequently, according to 
Omnes[7], we have to learn to live with two mutually excluding logics and nevertheless 
can have a true theory! 

It seems that Bohr [6] was the first who grasped this logical dilemma to its full extent 
by establishing his "principle of complementarity" , whereas others weakened or misin- 
terpreted (or even misunderstood) this original idea of Bohr (for a historical account see 
ref. [8]). For instance, Pauli "in his article for the Handbuch der Physik [9] called two 
classical concepts - and not two modes of desription - complementary, if the applicability 
of the one (e.g., position coordinate) stands in relation of exclusion to that of the other 
(e.g., momentum) in the sense that any experimental setup for measuring the one inter- 
feres destructively with any experimental setup for measuring the other. Pauli, as we see, 
in contrast to Bohr, ascribed complementarity to two notions which belong to the same 
classical mode of description (e.g., the particle picture) and not to two mutually exclusive 
descriptions" (ref '.[2], p. 369). 

An other example of distortion of Bohr's complementarity proposal is given by C.F. 
von Weizsacker who said " The complementarity between space-time description and the 
claim of causality is therefore precisely the complementarity between the description of 
nature in classical notions and in terms of the ip function" (ref. [2], p. 369). Thus, von 
Weizsacker in contrast to Pauli, interpretes correctly Bohr's complementarity principle, 
namely as referring to two mutually excluding "modes of description" (or "logics", as 
Omnes terms it [7]); but in contrast to Bohr, von Weizsacker obviously thinks that these 
two (mutually excluding) modes should refer to the classical and the quantum-mechanical 
description of nature. 

However, concerning that notorious concept of complementarity, we do neither join 
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here to the misunderstanding of Pauli nor to that of v. Weizsacker; but rather we cling 
to Bohr's original idea which says that the notion of complementarity must refer to two 
mutually excluding logical systems (or "pictures") which both are to be used for the de- 
scription of quantum (!) systems, i.e. the probabilistic point-particle picture in opposition 
to the fluid- dynamic wave picture of the elementary matter. If Bohr's complementarity 
principle is understood in this way (i.e. more precisely: is understood in the sense of 
Omnes[7]), the logically coherent and causal space-time descriptions of classical particles 
and waves (e.g. "Classical Electrodynamics" [10J ) must be broken up for the transition 
to the quantum-mechanical description into two mutually excluding modes of description 
(or "logics") for the massive particles, namely into the probabilistic point-particle pic- 
ture with its lack of a physical mechanism for the local but "indeterministic" collapse of 
the probability distribution and into the fluid-dynamic wave picture with its unability 
to predict the instantaneous statistical correlations between widely separated regions of 
three-space (see the EPR phenomena [4,5]). This particle- wave duality is in the first line 
thought to refer to the quantum mechanics of the classical massive particles, whereas the 
classical massless fields are left untouched in the first step of quantization. In this context, 
one could perhaps mention in favour of RST that this theory provides a dynamical mech- 
anism for the (otherwise only kinematically describable) exchange phenomena! These are 
brought forth by the exchange potential £> M as that part of the bundle connection 
which is responsible for the non-Abelian character of the theory. 

Thus, each of the two mutually exclusive quantum fragments (i.e. probabilistic particle 
vs. fluid-dynamic wave picture), being left behind from the classical causal space-time 
description, suffers from a certain deficiency if viewn from the side of the original classical 
approach. And consequently one will hesitate to attribute to one or the other of the two 
competing quantum approaches the status of completeness (being mostly claimed in favour 
of the probabilistic particle picture alone, i.e. the "conventional quantum theory" [1]). 
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2. Relativistic Schrddinger Theory 



Accepting now the necessity of two mutually excluding logical systems (albeit only 
as the preliminary state of the art), one might think that these two pictures have been 
elaborated in a comparably symmetric manner during the past development of quantum 
theory. However, it seems to us that just the contrary did occur: the generally accepted 
and almost exclusively applied form of quantum theory is based upon the probabilistic 
particle picture (i.e. the "conventional quantum theory"), whereas a fluid-dynamic de- 
scription of quantum matter has not been tried at all, apart from a few isolated (but in 
the meantime forgotten) attempts. Such a situation may now be taken as sufficient mo- 
tivation in order to establish a fresh fluid-dynamic approach to quantum matter which is 
able to counterbalance the conventional probabilistic particle logic, from both the mathe- 
matical and physical point of view. As such a competitor of the fluid-dynamic type, there 
has recently been proposed the Relativistic Schrddinger Theory (RST), see the precedent 
paper [11] and the other papers cited therein. Indeed, this fluid-dynamic theory differs 
from its probabilistic competitor (i.e. the conventional probabilistic theory) in both the 
mathematical and physical respect: (i) the many-particle systems are described in RST 
by the Whitney sum of one-particle fibre bundles, not by the tensor product of one-particle 
Hilbert spaces (as in the conventional theory); and (ii) the RST wave functions if) (as the 
sections of complex vector or spinor bundles) are used in order to construct physical den- 
sities of charge, current, energy-momentum etc., not in order to construct probabilities 
as in the conventional theory! 

However, the essential point with such two mutually excluding competitors is now that 
they both are expected to describe the same quantum world; and therefore there must 
exist certain intersecting domains of application for which both proposals make definite 
predictions. Clearly, if both theoretical frameworks should pretend to the same physical 
relevance, they are required to produce identical numerical predictions in that intersecting 
domain of applications. 
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3. Test Case: Non-Relativistic Positronium 



One such common field of competition certainly refers to the energy (E) of bound sys- 
tems. This physical quantity emerges in the probabilistic particle picture as the eigenvalue 
of the Hamiltonian H of the considered system 

Hip = Eijj . (LI) 

On the other hand, the fluid-dynamic character of RST lets emerge the total energy 
of any RST field configuration as the spatial integral of the total energy density ^T 00 (r) 

E T = [ d 3 f^T 00 (f) (1.2) 



where the energy density ( - T - ) Too(r) is the time-component of the total energy-momentum 
density ^T^. Now it has already been demonstrated explicitly for the (relativistic 
and non-relativistic) hydrogen atom that the particle-like energy E (II. ip is numerically 
identical to the fluid-dynamic RST energy E^ (1-2); and this identity holds for all bound 
one-particle systems [T2] . 

Therefore the desired comparison of the RST and conventional predictions must refer 
to the many-particle systems which receive a quite different treatment by each of both 
theoretical frameworks. Perhaps the simplest one of this non-trivial class of test cases is 
the non-relativistic positronium system which consists of two oppositely charged Dirac 
particles of the same rest mass M. The corresponding (non-relativistic) Hamiltonian H 
( II. ip of this two-body arrangement reads 

H = + -4- - ^ 6 ' , (1.3) 
2m 2M |r e -f p | 

and the conventional Schrodinger equation (II. ip of this two-body problem can be separated 
into the center-of-mass motion and the relative motion of the electron and positron (see 
any standard textbook about quantum mechanics , e.g. ref. [13]). The conventional energy 
spectrum of the relative motion is essentially the same as for the ordinary hydrogen 
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problem (with infinitely heavy proton) where merely the mass M of the electron must be 
replaced by the reduced mass M/2; i.e. the conventional non-relativistic energy spectrum 
of positronium is given by 



where n c = 1,2,3,4,5... is the principal quantum number of the conventional theory. 
On the other hand, a first tentative RST calculation on the basis of a two-parameter trial 
amplitude has reproduced this conventional spectrum (1.4) up to some few percents of 
deviation [TJl [33]. Thus, a first crucial test for the physical relevance of the RST predic- 
tions may be referred to the accuracy of the reproduction of the conventional spectrum 
(1.4) when better approximation techniques are applied in RST; and this is just the main 
goal of the present paper. 

However, as will readily become evident, the RST treatment of positronium is much 
more complicated than is its counterpart in the conventional theory, even in the non- 
relativistic domain. But this complication is not an intrinsic feature of RST but is rather 
a consequence of the (purely technical) fact that exact solutions of the RST eigenvalue 
problem are very hard (or even impossible) to obtain. The RST eigenvalue problem 
itself is a simple matter from the purely conceptual viewpoint since it consists merely in a 
simple (albeit coupled) system of two equations, i.e. the Schrodinger-like equation (III. 34) 
and the Poisson equation (III. 36). But despite its structural simplicity, this system needs 
considerable work in order to construct approximate solutions (see below). In contrast 
to this, the conventional eigenvalue problem (1.1) plus (1.3) admits exact solutions which 
can easily be worked out because the internal (i.e. relative) motion agrees with the 
standard hydrogen problem whose solutions can be read in any textbook of elementary 
quantum mechanics. But in view of such a fortunate situation with the conventional 
theory, on should not forget the fact that the conventional eigenvalue problem (1.1) plus 
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(1.4) 
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(1.3) cannot be deduced from a generally valid relativistic quantum mechanics for many- 
particle systems (not to be confused with quantum field theory!). It is true, the well-known 
Bethe-Salpeter equations are intended to represent such a relativistic quantum mechanics 
for many-particle systems, but they seem to be afflicted by many deficiencies and are 
therefore not generally accepted [151 HE] • Presumably, these deficiencies in the relativistic 
domain must necessarily emerge as a consequence of the probabilistic nature of those 
equations. 

On the other hand, RST provides a consistent logical framework for relativistic N- 
particle systems where the particles interact via the principle of minimal coupling which 
then identifies RST as a relativistic gauge field theory (being Abelian for different par- 
ticles and Non- Abelian for identical particles). Within such a framework, one can easily 
deduce a well-defined non-relativistic limit of the relativistic wave equations but the cor- 
responding non-relativistic interaction potential depends now on the quantum state of 
the two-body system (in contrast to the state-independent Coulomb potential (11.31) ); and 
additionally the interaction potential is now also anisotropic (in contrast to the spheri- 
cal symmetry of the conventional Coulomb potential (1.3)). This entails that we have to 
solve simultaneously both the wave equations for the material particles and the gauge field 
equations for the interaction potentials! As a consequence, such simple exact solutions 
as for the conventional Coulomb problem (I.l)-(1.3) are not available for the RST case; 
and therefore one has to put much effort in the development of appropiate approximation 
techniques. This will require a large part of the paper. 



4- Spherically Symmetric Approximation 

But here we can resort to the fact that the coupled matter and gauge field equations 
may be obtained by extremalizing the RST energy functional -E[t], both in the relativistic 
and non-relativistic situations. This provides us with the possibility to invent trial config- 
urations with a certain number of variational parameters whose values become then fixed 
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through extremalization of the corresponding value of the energy functional Em on just 
those trial configurations. Moreover, this method allows us to try spherically symmetric 
configurations in order to avoid the necessity to deal with anisotropic interaction poten- 
tials. Of course, one might suppose here that such a spherically symmetric approximation 
could eventually produce unacceptable large deviations from the true RST result. But 
amazingly enough, the contrary seems to be true: by use of hydrogen-like trial amplitudes 
of SO (3) symmetry one gets already in lowest approximation order RST predictions for the 
non-relativistic positronium spectrum, which come close to the conventional prediction 
(1I.4|) up to 10% or less [HJ E], even for the highly excited states. 

From these tentative results we conclude that the RST predictions could eventually 
come still closer to their conventional counterparts if one would admit anisotropic trial 
configurations and would also apply trial amplitudes with a larger number of variational 
parameters. Thus, there naturally arises the question which of both effects is the dominant 
one: Is it more effective to improve the spherically symmetric approximation by resorting 
to trial amplitudes with more than two variational parameters; or is it more promising to 
take into account the anisotropy of the interaction potential? Or are both effects of the 
same order of magnitude? 

The answer will turn out in favour of the anisotropic corrections which do diminish 
the RST deviations up to 50%; whereas the higher-order corrections of the spherically 
symmetric type amount to merely 0,2%. 

These results are worked out now in the following arrangement: 

A. RST Fundamentals 

In Sect. II, the RST fundamentals are briefly presented to that extent which is neces- 
sary for the subsequent discussion of the positronium level system. Here the point of de- 
parture is the coupled system of matter equations (IH.ip and gauge field equations OII.20p . 
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By splitting off the usual exponential time-factor from the Dirac four-spinors ip a {f,t), 
cf. (IH.lOj) . one gets the mass eigenvalue equations which are written here in Pauli form, 
i.e. in terms of 2-spinors <p±{r), see equation (III. IT)) below. The corresponding time- 
independent form of the gauge field equations turns out as the ordinary Poisson equation 
for the electrostatic interaction potential ^Ao(f), see equation (II 1. 2 7 j) below. Thus the 
combination of the mass eigenvalue equation (III. 171) and the Poisson equation (III. 271) just 
represents the RST eigenvalue problem for positronium. 

At this stage, the eigenvalue system is still fully relativistic with the mass eigen- 
value M* to be determined in one step together with the Pauli eigenspinors ^<p±(r) and 
the gauge potential ^A (r). However, these positronium eigenvalue equations do not 
represent the most general situation, for (i) the magnetic effects are neglected elec- 
trostatic approximation) and (ii) the exchange interactions are also omitted since this 
kind of force can be active only for identical particles. Of course, this absence of the 
magnetic and exchange interactions does considerably simplify the eigenvalue problem 
which, however, is still left too complicated in order that exact solutions could be found. 

B. Spherical Symmetry 

But fortunately, there does exist an energy functional Em (111.32ft whose extremal 
equations do just coincide with those coupled mass eigenvalue and gauge field equations 
which define the RST eigenvalue problem. The corresponding variational principle (111.400 
is the RST principle of minimal energy which can now in Sect. Ill be exploited for 
obtaining approximate solutions for the considered positronium eigenvalue problem (~> 
spherically symmetric approximation) . This approximation method works as follows: 

In order to explicitly solve the mass eigenvalue equations one has in any case to 
parameterize the Pauli spinors by means of both the spinor basis and the associated 
spinor components. The RST hypothesis for this parameterization procedure is now that 
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the fermionic or bosonic character of the bound system as a whole is present already 
in any of its constituents. Since positronium (as a whole) is a boson, both constituents 
(electron and positron) do therefore also adopt bosonic character, i.e. the constituent 
Dirac four-spinors must be decomposed with respect to a bosonic spinor basis which then 
turns out to be double-valued, see equation flIII.5j) below. 

But once the desired decomposition of the Pauli spinors ^(f±(r) into its compo- 
nents {p) ll±, ^S± (i.e. the ' 'wave amplitudes" ) has been attained, see equations (IIII.9a[) - 
(lIII.9bj) below, then one can write down the eigenvalue equations in terms of these wave 
amplitudes, (1III.15aj) - (1III.16j) . This set of equations is now the point of departure for 
separating off the angular part of the wave amplitudes; however, this separation process 
can be performed exactly only if the electrostatic gauge potential ^A (f) is spherically 
symmetric. But on the other hand, this symmetry for the gauge potential can not be 
exactly realized, if the angular separation process for the wave amplitudes is required to 
be exactly feasible, as described by equations flIII.17ap - (IIII.17bj) below. The reason for 
this is that the gauge potential ^Ao(r) must be anisotropic in any case even if the wave 
amplitudes are adopted to exactly obey the spherical symmetry, see the Poisson equation 
(1III.16p . Therefore we have to put in by hand the spherically symmetric approximation 
which consists just in the postulate that the interaction potential ^A (f) (but not the 
wave amplitudes themselves) be spherically symmetric, see equation flHI.lj) below. But 
once this approximation assumption is accepted, it becomes a rather straightforward pro- 
cedure to get the desired approximate energy spectrum (for the present purpose we are 
satisfied with the non-relativistic situation). Namely, after the selection of a nearby trial 
amplitude $(r) (see equation (1III.50I) below) and computation of the associated gauge 
potential ^A (r) (App.A) one substitutes this trial configuration {$(r), ^A (f)} into 
the non-relativistic energy functional Er$i (1III.41P and then one minimalizes the resulting 
energy function E^ Iv l(/3, u) flIII.551) with respect to the two variational parameters (3 and v 
which are contained in the chosen trial amplitude $(r). 
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This minimalization process yields then the energy spectrum displayed in table 1 on 
p. EH from where it is obvious that these RST predictions for the positronium energy levels 
deviate from the conventional predictions E conv (jL4j) mostly by less than 10%. Indeed, 
this is a somewhat amazing result in view of such a simple trial amplitude $(r) 50[) 
with only two variational parameters. Moreover, the deviation from the conventional 
predictions E corw decreases with increasing principal quantum number n c , so that for high 
quantum numbers n c < 100 one arrives at deviations by less than 4%, see table 1 on 
page HH 

Clearly, such a result must evoke now some urgent questions: 

i) Is that RST deviation of some few percent an intrinsic feature of RST (and therefore 
not improvable); or is it merely due to the applied approximation method? In the 
latter case, the exact RST predictions could possibly coincide with their conventional 
counterparts (11.40 . 

ii) If the deviation is due to the approximation method, is this deviation then an 
inevitable consequence of the use of the spherically symmetric approximation or is 
it merely due to the selection of a very rough trial amplitude $(r) (lIII.50j) with too 
few variational parameters? In the latter case it would be worth while to try 0(3) 
symmetric functions with more than two variational parameters. 

iii) If the deviation is due to the use of spherically symmetric trial configurations, can 
then the RST predictions be improved by resorting to anisotropic gauge potentials? 

These questions will receive the following answers by the subsequent elaborations: 

(i) the former 10% deviations jTT] from the conventional results ( II .4p can be pressed 
down to roughly 3% by resorting to the higher orders of the chosen approximation 
method. This supports the expectation that even higher orders of approximation 
will shift the RST predictions further towards their conventional counterparts (11.4|i . 
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(ii) the RST predictions do receive an improvement of merely 0,2% if the trial configu- 
rations are restricted to obey the spherical symmetry 

(iii) The reduction of the deviations from (roughly) 10% to 3% is due to the use of 
anisotropic interaction potentials which therefore provide the possibility for even 
further improvements. 
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II RST Eigenvalue Problem 

From a more philosophical viewpoint, the subsequent relativistic mass eigenvalue prob- 
lem is to be conceived as the RST counterpart of the non-relativistic Schrodinger energy 
eigenvalue problem, i.e. the eigenvalue problem of the Hamiltonian H. The latter problem 
appears as the time-independent specialization of the general time-dependent Schrodinger 
equation when a certain factorization ansatz for the wave function is adopted so that the 
usual exponential time factor can be splitted off. The solutions of the remaining time- 
independent equation describe then the groundstate and the excited states of the consid- 
ered bound system. The same logical arrangement does apply also to the corresponding 
RST eigenvalue problem where, however, the gauge field equations for the interaction 
potentials of the material particles must complete the matter wave equations! It should 
be obvious that such a highly interactive system will necessitate to solve simultaneously 
for both the matter fields and the interaction potentials. Since exact solutions of such an 
intricate system are very difficult to obtain we resort to an approximation method. This 
consists in adopting some physically plausible trial function which contains a certain set 
of variational parameters serving for extremalizing the RST energy functional. Indeed, 
the extremal equations of that energy functional are identical to the system of coupled 
matter and gauge field equations forming the RST eigenvalue problem. 

Obviously this program requires to discuss the following items in due order: (i) wave 
equations for the matter subsystem, (ii) equations for the gauge field subsystem, (iii) 
stationary bound systems, (iv) energy functional and principle of minimal energy, and 
(v) spherically symmetric vs. non-symmetric trial functions. 
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1. Matter Subsystem 



Quite generally, the dynamics of the spinning quantum matter is described in RST by 
the N-particle Dirac equation 

ihcT^V^ = Mc 2 ^ . (II. 1) 

Here, for the considered two-particle system (i.e. positronium), the two-particle velocity 
operator is the direct sum of the ordinary Dirac matrices 7 M 

r M = (-7 M )©7/, (n.2) 

where the minus sign refers (by convention) to the "first" particle (i.e. the positron). 
Furthermore, the two-particle wave function \1/ is adopted in RST to be the Whitney sum 
of the two one-particle wave functions ip a (a = 1,2) 

(r, t) = Vi(r, t) © ip 2 (r, t) (II.3) 

which says that each of the two particles does occupy a well-defined one-particle quantum 
state. 

The gauge-covariant derivative (T>) emerging in the basic equation ( III.ip is defined as 
usual in the gauge theories, i.e. 

V l $ = dJ$ + A l $ , (II.4) 

where the gauge potential takes here its value in the Lie algebra of the structure group 
U(2) and may therefore be decomposed as follows 

2 

A ti = Y, A \ T - + B ^ - B,x ■ (IL5) 

a=l 

Strictly speaking, this general form of the gauge potential does apply for a system of 
identical particles which are subjected to the exchange forces (being described by the 
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exchange potential B^). But since the presently considered positronium system consists 
of two different particles (i.e. positron and electron), the exchange potentials B^.B^ as 
the components of A^ with respect to the exchange generators x, X must vanish (B^ = 0). 
Thus we are left alone with the purely electromagnetic interactions which are described 
by the electromagnetic potentials A a M (a = 1, 2) as the components of with respect to 
the electromagnetic generators r a : 

2 

a^J2 A v<> = {1)A » • ri + {2)A » • r2 • ( IL6 ) 

a=l 

And finally, the mass operator Ai in equation (III. II) can be taken to be proportional to 
unity (1) for the present situation where both particles have the same rest mass (M, say) 

M = M ■ 1 . (II.7) 

In order to further proceed towards the time-independent eigenvalue equations, one 
first decomposes the original two-particle equation (IH.lj) into two one-particle equations 
for the one-particle wave functions ipiif^t) and if)2(r,t): 

ifvc^D^x = -Mc 2 ^ (II.8a) 
ihcYD^ 2 = McV 2 • (II.8b) 

Here, the gauge- covariant derivative (D) for each wave function ip a (a = 1,2) is defined 
as follows: 

D^ = ^-i (2) VV>i (n.9a) 

D li <fo = d li ik-iWA ft -<fo (II.9b) 

which expresses the fact that the two different positronium constituents are not able to 
feel the exchange forces (~* B^ = 0) but are subjected exclusively to the electromagnetic 
interactions! 

Next, one has to face the problem of time-dependence of both wave functions ip a (r, t) 
where one naturally will expect that the relevant physical objects must be time-independent 
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for a stationary bound system. In this sense, one tries the usual factorizing ansatz for the 
one-particle wave functions ip a {f,i) (a = 1,2) 

/ M c 2 \ 

xf> a (r,t) = exp l-i-±-t\ .^ a (r) , (11.10) 

where the mass eigenvalues M a are to be determined just by solving the mass eigenvalue 
problem below. Furthermore, one observes the fact that the Dirac four-spinors i[> a (r) may 
be conceived as the Whitney sums of Pauli two-spinors ^ a '(f±(r), i.e. we put 

^ a (r)= ( a V+(r)© (a V-(r) ; (11.11) 

and for these Pauli spinors one deduces from the original Dirac equations (III.8al) - (III.8bl) 
the following mass eigenvalue equations [17] 

iff- V^ViW + (2) ^o(r) ■ (1 V T W = ±M + M * c- (1 VtW ( n - 12 a) 

2<T- V (2 V±(f) + (1) A (f) ■ ( 2 Vt(^) = - M *± M C . (2)^( r -) . (II.12b) 

Here, the common rest mass is denoted by M, cf. (1X171) . and the common mass eigenvalue 
by M* (= —Mi = M%). Furthermore, we have also neglected the magnetic potentials (i.e. 
putting A a (r) = electrostatic approximation) . An essential point with this neglection 
of magnetism aims at the well-known ortho/para dichotomy of positronium [TT] . Namely, 
it should be self-evident that the physical difference of ortho- and para-positronium be- 
comes annihilated through disregarding the magnetic interaction energy, see the discussion 
of this in ref . [TI] . But observe on the other hand that the present mass eigenvalue system 
flII.12aj) - flII.12bj) is still of truly relativistic nature, though it is written in terms of Pauli 



2-spinors, not in terms of Dirac 4-spinors! 

A further crucial point of the eigenvalue system fill. 12a]) - fill. 12bl) concerns the plausible 
assumption that, from symmetry reasons, both particles (i.e. electron and positron) will 
occupy the same physical quantum state (for a more detailed discussion see ref. [TTj ) . A 
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nearby consequence of this assumption is that the Dirac four-densities ^k^(r) must be 
the same for both particles 

(%(f) = (%(f) = <P\(f) , (11.13) 

provided the two particles do combine to the para-state (for the case of ortho-positronium 
see ref.s [Hi E]). On the other hand, the Dirac densities ^ko(r) read in terms of the 
Pauli spinors ^Lp±{r) ( III. lip 

^k Q {f) = (a V+(r) (a V+(*0 + (a V-(^) (a V-(r) 
and similarly for the Dirac currents k a {f) 

k a (f) = ( a Vi-(r)^ (a V-(0 + (a V-(^)^ (a V+(^) • (11.15) 

Therefore, the physical equivalence of both one-particle states allows us to parametrize 
both states by only one set of Pauli spinors ^(p±(r) which then lets appear both original 
spinors ^(p±(r) and ™(p±(r) as follows: 

&<p + (r) = (1 V+(r) = % (k' a) {2) <f+(r) (II.16a) 

(p)^_(f) = ^V-(r) = i (2 V-(r) (II.16b) 

k r , 

k 



l^pl 



By use of the usual spinor algebra one can easily verify that the algebraic require- 
ments (III.13P are actually satisfied through the present arrangements flII.16a[) - (lII. 16b|) : 



and additionally the two equations (III. 12a,|) - (III. 12b|) for the individual spinors ( 1 '<p±{f) 
and ^ip±{r) become reduced to one equation for the para-spinor ^tp±(r): 

iff- V (p V±(r) - (p) ^o(r) ■ (p VtW = M * f M c " (P Vt(0 ■ ( IL17 ) 

Here, both electrostatic potentials ^A (r) and ^A (r) have also been identified up to 
sign 

w A (f) = - {2) A (r) = (p) A (f) (11.18) 
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because they are generated by the same Dirac density ^k (r) 

(P)jfc (f) = ^k (r) = (2) A; (f) . (11.19) 

The precise form of the link between the potential ( p ) Aq (r) and the Dirac density *- p ) ko (r ) 
must now be deduced from the original gauge field equations. 



2. Gauge Field Subsystem 



The gauge field counterpart of the basic matter field equation flll.lj) is the (generally 
non-Abelian) Maxwell equation 

V^ u = -Ama s J v (11.20) 

with the curvature T^ v of the bundle connection being defined as usual in the 

gauge field theories 

2 * (11.21) 

o=l 

However, for the present situation of different particles the exchange potential van- 
ishes identically (B^ = 0) and consequently the theory becomes Abelian with the two 
field strengths F a ^ u being defined in terms of the residual electromagnetic potentials A a ^ 
through (a = 1,2) 

F\ v = V^A\ - V V A\ . (11.22) 

Furthermore, since we are presently satisfied with the electrostatic approximation, 
the relativistic relation (1H.22|) reduces to a simple gradient link of the electric field 
strengths E a {f) to the electrostatic potentials ^A (f) 

E a (r) = -V^A (r) . (11.23) 
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From the same reason, the original Maxwell equations fill. 20 j) do reappear as simple source 
relations 

V • E a = 4vra s • ^ Jo (r) (11.24) 

where however the Maxwell densities j a „ = { ^jo, —ja} are related to the Dirac densi- 
ties k ail through 

j\ = k lfl = {^k ,-h} (II.25a) 
j\ = -k 2 » = {- {2) k ik 2 } . (II.25b) 

Thus the source equations (111.24j) ultimately appear as the well-known Poisson equations 

A «A = -4vra s (1) A; (II.26a) 
A {2) A = 47ra s {2) k , (II.26b) 



which however are contracted to only one equation: 
A (p) A (f) = -47ra s • (p) fc (f) 
= — 4Tra* 



(11.27) 



(P VV (r) (p V+(r) + (p V-(r) (p V-(r) 

namely as a consequence of the former identifications (1II.13j) and fill. 18 j) . 

Summarizing, the RST eigenvalue problem for positronium in the electrostatic approx- 
imation consists of the mass eigenvalue equation (1II.17P in combination with the present 
Poisson equation (111.271) . It is true, this is a closed system for the Pauli spinors ^(p±(r) 
and the interaction potential ( p ^Ao(r) but it demands the specification of certain boundary 
conditions. Clearly, for the matter fields one demands their sufficiently rapid vanishing 
at infinity (r — > oo), i.e. 



lim (p V±(r) = lim {p) k {r) = 



(11.28) 



and moreover one wishes also to have thereout the asymptotic Coulomb form for the 
interaction potential ^Ao(r) 

lim (p) A (f) = — . (11.29) 
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Both conditions (1II.28l) - (lII.29j) lead us to the standard solution of the Poisson equation 
(1X1271) . i.e. 

(p) A,(r) = a s [ d 3 f'^Q , (11.30) 

i 7" 7" 1 1 

where the normalization condition for the Pauli spinors 

fjpf (p)k (r) = 1 (11.31) 



(cf. (III. 14ft ) actually ensures the asymptotic Coulomb form ( 1II.29j) of the potential ^Ao(r) 

(HEP. 



3. Energy Functional 

Even if one would be able to find exact solutions of the eigenvalue equations fill. lTj) plus 
(III. 271) . with both constraints (111.29 j) and (III. 311) being obeyed, one nevertheless would be 



forced to face the problem of the energy content carried by that elaborated solutions. In 
other words, one needs an energy functional (J5[t] , say) whose value upon the constructed 
solution yields its physically relevant and observable energy. Surely, the wanted energy E T 
cannot be identified with the mass eigenvalue M*c 2 because this quantity refers separately 
to each of both particles and therefore is a one-particle quantity whereas we would like to 
know the total energy of the interactive two-particle system. It would also not help to 
take 2M*c 2 as the desired energy Et because M*c 2 alone does already contain the whole 
electrostatic interaction energy which then would be counted twice. 

The solution of this energy problem has been already worked out in some precedent 
papers [TT] and for the sake of brevity it may suffice here to simply quote the result: 

E [T] = E™ + 2A D • N D + A§> • . (11.32) 

This general form of the energy functional says that there are essentially two parts which 
are equipped with a rather different meaning: the first part is a collection of the 
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truly physical energy contributions 

= 2Mc 2 -Z 2 p + A ( p )T kin + 4 e) (11.33) 

while the second part does refer to the constraints. 

Here the first constraint (Nd) refers to the wave function normalization (III. 3 1 [) : 



N D = J d 3 f ( (P V+ (P V+ + (p V- (P V-) -1 = 0, (11.34) 

and Ad is the associated Lagrangean multiplier. Similarly, the second constraint (Nq^) 
does refer to the Poisson identity 

n£ ] = - M (e) c 2 = (11.35) 

with the gauge field energy being defined through 

4 e, = -^/^(£ p .i p ) (11.36) 

and its mass equivalent (M^c 2 ) through 

M (e) c 2 = -he [ d 3 r {p) A (f) • {p) fco(f) . (11.37) 



The Poisson identity ( 1II.35P itself is an immediate consequence of the Poisson equation 
(III. 271) and therefore is automatically satisfied not only for the exact solutions of the 
RST eigenvalue problem but also for its approximate solutions (see below) if only the 
approximate potential ^Ao(r) is linked to the (approximate) Dirac density ^ko(r) via 
the Poisson equation (III. 271) . Thus, if one takes care of satisfying simultaneously both 
constraints fill. 341) and (III.35P the energy functional Em fill. 321) becomes reduced to its 
physical part E^ (III. 331) . 

This latter part (III. 33|) represents the proper physical content of the energy func- 
tional Em and consists of the (renormalized) rest mass energy, the kinetic energy of both 
particles and the gauge field energy 

£r } (|IL36|) . The mass renormalization factor Z a is 

given by 



Zl = / d 3 r ( ( p V! (p V+ - (l V- (P V-) , (H.38) 
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and the one-particle kinetic energy ^T^ in reads 

V+(r) (a* 



(p),„t ^W^. v) (p V-(r) + (p V-(^) (<?• V) (p) 



. (11.39) 



After all, it is a nice consistency check to carry explicitly through that standard recipe of 
variational calculus as it is required by the principle of minimal energy 



5E\ 



[T] 



. 



(11.40) 



Or in other words, the extremalization of the energy functional Em (111.320 with respect to 
the Pauli spinors ^<f ±, ^ P V± yields just the mass eigenvalue equations (111.170 as the first 
part of the RST eigenvalue problem; and analogously the extremalization of Em (III. 321) 
with respect to the electrostatic potential ^A (f) lets one recover the Poisson equation 
(III.27I) . provided both constraints (III. 34j) and (III.35P are duly regarded. 
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Ill Spherically Symmetric Approximation 



From the preceding presentation of the RST eigenvalue problem it should be evident 
that it is (almost) impossible to construct exact solutions of such a peculiar kind of eigen- 
value problem. On the other hand, it is perhaps not even necessary to know those exact 
solutions since their main features could possibly be grasped already by studying certain 
approximate solutions. A nearby simplification as the basis for such an approximative 
approach refers to those anisotropic effects which are invading the calculations via the 
spin degree of freedom inherent in the Dirac equations. In order to get rid of this type 
of complication, one may neglect the anisotropic influence of the spin on the electrostatic 
interaction potential and may assume that the interaction potential ^A (r) is nearly 
spherically symmetric, i.e. we put 



where this assumption of spherical symmetry needs not be applied to the wave amplitudes 
of the particles themselves. For such a first estimate of the positronium spectrum (on the 
level of accuracy of the conventional treatment (11.4)1 ) it is also not necessary to stick to the 
relativistic formulation of the theory, but rather one may again be satisfied with the non- 
relativistic limit. Thus we will first eliminate the explicit spin anisotropy effect from the 
eigenvalue problem and afterwards we will resort to the non-relativistic approximation. 



The manifest spin degree of freedom may be eliminated by first parametrizing the 
Pauli spinors ^ P V±(^ through the wave amplitudes ^TZ±(r) and ^5±(f) as the com- 
ponents of ^ip±(r) with respect to a selected spinor basis; next, one separates off their 
angular parts and finally one impresses a rigid link on the remaining radial parts of both 



(p)Ao(f) => [p U (r) 



(III.l) 



r = | \r\ 




1 . Double- Valued Spinor Fields 
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amplitudes ^1Z±(r) and ( p >S±(r). This procedure yields (in its non-relativistic version) 
a scalar eigenvalue problem of the Schrodinger type, see equation (1III.34j) below. 



The selection of an appropriate spinor basis starts with the standard eigenspinors (j ,m 
of total angular momentum J (= L + S) in two-dimensional unitary space 

>cr=j(j + i)^ 2 -cr (Hi.2a) 
l 2 d m = i{i+\)h 2 -cr (m.2b) 

I 2 cr = + ■ cr = \ Q + 1) fr 2 ■ cr 

j z c,r = mh-ci m ■ (in.2d) 

For the sake of simplicity, one may prefer here / = 0, 1 for the orbital part of angular 
momentum so that {Q' 2 } could be a basis for the "positive" Pauli spinors ( p ^ + (f) and 
similarly {Q' 2 } a basis for the "negative" spinors ( p )y>_(r). Now one can show that 
(for para-positronium (III. 16a,|) -( 1II. 16b[) ) the action of the z-component J z of total angular 
momentum J must annihilate both spinors ^<p±(r), i.e. 

j(+) <P)<p + (f) = (III.3a) 
jH (p V-W = 0) (III.3b) 

see ref.s (TTj [T7]. Or in other words, the action of the total operator J z (= © Jz ) 
annihilates the one-particle Dirac spinor field tp p (= ^ P V+ © ^ P V-); i- e - 

M P = o • (IH.4) 

But clearly, such a one-particle state must be of rather exotic nature because an 
ordinary ( "fermionic" ) Dirac eigenspinor of J z has always half-integer eigenvalue; and 
consequently such a strange ( "bosonic" ) state as ( 1III.4P must own some unusual property. 
Indeed this concerns its uniqueness in the sense that the spinor field ip p (IIII.3aj) - flIII.4j) is 
double-valued [HI [T7] 

^M,0 + 27r) = -^(r,tf,0) (III.5) 
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where r, are the usual spherical polar coordinates of flat three-space. On the other 



- ±- - it- 

hand, both basis systems {C 2 ' 2 } and {Q' 2 } flIII.2ajl - flIII.2dll are unique over three- 



space; and thus if one would decompose such a spinor-like %fj p (1III.5P with respect to a 
standard basis, the corresponding components would have to be non-unique scalar fields 
over three-space. However, our option is here just the other way round, namely to choose 
certain basis systems and which themselves carry the double- valuedness (IIII.5I) 
alone so that the components of ip p can remain unique scalar fields, i.e. 

4 ±} (r, 0, + 2tt) = -4 ±} (r, 0, 0) (III.6a) 
(r, + 2tt) = -oo[ ±] (r, 0) . (III.6b) 



Now our choice of basis system looks as follows: 





■ 

— 7 — 

= e 2 


i i 
So 




• 

= e l 2 ■ 


1 1 

1- 2 ' 2 




= e 2 


1 1 

CP 




• 

= e l 2 ■ 


1 1 

<"-2' 2 
Si 



(III.7a) 
(III.7b) 
(III.7c) 
(III.7d) 



so that this w-basis becomes annihilated by J^: 



U M = o (III.8a) 
jH w ; ±} = . (III.8b) 

The double- valuedness (lIII.6aj) - (lIII.6bl) of the chosen basis system transcribes now to 
the Pauli spinors ( p 'ip±(r), provided one adopts the corresponding components ^ P 'TZ±, 
( pS) S± to be unique. Accordingly, the desired decomposition reads 



(p) 



v+ (f) = (P)ft + (r) • J +) + (p) <S+(r) • co^ (III.9a) 



{p V-(r) 



i { (p) ft_ (r) ■ c^ +) + (p) S_ (f) ■ } (III.9b) 
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with unique scalar components 

(p) ft±(r,tf,0 + 27r) = {p) ft±(r,tf,0) (111.10a) 
(p) <S ± (r,tf,0 + 27r) = (p) S±(r,tf,0) . (III. 10b) 



2. Eigenvalue Equations for the Wave Amplitudes 

Further information about the components ( pS, lZ±(f) and ( p ^S±(r) can be gained by in- 
spection of how the annihilation process flIII.3a[) - (1III.3b[) can be realized in detail. Indeed, 
a straightforward calculation yields 

( P ) Mrl = {Lz (p)^ + ).4 + ) + (L 2 Ms + )-ujT ) + (p) ^ + -(Ji +) 4 +) )+ 

(111.11) 

and analogously for Wyj_(f). Observing here the annihilation relations for the w-basis 
(lIII.8al) - (IIII.8bj) it is clear that one has to demand 

L z (p) TZ+ = L z (p) 5+ = (III. 12) 

in order to have the eigenvalue equations (1III.3a[) - (1III.3b[) for angular momentum satis- 
fied. However, the latter demand (1111.120 can trivially be satisfied by letting the wave 
amplitudes ^TZ±(r), ^S±(r) depend exclusively upon r and but not on <p. Thus one 
puts 

(p)ft ± (f) ^R±{r,$) (III.13a) 
V>S±(t) => (p) S±M) . (III. 13b) 

Furthermore, it turns out that the equations for the wave amplitudes become consid- 
erably simplified if one resorts to a further transformation: 

{p) R ± (r, 0) = (P) R± ^ g ) (III. 14a) 

V r sin § 

U>S ± W)= iP) f ±M . (III.14b) 

V r sin i? 
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Namely, the relativistic Pauli equations ( 1II.12aj) -( 1II.12b|) are then recast to the following 



relatively simple eigenvalue equations for the new amplitudes ^R±, ^S± p] 

(III.15a) 
(III. 15b) 
(III.15c) 





dr 


+ + 

r d"d 


- (P) A • 




M + M* 
= ft C " 


(p) i?_ 




dr 


ld^R + 
r dd 






M + M* 
= ft C " 


(p)£_ 


d^R_ 
dr 


+ - • Cp)«_ 
r 


ld^S_ 
r dd 


+ (p) A) • 




M-M* 
= ft C " 


( p )£ + 


d^S_ 
dr 


+ 1 • ( p )S_ 
r 


id®iL 

r dd 


+ {P) A 


{p) s+ 


M - M* 
= ft C " 


(p) 5+ 



(III.15d) 

In order to deal with a closed system of equations for the new amplitudes, one substitutes 
the Pauli spinors (jIII9aP -( lIII.9b|) into the Poisson equation ( 1II.27P which lets reappear 
this equation now in the following form 

(P)»2 i ( P )C2 , (p)62 , ( P )C2 

A ^A = -a s K++ b++ q - + ^ . (111.16) 

r sin 17 



The coupled system 15a[> - 16j) represents now the RST eigenvalue problem for 



para-positronium in the electrostatic approximation; but despite its simplified form it is 
hard (or even impossible) to elaborate the exact solutions hereof. 



3. Spherical Symmetry 

In this situation, the approximative assumption of spherical symmetry (IIII.ip allows 
further progress. Namely, it is just this assumption which suggests to try a product ansatz 
for the wave amplitudes of the following form [TT] 

^R ± (r,$)= 'R ± {r) ■ f R (0) (III.17a) 
^S ± (r,$)= 'S±(r)-f s (#) . (III.17b) 
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Indeed, substituting this ansatz into the mass eigenvalue equations (1IH.15atlHI.15dl) allows 
one to separate the variables so that there emerges one subset of equations for the radial 
functions 'R±(r), f S±(r) and one subset for the angular functions and fs($)- Here, 

the angular equations introduce the quantum number l p of orbital angular momentum in 
the following way: 



dfs(#) 
d$ 



l p ■ f s (0) (III.18a) 
-Wfl(0) , (III.18b) 



where a second possibility does exist which, however, has merely the sign of l p reversed 
(i.e. l p =>■ —lp). The separated angular problem flIII.18aj) - flIII.18bj) is so simple that its 
solutions are immediately evident, namely either 

f R (&) = cos(/ p • 0) , f 8 (#) = - sin(/ p ■ #) (111.19) 

or 

f R {&) = sin(/ p • t?) , f s {0) = cos(/ p • 0) , (IH.20) 
with the quantum number l p adopting (half-) integer values 

ip = 0, {\l 1,(^,2,(^,3, ... (111.21) 

This postulate of (half-)integrity follows from the demand that the Dirac current k p re- 
mains finite on the z-axis [d = 0, 7r). For the exclusion of the half-integers see below. 

After the angular part is split off in that way (1III.17al) - (1III.17bl) . there remains a purely 
radial problem for the determination of the mass eigenvalue M*. Indeed, substituting the 
separation ansatz (lIII.17ap - (lIII.17bl) back into the original eigenvalue equations 15a|) - 
( HII.15dl) and using also the angular equations (IIII.18al) - (IIII.18bl) yields the following set 
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of coupled radial equations: 

d'R + (r) l p 
dr r 
d'S+(r) L 



dr 

d'R-(r) 1 ,„ , 



dr 



's + 


(r) - [p U ( 


r) ■ 




(r) = 


M + M* 

n c ' 


'R-(r) 


(III.22a) 


'R + 


(r) - [p U 


(r) 




(r) = 


M + M* 

n c ' 


f S-(r) 


(III.22b) 




(r) + t p U ( 


r) 


ft* 


(r) = 


M — M* 
ft C ' 




(III.22c) 




(r) + [ p L4 


(r) 


'S + 


(r) = 


M-M* 
h C ' 


f S + (r) . 


(III.22d) 



dr r r 

Surely, it appears somewhat strange here that one and the same eigenvalue M* enters 
all four equations (IIII.22a|) -( ITTL22d|) : but this may be understood clS cL hint at the cir- 



cumstance that there exists a rigid link between the four variables 72±(r), 'S±(r). And 
indeed, the following identifications 

'R+(r) = 'S+(r) = <§+(r) (III.23a) 

(R_(r) = &_(r) = <£_(r) (III.23b) 

recast this system of four equations to only two equations 

- l ± ■ $+(r) - [ p U (r) ■ <£_(r) = c • $_(r) (III.24a) 
dr r n 

+ k±l . $_( r ) + [pU ( r ) . $ + ( r ) = M ~ M * c- $ + (r) . (III.24b) 
dr r a 

This spherically symmetric eigenvalue problem may finally be completed by the cor- 
responding Poisson equation 



which is to be conceived as the spherically symmetric simplification of flIII.16j) . Notice here 



also the important fact that the quantum number l p does not explicitly enter the source 
terms of the spherically symmetric Poisson equation (1III.25j) ! The standard solution of 
this type of equation reads formally 

/fg»+M+f-W ( ,n. 26) 
8 J r' \\r — r '\\ 
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and approaches the Coulomb potential at infinity (r — > oo), i.e. 

lim [p U (r) = — , (111.27) 



namely just on account of the normalization condition fill. 31 j) which reads in terms of the 
wave amplitudes $±(r) 



dr'r' <§ 2 Jr') + $i(r f )\ = 1 . (111.28) 



Clearly, such a requirement can only be satisfied if the amplitudes $±(r) obey the usual 
boundary condition 

lim <£ ± (r) = . (111.29) 

r— >oo 

However, a somewhat more critical point concerning the boundary conditions refers 
to the behaviour of the angular functions /r($) and fs($) (1III- 19[) - (1III.20[) on the z-axis 
(i.e. $ = 0,7r). From generally accepted arguments one would think that all physically 
relevant objects of the theory should be non-singular everywhere. In the present context, 
the crucial object is here the Dirac current k p which for the presently considered kinematics 
of the Dirac spinor field ip p = ( - p V+ © ^ P V- encircles the z-axis 

Mr)= (p) fyM)-e* (111.30) 

with the azimuthal component being given by 

» ttM ) = i i|f>i - |r) ." [(i ''V Hl . (111.31) 
2rrr sin v 

Therefore, if one wishes to keep this current component finite on the z-axis then one is 
forced to demand that l p can adopt only (half-)integer values, cf. (jIII.2ip . Furthermore, 
if one demands that ^'k^ should be distributed symmetrically along the z-axis, i. e. 

(P)^(r,0) = (p) ^(r,7r) Vr , (111.32) 

then one can admit only integer values of l p 

l p = 0, 1, 2, 3, ... (111.33) 
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(For a sketch of for integer and half-integer l Pl see fig. 1 of ref. [11]) . Henceforth we 
will prefer the present quantization prescription (1III.33j) over the former one [TT]). 

But with the eigenvalue and Poisson equations flIH.24aP - flIH.25p being fixed now to- 
gether with all the boundary conditions, one could in principle look for the solutions of 
this eigenvalue problem and determine the corresponding binding energies. However, it 
seems meaningful to first resort to a further (but ultimate) approximation. 



4- Non-Relativistic Approximation 

Besides the electrostatic and spherically symmetric approximations one can now take 
into account a final simplification, i. e. the non-relativistic limit. The point with this 
approximation is namely that the relativistic corrections may well be of the same order 
of magnitude as the magnetic effects; and therefore it seems somewhat inconsequent to 
neglect the magnetic effects on the one hand but on the other hand to cling to a fully 
relativistic formulation. Thus one should be willing now to pass over to the non-relativistic 
approximation of the relativistic eigenvalue system (IHI.24ap - flIH.25p . Here, the two mass 
eigenvalue equations (1III.24a|) -( 1III.24bT) contract to only one equation for the "positive" 
Pauli spinor $+(r) which for the sake of brevity is simply termed as $(r) and has then 
to obey the following Schrodinger-like equation: 

- m + l^} + " ' *w - E - ' *« • c^> 

The Schrodinger eigenvalue E* emerging here is related to the former relativistic mass 
eigenvalue M* by 

E, = (M* - M)c 2 , (111.35) 

and the meaning of l p as a quantum number of orbital angular momentum becomes now 
evident. In order to have the non-relativistic eigenvalue system closed, the Schrodinger 
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equation ( 1III.34|) must be accompanied by the non-relativistic version of the Poisson equa- 
tion fjIIL25|l . i.e. 

^+^r u °«=-2 a -^ i - < iil36) 

Correspondingly, the non-relativistic approximation of the former solution ( 1III.26I) for 

[ p l4 (r) appears now as 

®Mr) - a ij^0^ (HI.37) 

and thus continues to obey the same Coulomb-like boundary condition flIII.27j) . namely 
on account of the non-relativistic form of the normalization condition flIII.28j) 

- J drr $ 2 (r) = 1 . (111.38) 

In this way, one finally has arrived at a well-defined non-relativistic eigenvalue problem. 

But even if one were able to solve exactly this radically simplified problem, one would 
still be left with the problem of the energy (Et, say) being concentrated in this field 
configuration. This is the question of the right energy functional Em whose value upon 
the solutions of the eigenvalue problem is required to not only yield the desired energy 
Et but whose extremal equations are also required to just coincide with the eigenvalue 
equations: 

5E [T] = . (111.39) 

This principle of minimal energy can then also be used in order to construct approximate 
variational solutions of the eigenvalue problem; namely by guessing some trial configu- 
ration with a certain number of variational parameters {/3;i>i,V2, ■■■Vk}- The value of 
the non-relativistic version Em of the relativistic energy functional Em upon this trial 
configuration yields then an energy function E^ Iv ^(/3; u±, V2, ■■■ Vk) whose (local) extremal 
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values are determined by the vanishing of all the first-order derivatives 

dE^(f3;u u u 2 , ...v k ) 



df3 

dmW(f3;is l7 v 2 , ...u k ) 
dvi 







(111.40) 



dE^(f3;u u u 2 , ...u k ) =Q 
dv k 

The corresponding extremal values e£ do then specify the energy levels of the bound 
system, albeit only in the non-relativistic, electrostatic and spherically symmetric approx- 
imation! 

The desired energy functional E[tj is the non-relativistic version of the former E\^\ 
(III.32p - flII.33p and has already been deduced in full detail in the preceding paper (see 



equation (VI.50) of ref . [II]) 

E [T] E m = 2E Un + E§ + 2A S N* + ^ ■ (111.41) 

Here, the first term represents the kinetic energy for the two particles where the one- 
particle energy E^in is the sum of the radial part ^E kin and the longitudinal part ^E kin , 
i. e. 

£ ki n = (r) £ kin + w E kin , (111.42) 
with the radial part being given by 

1 r drr f*wV (,,,.43) 

km 2M 2 J y air J y ' 
and similarly the longitudinal part by 

W Ek . n = —l 2 -- f° - $ 2 (r) . (111.44) 

km 2M p 2 L r V ; V ; 
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Next, the energy content of the spherically symmetric potential H4 (r) reads 



and thus the physical part of the functional E[$] ( lIII.41j) consists just of these three energy 
contributions flIH.43jl - pi.45jl . 



The residual two terms in dIII.41[) are constraints with Lagrangean parameters As and 
. The first < 
amplitude $(r) 



Aq\ The first constraint refers to the normalization flIII.38j) of the non-relativistic wave 



poo 

71 f ■ \2, 



Ns = -/ drr$\r) - 1 = , (111.46) 



o 



and the second constraint is the Poisson identity 

Ng = Eg 1 " M [o] c 2 = (111.47) 
with the mass equivalent M^c 2 due to the gauge field energy being given by 

poo 

M [cl c 2 = he I drr [p U (r) • $ 2 (r) . (111.48) 

2 Jo 

A nice check of all these assertions consists now in deducing the present eigenvalue 
system flIH.341) plus 36[) from the non-relativistic version (SKm = 0) of the principle 
of minimal energy (IIII.39[) by means of the standard variational techniques. Thereby the 
Lagrangean parameters turn out as 

A s = -E* (III.49a) 
Ag } = -2 . (III.49b) 

This calculus of variations may now be exploited in order to compute approximately the 
non-relativistic energy spectrum of positronium. 



40 



5. Spherically Symmetric Spectrum 

In order to present such an example for the energy spectrum due to the spherically 
symmetric approximation, one may select as trial amplitude $(r) the following function: 

$( r ) = $, r u e-P r , (111.50) 



where the normalization condition f lIII.46j) fixes the constant to 



* Ti T(2z/ + 2) v ; 

This trial amplitude (1IH.50|) has only two variational parameters (i.e. /3 and u), and 
therefore one cannot expect that the corresponding approximate energy spectrum will 
turn out to be more accurate than up to some few percent (< 10%, say). Nevertheless 
this order of accuracy should be sufficient in order to estimate the magnitude of the 
anisotropic corrections which, plausibly, will be found to be somewhat smaller than the 
error induced by just the spherically symmetric approximation. In any case, once the trial 
amplitude $(r) is fixed, one can proceed to determine the value of the non-relativistic 
energy functional Em upon this trial configuration which then yields the corresponding 
energy function E^ Iv ^(/3, v) through the following steps: 

(i) Calculate the kinetic energy (lIII.42p - flIII.44p by means of the selected trial 
function $(r) (1III.50p and find this quantity to be of the form 

E kin = ^- (2a B /3) 2 • s^M (HI.52) 
2a B 

h 2 

(0b = • • • Bollr radiuS ) 

with the kinetic function £kin(^) being given by 

1 (\ I 



^) = ^TT( i + ^). (HI.53) 
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Solve the Poisson equation (1III.36j) for the spherically symmetric potential H4o(r) 
(App.A) and compute by means of this result the eloctrostatic gauge field energy 
(1III.45P which then emerges as follows 

4 6] = -— (2/3«b ) • e pot (u) . (HI.54) 

Here the potential function e po t{v) may be expressed in two alternative ways on 
account of the Poisson identity (lIII.47p . see equation f]A.28|) of App.A. (For a 
sketch of the spherically symmetric trial potentials [ p i4 (r) (1III.37p . see Fig. A. I). 

Substituting both energy contributions E^ in (1III.52I) and Eft (1III.54j) back into the 
energy functional Er$i (lIII.4ip yields a certain function E^ Iv ^(/3, u) of the two vari- 
ational parameters (3 and z/, see equation (IA.22I) of App.A: 

2 

E [iv] (/3; u)^—\ (2a B (3) 2 ■ £ki » - (2a B (3) ■ e pot (v)\ . (111.55) 

(Observe here that both the normalization condition (1III.46P for the wave amplitude 
$(r) and the Poisson identity (1III.47P are satisfied, so that the energy function 
E^ ((3, u) due to E[t] (jIII.411) consists solely of the kinetic energy and the gauge 
field energy). 

The extremalization process (lIII.40p can now be carried through for the obtained 
function Et Iv l(/3, u) flIII.55p . Here it is most convenient to perform this in two steps: 



After the extremalization with respect to /3 (cf. (IA.23aj) ). there remains to look for 
the minimal values (Et(^*), say) of the reduced energy function Et(z^) 

-2 



Et(*0 = • S p {v) (III.56a) 



Spiy) * £ ^f\ ■ ( IIL56b ) 

Such a minimalization process may be achieved by means of an appropriate com- 
puter program. The conventional predictions fjl.4j) would be obtained if the spectral 
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function S p (u) took its maximal values at 



(111.57) 



[n. 



i, 2, 3, 



Of course, this dream result cannot be attained by use of such a simple trial am- 
plitude as the selected $(r) flIII.50j) : but nevertheless the present approximation 
procedure yields predictions Et(z^*) which come amazingly close to the conventional 
predictions -Econv (0|, see the subsequent table 1. 



43 



n (— 1 4- 1} 

n p {— L p -i- i) 


/?( n ) r P vi (TrU) 


¥rr>(iJ- n h [pVl flTTT ^(\f\ \ 


H 

Is* 


BioL-E T (4 nl )ro/i 

F (n) K°J 


1 


/■> A i~\ A A 

-6.80290 . . . 


-7.23055 . . . 


A A A A A 

—0.2049 


-6.3 


2 


— 1. 70072 . . . 


-1.55087. . . 


1.7942 


8.8 


3 


-0.75588 . . . 


A A A A 1 A 

—0.66914 . . . 


3.7528 


11.5 


4 


a A r\ v -i o 

—0.42518 . . . 


-0.37366 . . . 


5.8740 


12.1 


5 


A A A 1 A\ 

—0.27212 . . . 


A c\ e~\ A A A 

-0.23906 . . . 


8.1307 


12.1 


6 


-0.18897. . . 


A -1 AAA t—T 

-0.16637. . . 


-| A h* A A A 

10.5044 


12.0 


10 


a a /"> a a o 

-0.06803 . . . 


A A A A A A 

-0.06069 . . . 


A A A V O A 

20.9538 


10.8 


15 


A AH AA 1 

—0.03024 . . . 


a a r\ r?o p* 

—0.02735 . . . 


35.7017 


9.5 


20 


—0.01701 . . . 


—0.01555 . . . 


51.9196 


8.6 


25 


a a -i aaa 

—0.01088 . . . 


A A 1 AAA 

—0.01003 . . . 


a a o r* o o 

69.3583 


7.9 


30 


A A A A 

—0.00756 . . . 


A A A rtA -I 

—0.00701 . . . 


O AAA ^7 

87.8627 


7.3 


35 


—0.00555 . . . 


A A A f-l n 

—0.00518 . . . 


107.3506 


6.8 


40 


A A A A C\ T* 

—0.00425 . . . 


A A A A A A 

—0.00398 . . . 


127.8215 


6.4 


45 


a a a o o r* 

—0.00336 . . . 


A A A A 1 C 

—0.00316 . . . 


149.3660 


6.1 


50 


—0.00272 . . . 


a a a c\ p* r* 

—0.00256 . . . 


172.1617 


5.8 


60 


f\ c\f\ i on 

—0.00189 . . . 


f\ f\f\ 1 *~7C\ 

—0.00179 . . . 


222.5187 


5.3 


70 


-0.00139... 


-0.00132... 


280.9499 


4.8 


80 


-0.00106... 


-0.00102... 


348.2861 


4.5 


90 


-0.00084... 


-0.00081 . . . 


422.8567 


4.1 


100 


-0.00068... 


-0.00065... 


501.4648 


3.8 



Table 1: Conventional Predictions and Spherically Symmetric 

Approximation 

The minimalization of the RST energy function E T (u) f lIH.56aj) - (IIH.56bj) with respect 
to the second variational parameter v occurs at the values (fourth column) and yields 
the corresponding RST prediction Et(za!> ) (third column). The relative deviation from 
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the conventional prediction EcoL (11.4)) is displayed in the last column; this deviation de- 
creases from (roughly) 10% at low quantum numbers n p < 10 up to 4% for high quantum 
numbers (n p ~ 100). Such a decreasement of the relative deviation is necessary in order 
that no cross-over with the conventional predictions E^mr can occur (see Fig. III. A). The 
negative deviation of the groundstate (n p = 1, first line) is presumably caused by the 
weak singularity of the trial amplitude $(r) [TTJ) (here = —0,204. . ., see Fig. 4b of 
ref. [H]). Such a singularity signals a certain crowding of the electric charge in the vicin- 
ity of the origin (r = 0) with a non-zero Dirac current k p on the z-axis, see both Fig.l 
and equation (VI. 30) in ref. [TTJ. Consequently, magnetism may not be neglected for this 
situation very close to the origin, and the presently used electrostatic approximation is 
assumed to break down for the groundstate (n p = 1) from this reason. On the other 
hand, the excited states (n p > 1) have vanishing Dirac current on the z-axis so that no 
magnetic excess force can occur and the electrostatic approximation allows for sufficiently 
realistic predictions. 

The average deviation (A[t], say) due to the presently considered lowest order of the 
spherically symmetric approximation may be defined through 



where the individual deviations A T (z/* ) are specified by the last column of table 1, i.e.: 



Thus the average deviation of the present predictions is found as Apr] = 7, 7%. This will 
become halved through regarding the anisotropic corrections (Sect.V). 




100 



(111.58) 




'conv 



(111.59) 
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Fig. III. A Energy Spectrum in the Spherically Symmetric 

Approximation 



The minimal values Et(z4 ) solid lines) of the energy function Et(V) (1III.56a|) are 
always unambiguously related to the corresponding conventional predictions Eml_ v fll-4j) 
(broken lines) because the relative deviations do not amount to more than (roughly) 10%, 
see also table 1 on p. SU Therefore there occurs no crossing-over of the conventional and 
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RST predictions which should provide a solid basis for estimating the magnitude of the 
anisotropic corrections in Sect. IV. The latter naive type of corrections turns out to be 
of the order of (roughly) 0.4% for the groundstate, see equations f1IV.45[) - (1IV.46|) below, 
but would deteriorate the present "isotropic" RST predictions by (roughly) 8% already 
for the first excited state, see App.C. These unrealistic large anisotropy corrections for 
the excited states necessitates then a reconsideration and improvement of the applied 
approximation method (Sect. V). 
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IV Anisotropic Corrections of Gauge Potential 



The necessity to go beyond the spherically symmetric approximation has already been 
stressed in the introduction (Sect. I). But in view of the preceding difficulties with obtain- 
ing exact solutions of the spherically symmetric type one will not expect that anisotropic 
exact solutions could ever be found. On the other hand, an estimate of the magnitude of 
the anisotropic effects would be highly desirable and we will try to attain this now at least 
partially, namely by expanding the formal solution (111.30}) for the electrostatic potential 
^Ao(r) with respect to the longitudinal variable $ of the spherical polar coordinates 
(r, (f)). This means that we relax the spherically symmetric approximation (IIII.ip to the 
following weaker condition: 

(p)Ao(f) => {p UoM) . (iv. l) 

Here the "^-dependent potential ^Ao(r, i?) needs not be the exact value of the integral in 
( 1II.30P but we will be satisfied with some approximative value in lowest order due to a 
certain expansion with respect to the anisotropy effect. 

In any case, the latter effect will not be fully taken into account since we still adopt 
that form of the wave amplitudes ^TZ±, &'S± as it results from the combination of the 
transformations fllll.13ap - flIII.13bp . f llll.14ap - flIH.14bp . f llll.17ap - flIH.17bp . f llTU^ - fliTl^O]) 
and flIH.23ap - (lIH.23bp . i.e. in the last end: 



(P)^ ±(r1 ^ <M0 •/*(*) (IV 2a) 

V r sin $ 

«S ± (f) =► * ±( , r) • fM . (IV.2b) 

V r sin $ 



Observing here the solutions f lIH.19P -f lIH.20p for the angular functions /r($) and fs($) 



one finds for the Dirac density ^ko(f) fill. 191) 



<f) 2 (r) -L Cf) 2 (r) 

(p)A; (f) =► { p ^ (r, d) = + f ~ { ] . (IV.3) 

47rr sm v 

If this result is substituted into equation (III.30p . the electrostatic potential ^A (r) reap- 
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pears as { p ^4o(r, $) in the following form 

^M^/^HMiHM, (IV . 4) 

J Anr'smv' \\r — r'\\ 
or in the non-relativist ic approximation ($ + (r) =>■ $(r), $_(r) =>• 0) 

*yMr,*) = ±[ ff „||! V i,| r (IV.5) 
47r y rsmu ||r — r || 

From here it is seen that the anisotropy of the wave amplitudes themselves is taken 
into account only to that extent as it is assumed also in the spherically symmetric approx- 
imation for the potential H4o(r) (1III.37j) ! Thus, the main anisotropic effect is thought to 



be due to the formal modification of the three-volume element in equation (1IV.5j) 



d t 

d 3 r =^ = rdrddd(t) , (IV.6) 

r sinw 



where the denominator is a consequence of the former transformation pII.14aP -f lIII.14bp 



of the wave amplitudes. Accordingly, our step beyond the spherically symmetric approx- 
imation consists now in computing the potential ^Aq{t^) 

{p U (r,tf) = ^ j dr'dd'dcf)' ^ (IV.7) 

in an approximate way (i. e. in lowest order beyond the spherically symmetric approxi- 
mation (IHI.lj) ) and then using the result for calculating the value of the energy functional 
Em (jIII.4ip on the present trial configuration {$(r), ^A (r, $)}. The corresponding ex- 
tremal values of E[$] within this class of trial configurations do then determine again the 
approximative RST energy levels of the bound system which must be opposed to both the 
conventional predictions and to the spherically symmetric RST predictions of Sect. III. 
Anticipating the result of this comparison, we will find that these anisotropy corrections 
are lowering the groundstate energy (first line of table 1 on p. |44"P by 0.4% so that the 
anisotropy correction appears to be actually a negligeable pertubation of the spherical 
symmetry. However, this pleasant result is readily shown to hold only for the groundstate 
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(App.B) and even here is spurious so that we are forced to develeop a more realistic 
pertubation procedure (Sect.V). 



1. Anisotropic Potentials 

The wanted method of expanding the electrostatic potential ^A (r, ■&) (IIV.7I) with 
respect to the magnitude of anisotropy will (in the simplest case) consist in naively ex- 
panding the denominator as follows: 

1 1 rr' 3 r 2 r' 2 



(r • r ') H - (r • r ') 



/\2 



r' 2 ) 2 



( r 2 + r /2)2 ( r 2 + r /2)f v 2( r 2 + 

+ ^^^ r (r^')^^— ^ r (f • f ') 4 + ... (IV.8) 

Z _|_ j,l2^ 2 o (j~2 _|_ ^>/2^ 2 



(f = 7^77 = - ||f|| 2 = ||f'f = l) 



Obviously, such a naive expansion would induce a corresponding expansion of the potential 
{p^Aq^^) of the following product type: 

{p} A)M) = {p U,(tf) • H4,(r) + {p U M (tf) ■ {p U„(r) + H4 m (tf) ■ {p l4 M |(r) 

+ {p U,v(^) • {p} A,v(r) + H4 v (tf) • {p} A/(r) + ... (IV.9) 

Here, the angular factors are found to look as follows 

H4„(tf) = / (r-r') = (IV.lOb) 

{rt Ai M = ~ / ^J£{?.?f = 1tt(1 + cos 2 0) (IV.lOc) 

Z J 47T lb 

H4 |v(l? ) = - / (r • r') 3 = (IV.lOd) 

2 J 47T 

{pU v (0) = — / ^^(r-r') 4 = — 7r{cos 4 + cos 2 0sin 2 ^ + -sin 4 0} (IV.lOe) 
v ; 8 7 4tt ; 128 1 8 1 V ; 
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Evidently, these angular factors do not display the expected tendency of decreasing 
with increasing order of the expansion terms; therefore such a decreasing tendency must 
arise now with the radial factors, otherwise the series expansion (1IV.9P of the potential 
{p^ofV^) would not be convergent. 

The radial factors emerging in the expansion flIV.9j) are found to look as follows: 



H4,(r 
{p U M (r 
H4 IM (r 
{p U, v (r 
{p U v (r 



a s ■ r 



a K ■ r 



a K ■ r 



dr' 



[r 2 + r 



>2\ 



,12 



dr 1 



4> 2 (r') 



(r 2 + r 



I2\ 



oo ,3 . $ 2^ 

dr 



o 



o 



/^.2 _j_ ^./2^ 2 

, r /4 • <j> 2 (r') 



oo r /5 . 

dr' 







^2 _j_ ^./2^ 2 



(IV.lla) 
(IV.llb) 
(IV.llc) 
(IV. lid) 
(IV.lle) 



Observe here that both potentials H4 M (r) (llV.llbp and {p U, v (r) f HV.il dp are irrelevant 
because they become eliminated by their vanishing angular pre- factors {p U M (i?) fHV.lObp 
and ^A\\/{p) (HV.lOdp . Concerning the relative magnitude of the remaining potentials 
H4,(r) f llV.llap . {p U m (r) fllV.llcp and {p U v (r) fllV.lle|> one could of course explicitly 
compute these potentials by resorting to the former trial amplitude $(r) f lIH.50p which 
then would immediately yield the desired relative magnitudes. However, one may cir- 
cumvent those unwieldy integrals by establishing a certain differential link between these 
radial auxiliary potentials pV.lla[) - (IIV. lle|) and then adopt a plausible guess for the 
lowest-order potential ^A\(r). This guess can then be taken as the point of departure 
for calculating the higher-order potentials ^A\\\(r), ' p Uv(r), ... just by means of that 
established differential link (see below). 
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Concerning now the relative magnitude of the partial potentials (r, i?) = ^A\ (i?) • 
u>}A\(r) etc., cf. ( 1IV.9|) . it is obviously sufficient to consider their radial factors ^A\(r) 
etc., cf. pV.llap - fllV.llep . Here, the first one (i.e. H4,(r) (llV.llap ) is singled out by 
its boundary values for r = 0, oo. At the origin (r = 0) one finds 



{p U,(r) 



r=0 



n s / dr'$ 2 (r') 
'o 



(IV.12) 



i.e. for the spherically symmetric trial function $(r) flHI.SOp 

1 



I V J 



7T 



r=0 7T ' ' 2U + 1 

Combinig this with its angular associate ^A\ ($) (lIV.lQaP yields then at the origin 

2/3a s 



(IV.13) 



{p U,(r,tf) = [ {p U,(tf) ■ {p U,(r)] 



which just agrees with the value of the spherically symmetric approximation H4 (r) at 
the origin (r = 0), cf. the equation flA.101) of App. A. On the other hand, this value 
flIV.12j) of the first auxiliary potential ^A\(r, $) (or its specialization flIV.14j) . resp.) at 
the origin (r = 0) actually coincides with the corresponding value of the "exact" solution 
H4 (r,tf) fHV7|) : 

{p U (r,tf) = — / dr'dtfdfi <§ 2 (r') = - a s f dr'$ 2 (r') = - ■ {p U,(r = 0) , 
r=o Ait J 2 ,/ 2 

(IV.15) 

see equation (IIV.12p . Thus the conclusion is now the following: since already the first 
auxiliary potential { p U,(r,tf) = H4[(f?) ■ H4,( r) in the expansion fHV^P for H4 (r,tf) 
adopts the correct value of ^Aq(t, $) at the origin, all other auxiliary potentials must 
necessarily vanish at the origin (r = 0); or expressed in terms of the radial potentials 
(llV.llbp -(UV.llel) 



(IV. 14) 



{?} A 



in ir, 



r=0 



{p U v (r) 



. 



r=0 



(IV. 16) 



A similar conclusion does also hold for the behaviour of the radial auxiliary potentials 
at infinity (r — > oo). First, observe here that the "exact" anisotropic solution (IIV.7P 
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adopts the Coulomb form for r — > oo 



poo 

lim {p U (r, 0) = — ■ - / dr' r'¥(r') = — , (IV.17) 
r->oo r 2 n r 







provided the non-relativist ic wave amplitude $(r) obeys the normalization condition 
(1III.38I) . But on the other hand, the first-order radial potential ^A\(r) (liV.llaj) has 
just the same limit form (apart from the angular factor (lIV.lQap ): 

lim {p l4,(r) = -• — , (IV. 18) 

r-¥oo 7T T 

and therefore the expansion (1IV.9I) tells us that the higher-order potentials must vanish 
faster than - at infinity: 

lim {p U m (r) = lim {p U v (r) = . . . = o(\) . (IV.19) 

? — >oo ? — >oo T 

Since both the second potential ^A\\[r) (lIV.llb|) and the fourth potential ^A\y(r) 
( llV.lldj) become annihilated bv their angular counterparts ^ p ^An (lIV.lQbj) and ^Aw/i^fi) 



(lIV.lQdj) . it is sufficient to look for the announced link only between the non-zero poten- 
tials H4,(r) pv.llap . {p U m (r) fllV.llcjt and ^A w (r) pV.llej) . For this purpose, one 
differentiates twice the first potential ^A\(r) fllV.llaj) and thereby finds immediately for 
the third potential {p U m (r) fllV.llcj) 

H4,„(r) = - l - r ^- 2 [v H4,(r)] . (IV.20) 

A similar procedure yields for the fifth potential ^Ay{r) fliV.llej) 



{P ^v(r) = -4r4^(- 3 - {P HiiM) 



35 dr 



(IV.21) 



etc. The importance of this procedure is now that one can guess some reasonable approx- 
imation for the lowest-order potential ^A\{r) which is manageable more conveniently 
than the integral on the right-hand side of equation fllV.llaj) but which on the other hand 



comes sufficiently close to the value of that integral. If this guess for ^A\{r) is substi- 
tuted in equation flIV.20p . one obtains the corresponding approximative expressions for 
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the higher-order potentials ^A\\\(r) and ^ p ^4v(V) whose relative magnitude can then be 
estimated. 



2. Simplified Groundstate Potentials (v = 0) 

For the purpose of a simple demonstration of our approximation method, a prelimi- 
nary choice for the lowest-order potential ^A\(r) could perhaps consist in the spherically 
symmetric potential ' p l4o(r) ( (IA.2j) plus (IA.5j) for v = 0) which has been treated exten- 
sively in the precedent Sect. Ill (and also in App. A, where the trial amplitude $(r) 
( 1III.50P does underlie). Defering again the purely technical subtleties to App. C, we can 
be satisfied here for the sake of demonstration with the special case v = where the 
selected trial amplitude $(r) ( 1III.50j) becomes especially simple because there remains 



only one variational parameter 

$(r) $*e^ r (IV.22a) 
$, = J- (2/3) . (IV.22b) 

V 7T 

The corresponding spherically symmetric solution [ p i4 (r) of the Poisson equation 
( 1III.36I) looks as follows: 

bU (r) = — (1 - e~ 2/3r ) . (IV.23) 
r 

Thus the asymptotic Coulomb form (1IV.18I) of the first auxiliary potential ^ P ^| (r, $) is 
safely adopted by the present spherically symmetric approximation (1IV.23j) . Furthermore, 
the value of the "exact" anisotropic solution ^A (r, $) (1IV.7[) at the origin 

{p} A (r,$) =a s - dr'$ 2 (r') = 2(3a s (IV.24) 
r=0 2 J 

also does agree with the corresponding value of that spherically symmetric approximation 
H4o(r) flIV.231) . Summarizing, the boundary values of the "exact" anisotropic potential 
^ p ^ ( r )^) P"V-7p are truly reproduced by the present spherically symmetric candidate 
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' p l4o(r) (1IV.23I) ; and this justifies to preliminarily adopt just that potential ' p L4o(r) as the 
starting point for generating the series of higher-order approximations via the differential 
link (jIV.20p - fllV.2ip . i.e. we try 

{p>A,(r) =► - • [p U (r) = -— (1 - e- 2/3r ) . (IV.25) 

7T 7T r 

But once such a convention concerning the starting potential ^A\(r) has been at- 
tained, one can now substitute this into the differential links (IIV.20p - flIV.2ip in order to 
calculate the higher-order potentials ^Um(r) and ^A\/(r), etc. Indeed, by this simple 
process of differentiation one easily finds the following results 

H4 m (r) A a s (3 2 r e" 2/3r (IV.26a) 
frU v (r) =► --^ (2/3a s ) e"^{^ - \ (2/3r) 2 - \ (2 fir)} . (IV.26b) 



Fig. IV. A below presents a sketch of all three radial potentials ^A\(r) (IIV.25P 



^ p ^4in (r) (1IV.26aP and ^A\/(r) (|IV.26b[) which clearly shows the hierarchy of their magni- 



tudes. Thus it is evident that the expansion (1IV.9P of the anisotropic potential ^ p 'Ao(r, i?) 
yields a (more or less rapidly) converging series of radial auxiliary potentials OlV.llap etc. 
which constitutes a kind of "expansion with respect to the magnitude of anisotropy" . For 
general variational parameter u, this effect is studied in App.B. The result is that for 
all values of the variational parameter v (> — |) the magnitude of the radial potentials 
{p}_A|(r), {p}A]\](r), ^A\/(r), •■■ is decreasing rapidly enough (see Fig. B.I) in order that 
one can tentatively restrict oneself to the first anisotropic contribution { p bl|[|(r, i?) for 
an estimate of the anisotropic effect upon the energy Et which is concentrated in those 
anisotropic field configurations. 
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Fig. IV. A Radial Auxiliary Potentials ^A\, ^Um, ^ P ^A V 

The potentials (lIV.25j) - (1IV.26bl) for v = are sketched in dimensionsless units, cf. 
equations (1B.6I) - (1B.10I) of App. B. For the present case v — there is a distinct hierarchy 
of magnitudes (^a < ^^a < '"So) which lets the anisotropy corrections for the 
simplified groundstate {y = 0) appear as small as only 0.4% (see the discussion below 
equation flIV.43j) ). For the excited states {u > 0), this hierarchy of magnitudes becomes 



flattened (see Fig. B.I on page 1 1 08[) and consequently the anisotropy corrections of the 
excited states do amount up to some unrealistic 10% and more, even for the first excited 
state (see the table on page II 14|) which necessitates to set up a more adequate series 
expansion than given by equations ( lIV.9j) - (1IV.llej) . see Sect.V. 
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3. Energy of the Simplified Anisotropic Groundstate [y = 0) 

After a choice of some trial amplitude $(r) (1IH.50|) has been made together with the 
subsequent determination of the corresponding anisotropic potential ^A (r, ■#) (1IV.7I) 
(to be used only up to a certain order of approximation) one can now tackle more 
concretely the problem of the total energy Et which is carried by such an anisotropic 
configuration. Here it may be sufficient (for a first rough estimate of the anisotropy 
effect) to restrict oneself to the inclusion of only the first anisotropic correction term 
{p U M |(r,$) = H4|m($) • {p U| M (r). By this restriction, one can represent the general 
approximation method in a rather transparent way without the numerical calculations 
becoming too tedious. But observe here that we have to admit a certain complication 
which consists in the fact that our adopted approximative potential (cf. (IIV.9j) ) 

{p} A {r,$) => <PU,(r,0)+ {p U| M (r,tf) = H4,(tf) . H4,(r) + H4 m (tf)- {p U m (r) (IV.27) 

is not an exact solution of the non-relativistic Poisson equation 

/ i $ 2 (r) 
AH4o(r,0 =-Ob , (IV.28) 

r smv 



i.e. the non-relativistic version of f lIII.16j) . As a consequence, the non-relativistic Poisson 
identity Nq = (1H.35|) will not be satisfied, in contrast to the normalization condition 
N$ = (1111.461) which merely fixes the normalization constant (see equation flIII.5ip ). 
Therefore the energy functional E[$] (IIII.41I) does not reduce to its proper physical part 



e2Y (= 2i?kin + Eft ) but must be kept somewhat more general: 

E [s] 2E kin + E R e} + \ { ° ] ■ N§ } . (IV.29) 

Thus the problem is now to look for the value of this functional on the class of trial 
configurations for the groundstate which by assumption is parametrized by only one vari- 
ational parameter (3 (because of putting v — 0). This will yield a function E$(/3) whose 
minimal value is just the desired total energy (E^, say) of the anisotropic groundstate; 
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and the crucial questions with this result are now (i) how close is this RST prediction 
E^ to the corresponding conventional prediction E^nv (11-41) . and (ii) how large is the 
anisotropic contribution (beyond the spherically symmetric approximation) included in 
this RST prediction ? 

For working through this program, one evokes the trial amplitude $(r) for the ground- 
state flIV.22aj) - flIV.22bj) in combination with the associated (but approximative) gauge 
potential {p U (r,$) f)IV.27|) . with the angular parts {p U,($) and {p U m (f9) given by 
equations (lIV.lQal) and (lIV.lQcl) and the radial parts ^A\(r) and ^A\\\(r) being given by 
(1IV.25j) and flIV.26aj) ; i. e. the (one-parametric) trial potential due to the (one-parametric) 
trial amplitude $(r) flIV.22ap reads ultimately 

{ p U Q (r,tf) =► ^.^ (1 _ e -2^ )+ 3 co§2 8_ as/3Ve _ 2 ^ 

2 tc r lb o7T 

Both these trial ansatze $(r) fllV.22aj) and {p U (r, i?) fllV.301) are now to be used in order 
to calculate any individual contribution to the energy function ]£$(/?) as the value of the 
energy functional E[$] (1IV.29I) on this trial configuration. 

The first contribution is the kinetic energy E^ n (IIII.42I) which itself is the sum of the 
radial part ^'E^ m (1III.43|) and the longitudinal part ^E kiQ (lHI.44p . The latter, however, 
must vanish for the groundstate because here we assume the quantum number l p of angular 
momentum flIII.21j) to be zero (=>• l p — 0). Thus the (one-particle) kinetic energy is of 
purely radial nature, cf. flIII.52p - flIII.53p for l p = and v = 

E kin ^^E kin = ^^ I drr[^} =^-(Pa B f (IV.31) 




Me 2 

The next energy contribution is the electrostatic gauge field energy £r } fHT36l) which 
reads in terms of our anisotropic potential ^Ao(r, (1IV-30I) 

4 e U4 c} = -^ J d s r\\^ {p} Mr^)\\ 2 = E} c} + E^ + E^ , (IV.32) 
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with the three contributions being found by explicit integration as (cf. App. C) 



E} e} = -hca s (3 = -e 2 /3 (IV.33a) 

£,l e} = ~ hca s (3 = -\ e 2 (3 (IV.33b) 
o 6 

E}° } = — — hca s (5 = — — e 2 (3 . (IV. 33c) 

111 160 160 v ' 

Here, the first contribution (IIV.33aj) is due to the spherically symmetric constituent 
h>U|(r, 0) (= {p U,(0) ■ t p U|(r)) of the potential fHVT9|) with the angular part h>U|(tf) be- 
ing given by (1IV. 10a|) and the radial part ^A\{r) by the postulate (II V. 251) . Therefore the 
first expansion term (r, •§) just coincides with the spherically symmetric approxima- 
tion H4 (r). Consequently, the present first contribution e} ( 1IV.33a|) must be identical 
to the field energy E^ ( 1III.54P of the spherically symmetric approximation, i. e. we must 
find for the sake of consistency: 

E } e} = -e 2 (3 = -2(3e 2 • £pot (0) = E§ . (IV.34) 

However, this consistency requirement 

2-e pot (0) = l (IV.35) 

is easily validated, namely simply by specializing the general result (IA.28|) of App. A 
down to v — which yields 

^ot(O) = 1 = £ ■ (IV.36) 

n=0 

But whilst the first contribution E\ e ^ (IIV.33aj) does refer to the spherically symmetric 
approximation, the other two contributions E\^ and (1IV.33bj) - (1IV.33cj) must repre- 
sent the anisotropic correction (E^ , say) 

£?W = E}. c} + E}.? = -(- + — ) e 2 (3 = -— e 2 (3 . (IV.37) 
v 6 160 ; 480 K ' 

However, the point with this anisotropic correction E^ flIV.37j) is now that it obviously 

cannot be considered as a negligibly small correction of the spherically symmetric term 
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E} e} fH\T34j) ! Thus the problem arises now how it may come about that the anisotropy 
corrections can be effectively smaller than the spherically symmetric contribution E^ 
flIV.34j) ! Otherwise one could in general not consider as the leading term of a mean- 
ingful perturbation expansion. Indeed, we will readily show that this effect is based upon 
a certain peculiarity of the principle of minimal energy. 

This question concerns now the last term of the energy functional Er$i on the right- 
hand side of equation f lIV.29j) . Namely, this contribution is non-zero whenever the Poisson 
equation is not satisfied by the selected trial configuration, i. e. more precisely: whenever 
the trial configuration { ^A (r, $(r)} does not obey the Poisson equation (IIV. 28 P for 
the gauge potential ^A (r, Since our selected approximation flIV.27j) for { p }Ao(r, i?) 
contains only the first two non- vanishing terms of the general expansion flIV.9j) . it is 
not linked to the chosen trial amplitude $(r) (IIV. 22ap - (II V. 22b]) via the Poisson equation 



( 1IV.28|) ; and therefore the anisotropic Poisson constraint (1II.35j) is violated. This then 
necessitates to explicitly calculate the constraint term Nq 

Ng } => = - M^c 2 (IV.38) 

which enters the energy functional E[$] (1IV.29|) as the last term on the right-hand side. 

Here, the gauge field energy E^ has already been determined by equations ( 1IV.32|) - 
(IIV.33cp so that we are left with the calculation of the (non-relativistic) mass equivalent 
M^c 2 (jTOTp . If our approximative potential ^A^r,!!)) (IIV.27[) . i.e. more precisely 
(IIV. 30 P for the present groundstate situation, is substituted therein, one finds that the 
mass equivalent appears as a sum of two terms 

M (c) c 2 M {c} c 2 = Mf o} c 2 + Mf,, } c 2 , (IV.39) 



where the first term contains the spherically symmetric approximation ^A\(r) (IIV.25P 
of the potential ^ p ^o( r ; $) an d yields the following result 

. . 2 /*°° 

Mf c} c 2 = - (-J he J drr {p U,(r) ■ <l 2 (r) = -hca s (3 = E{ e} , (IV.40) 
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cf. (1IV.34j) ; and furthermore the second contribution is due to the first anisotropic cor- 
rection flIV.26ajl of H4 (r,tf): 

2 



if,;>c 2 



(o_ \ ^ poo o 
— J he drr {p U| M (r) • <l 2 (r) = hca s (3. (IV.41) 
8 / Jo 32 



(For both results flIV.401) and (IIV.41I) the angular parts of the potentials are already 
integrated over). It should not come as a surprise that the lowest-order mass equivalent 
M^c 2 flIV.40p exactly agrees with the first contribution E^ (IIV.33aP to the gauge field 
energy E^ } fjnT32l ! The reason for this is of course that both quantities E^ (1IV.33al) 
and Mf e} c 2 fllV^Ol are based upon the spherically symmetric approximation [ p L4o(r) 
for the potential ^A {r^) an d also upon the same trial amplitude $(r) (IIII.50P with 
v = 0, so that the spherically symmetric Poisson constraint is satisfied in first order (i.e. 
E} e} -M.} e} c 2 = o). But the consequence of this peculiarity is now that the non-relativistic 
constraint term (IIV.38[) is built up exclusively by the anisotropic corrections: 

N§ } = (E} c} + E}° } + E}° } ) - (Mf 0} c 2 + M^ } c 2 ) (IV.42) 

{e} ^{e} _ Tfirfe}^ = ]]_ 
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E^ + E^-M^c 2 = -^hca s p. 



Collecting now all the partial results, i. e. the kinetic energy E kin (1IV.31I) plus the 
gauge field energy (lIV.33ap - (lIV.33cP plus the Poisson constraint term (II V. 42 p . and then 
substituting this back into the energy functional E[$] (IIV.29I) yields the value of this func- 
tional on the chosen trial configuration {$(r), ^Aq(t, $)} in form of an energy function 
(E$(/3), say) which is found to be of the following shape: 

E*(/3) = £-{(/3a B ) 2 -gl(/3a B )| . (IV.43) 

(Observe here also the value flIII.49bp for the Lagrangean parameter Xq*). It is now 
this result (IIV.43P which is well-suited for a first estimate of the influence of anisotropy. 
Namely, the spherically symmetric approximation does also produce an energy function 
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(i.e. EJ- w ^(/3, v), cf. (1A.22j) ). which for the special case v = looks as follows: 

2 2 

EP v ](/3, 0) = — {(2/3a B f ■ £ ki „(0) - (2/3a B ) • £ po t(0)} = — {(/3a B ) 2 - (/3a B )} 
a B a B 

(IV.44) 

where the special values of the functions £kin(^) and £ P ot(^) may immediately be read off 
from the equations (1A.27I) (=>- £kin(0) = |) and from f)IV.36|) (=>- £ po t(0) = |). Comparing 
now the anisotropic result (IIV.43|) to its present isotropic counterpart f)IV.44|) demon- 
strates clearly that the anisotropic corrections for the groundstate will amount to roughly 
4 parts in 10 3 ! More concretely, one determines the groundstate energy as the minimal 
value of these energy functions f lIV.43j) and (1IV.44j) and thus finds for the "anisotropic 
groundstate" ( 11 V. 430 



E* 1} = EJf3) = -4— ■ ( —) - -4— ■ 1-00417... (IV.45) 
1 in 4a B 1 480/ 4a B v ; 



and for the "isotropic groundstate" (1IV.44|) 



e 2 



E [IV1 (A0) = -J_ = £« v (WAG) 

min t:0 B 



e 2 



6.8029... eV) . 



4a B 

From these results it is now obvious that the influence of the anisotropy on the ground- 
state energy amounts to only 0.4%. It is true, our trial amplitude $(r) (IIV.22aP - flIV.22bj) 



with the fixation v = is not the optimal choice for the groundstate; if one admits non- 
zero v and carries through the isotropic (!) minimalization procedure, one will find the 
minimal value of the corresponding energy function at —7.23... [eV] for v = —0.2049... 
(see table 1 on p. HUand fig.s 4a-4b of ref. [H]). This anisotropic lowering (1IV.45|) of 
the groundstate energy below that simple isotropic RST prediction (IIV.46|) (which despite 
its naive deduction equals exactly the conventional prediction Econv (II.4p ) must now be 
tested for the whole positronium spectrum. But since the excited states require here the 
use of (at least) a two-parameter trial amplitude, we have to admit the more realistic 
trial amplitudes <£>(r) (IIII .50|) with v ^ in connection with an improved (!) expansion 
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with respect to the magnitude of anisotropy (see Sect. V). Namely, the point here is 
that the present anisotropic groundstate lowering (1IV.45j) is actually false, as pleasant as 
its smallmess (~ 0, 4%) may appear in the light of a pertubation expansion. Indeed, we 
will readily demonstrate now that an improved expansion with respect to the magnitude 
of anisotropy, not so naive as the above proposal flIV.9j) . will for the excited states yield 
anisotropic corrections in the range of some 6% as opposed to the present 0,4% due to 
our naive expansion (IIV.9|) . 
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V Anisotropic Corrections for the Excited States 



Since, for the treatment of the excited states (App.C), the expansion flIV.8p -f lTV.9p 
of the anisotropic gauge potential ^A (r, •&) flIV.5p in combination with the lowest-order 
postulate ( 1IV.25P has been revealed as being too naive, we have to develop now a more 
systematic way of expanding the gauge potential around that spherically symmetric ap- 
proximation H4 (r) (1IV.23P - The crucial point here aims just at the very concept of 
the spherically symmetric approximation, being termed as H4o(r). Naturally, one would 
associate to this approximation H4o(r) the property that it represents the average of the 
original potential ^A (r,i!}) on the 2-sphere, see equation (IV. 8 P below; or rephrased in 
other words: the average of the anisotropic constituent ( ^A an (r, d)) of ^A (r, d) on the 
2-sphere must vanish, see equation ()V.7j) below. But these requirements are not satisfied 
by our naive expansion (1IV.8P -f lIV.llep . Subsequently, we will iron out this preliminary 
deficiency of our anisotropy expansion and thus will obtain distinctly improved predictions 
for the positronium spectrum (cf. table 2 below, p. [771 to the former table 1, p. l4*4"p . 



1. Expansion around Spherically Symmetric Approximation 

To this end, we recall the integral representation ( 1III.37|) of the spherically symmetric 



approximation ^Ao(r) (for general v) and use this for the desired expansion of {p U (r,f?) 
around ' p L4 (r), namely by putting 

{p UoM) = [p U (r) + ( {p UqM) - [p U (r)) (V.l) 
a R f d 3 f' <f> 2 (r') a s f d 3 r' ( 1 n\ $ 2 fr n 





f d 3 f' 


( 1 


-a 


4tt J 




V sin i?' 





r \\r — r \\ 47T J r \smv' IJ \\r — r \\ 
Here, the second part is obviously to be considered as the truly anisotropic constituent 
H4 an (r,$), i.e. we define 

H4 an (r,tf) = {p U (r,i?) - [p U (r) (V.2) 
a s fd 3 f' f 1 tt\ $V) 
Atx J r' \sind' 2 J \\r — r'\\ 
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This anisotropic constituent of the gauge potential may now be expanded by use of equa- 
tion (1IV.8j) which yields a series expansion quite similarly as shown by equation (1IV.9)) 



{p U an (r,0) = {p U'(tf) • {p} A\r) + {p U"(tf) ■ {p} A u (r) + {p U'"(0) ■ {p U'"(r) 

+ {p U lv (0) • {p U lv (r) + {p U v (0) • {p U v (r) + ... (V.3) 



Although this looks very similar to the former expansion flIV.9j) . the individual terms 



are somewhat different especially concerning the angular pre-factors: 

™A*) = f—(^r--)=0 (V.4a) 

y ' J An \smtf> 2 J K J 

H4"(tf) = f — ( -V ~ t) (?' = ( v - 4b ) 



An 2 r 

{pU iv (^) = I I ^ ( JL - ^ (? • F') 3 = (V.4d) 



2 / 4tt Vsintf' 2 



{pU v (tf ) = ^ ( J— - -) (f . f ') 4 = — vr{cos 4 + cos 2 sin 2 + - sin 4 - — } 
V ; 8 7 4vr Vsin^' 2y v ; 128 1 8 15 / 

(V.4e) 



These angular factors { p U N (0) (N = I, II, III, IV, V, ...) have zero average over the 2- 
sphere, i. e. 

f {p} A n/^ = J_ sin ^ {pU n (0) = , (V.5) 

J An An J 

because of the commutativity of both angular integrations, i. e. 



J An K ' J An 2 J J An 



( 1 n\ 1 . . 

Therefore the average of the anisotropic potential ^A an (r, 0) ( IV. 3j) is zero: 

^ {pU an (r, 7?) = ; (V.7) 
An 
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and thus the spherically symmetric approximation H4o(r) (1III.37j) . as the leading term 
of the series expansion (1V.lj) - (1V.3|) . is revealed as the angular average of the general 
anisotropic potential ^A (r, (1IV.7j) . i.e. 



tPU (r) = J ^ WXo(r^) . (V.f 



Concerning the radial factors of the series expansion (IV.3j) . one observes on account 
of the vanishing of the first two angular factors ' p bl ($) (IV.4aj) and { p >A ($) flV.4b|) that 
the first active non-trivial potential is that of third order 

H4'"(r) = a s f dr'r'^—^P-= {p U|„(r) , (V.9) 
J (r 2 + r' 2 )§ 

cf. fliV.llcp . and the next non-trivial one to be considered here is of fifth order (cf. 

{p} A v (r) = a s [ dr'r' ^ ' = ^Ay{r) , (V.10) 

J (r 2 + r u ) 2 

cf. (liV.llej) . Thus the present expansion (IV. l|) -f TV73|) of the gauge potential is essentially 
the same as the naive one (IIV.9p . and differs from that former one only in the angular parts. 
Consequently, the hierarchical arrangement of magnitudes, as expressed in Fig. IV. A, 
should in principle apply also to the present expansion (1V.lj) -( TV\3|) . 



2. Energy Content of the Anisotropic Configurations 

The present expansion around the spherically symmetric approximation has also some 
useful consequences with respect to the anisotropic energy functional E[$] (1IV.29|) . The 



crucial point here is the circumstance that with the specification of the angular pre- 
factors t p t4 ($), { p }A u ($), ... one can integrate over these angular functions so that the 
principle of minimal energy becomes reduced to a purely radial problem for the remaining 
spherically symmetric fields H4o(r), ^A (r), $(r). This reduction to a spherically 
symmetric problem is carried through as follows: 
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The first consequence of the representation of the anisotropic potential ^ p ^o( r , as 
a sum of two terms, cf. (IV. lj) 



{p l4 (r,tf) = [p U (r) + {p U an (r,tf) , (V.ll) 
refers to its energy content flrOHjl 



4 6) 4 G} = / d'r ||V frU (r,^)|| 2 (V.12) 
which appears now as a sum of three terms 

4 } = 4 + ^ } + 4 e x } - (v.i3) 

Here, the first term is the energy content of the spherically symmetric approximation 

[p U (r) 



4 = -^-/rfvuv^own 2 



a s J \ dr 
and is for our trial amplitude $(r) 5Q|) explicitly given by 



' dr r (V.14) 



4 = -— (2a B |S)-e P otM (V.15) 



a B 

where the potential function e pot (u) is specified by equation (1A.28|) of App. A, see here 
also equation (IA.22I) and App. D of ref. [12] . Analogously, the second energy contribution 
Ei^ is the energy content of the anisotropic part ^ p ^y4 an (r, $) of the gauge potential and 
is given by 

£|n } = J d 3 r ||VfrU an (r,^)|| 2 . (V.16) 

This result says that the anisotropy energy is always negative and therefore must increase 
the binding energy! 
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And finally, there is a mixed contribution £^ which is built up cooperatively by both 
parts of the gauge potential 

E S = —^- J d 3 f (v [p] A (r))'(f {p} A an (r,^ . (V.17) 

But the interesting point with this last contribution is now that it actually vanishes 

E^=0. (V.18) 

This can easily be verified by simply observing the fact that the spherically symmetric 
approximation l p L4 (r) does not depend upon the polar angles $ and <p so that 

VH4 oW = ^Me,, (V.19) 
dr 

and consequently the mixed term depends exclusively on the derivatives in radial direction 
(~ e r ), i.e. we have by use of the anisotropy expansion (IV. 3j) 



mx ~ 2na s J dr dr lV Uj 



2na, 



o 



CO 



Jo dr dr J 

Jo " r dr J 



Thus the vanishing mean value ( IV. 6 j) of the angular factors { p }A N ($) just validates the 
claim (IV.18p . and consequently the gauge field energy becomes composed exclusively 
by the two quadratic terms (IV. 14j) and (IV. 16j) 



4 e> = 4 C] + E£ ■ (V.21) 

The intended reduction of the original angular-dependent variational problem to a 
purely radial problem must here be carried through only for this gauge field contribution 
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(lV.2ip because the isotropic part is already defined as a purely radial integral, 
cf. (1III.45P ; and for the anisotropic part E^ (IV. 16|) one integrates over the angular 
pre-factor '(i?) (IV.4cl) so that a purely radial integral over ^ p ^'"(r) is left over, see 
equation (1V.30I) below. 

Clearly, if one could solve exactly the Poisson equation (1IV.28I) for the anisotropic 
potential ^A (r,i!)), then the Poisson constraint (IIV.38P would not be a point of concern 
because it would disappear trivially from the energy functional E[$] (1IV.29|) . But since 
we will not find such an exact solution for the anisotropic gauge potential, we are forced 
to briefly inspect the Poisson constraint with reference to that splitting (IV. lip of the 
gauge potential into an isotropic and anisotropic part. Indeed, this splitting induces a 
corresponding splitting of (the non-relativistic version of) the original mass equivalent 
M^c 2 (1IL37P into a spherically symmetric part M^c 2 and an anisotropic part Min c 2 : 

= -he [ d 3 f {p U (r,?9) - (V.22) 

J Anr sin v 

= M [c] c 2 + M^ } c 2 , 



!{e} c 2 



with the isotropic part being given by equation (1III.48P and analogously the anisotropic 
part by 

Pc 2 = -— / H4 an (r,tf) • $ 2 (r) . (V.23) 

An J r sin v 

Concerning here the intended reduction to a purely radial problem, one integrates again 
over the angular part and is then left with a radial integral, see equation (1V.3ip below. 



As a further consequence of the splitting (IV. lip , the Poisson constraint Nq^ (1IV.38|) 
must then also be split up into two terms 

N§ } = Ng + N^> , (V.24) 

with the isotropic part Nq being given by 

Ng = - M [c] c 2 (V.25) 
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and similarly the anisotropic part N^} by 

Nl c n } = - M^c 2 . (V.26) 

Now the interesting point with this splitting (1V.24[) is the following: Since both isotropic 
parts ( 1V.14P and M^c 2 (1111.480 are built up by the spherically symmetric approxi- 
mation H4o(r) and this potential is an exact solution of the (albeit approximate) Poisson 
equation (IIII.36j) . the spherically symmetric part Wq flV.25j) of the general Poisson con- 



straint Ng' flg^D) is zero 



M 

Ng 1 = ; (V.27) 

and consequently the general Poisson constraint will consist exclusively of the prop- 
erly anisotropic part N^}: 

N& } Ntf . (V.28) 

whose vanishing must be considered separately (see below). Actually, this is the general- 
ized reformulation of the former result (IIV.40|) for the simplified groundstate. 

Thus, summarizing all the precedent partial results, the anisotropic energy functional 
Er$i flIV.29j) can be recast to the following 5*0(3) symmetric form (inclusive all con- 



straints) which then neatly displays the separation into the spherically symmetric and 
the anisotropic gauge field degrees of freedom: 

E m => E w = 2E kin + [E§ + Ag } • Ng] + {E™ + A§ } • N£>} + 2A S • N$ . (V.29) 

This result will readily become the point of departure for the approximate calculation of 
the "anisotropic" binding energy of the two-body system. But let us first remark here 
that the present functional E{$} f!V.29j) is to be understood in the sense that all angular 
factors are already integrated over so that we are left with a variational problem for the 
spherically symmetric fields $(r), f p l4o(r), (r), ^A w (r), ... The Poisson constraint 

term N§ is re-included here because this gives us the freedom to try (if we wish) also other 
configurations {$(r), H4 (r)} which eventually do not obey the spherically symmetric 
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Poisson constraint (1V.27j) . The same does hold also with respect to the "anisotropic" 



Poisson constraint N^} (1V.26j) . 



Observe now that all contributions to the anisotropic energy E{$} (IV.29|) are well- 
defined by our precedent assumptions: 

(i) the kinetic energy E^ in is given in terms of the spherically symmetric wave amplitude 
$(r) by equations fllll.42p - fllll.44p 

(ii) the "spherically symmetric" gauge field energy reads in terms of the SO (3) 
symmetric part H4 (r) of the gauge potential as shown by equation flV.141) 



(iii) the "spherically symmetric" mass equivalent M^c 2 is displayed by equation (1111.481) 

(iv) the "anisotropic" gauge field energy Ein (IV. 16|) becomes by the restriction to its 
lowest-order contribution ^A u \r), cf. (1V.3h 



320 a s J IV dr 



(v) the "anisotropic" mass equivalent Mll^c 2 (lV".23j) is found to be of the following form 

7r 2 r°° 

MW c 2 = he drr$ 2 (r) ■ {p U'"(r) (V.31) 

64 J 

where we are satisfied again with the first anisotropic correction ^A u \r). 



3. Principle of Minimal Energy 

With the {anisotropic!) trial configuration being specified now through the three 
spherically symmetric (!) fields $(r), H4o(r), ^ p ^'"(r) one can look for the extremal 
equations due to the functional E{$} (1V.29P with respect to just these three fields. This 
may be done by prescribing the functional form of the three fields and equipping this form 
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with certain variational parameters {(3, z/^} and then determining the extremal points of 
the corresponding energy function over the parameter space (see the former equations 
(lIII40p ). An other method would consist in leaving the three fields unspecified for the 
time being and deducing for them a coupled system of differential equations from the 
principle of minimal energy flIII.391) . After having solved this system, one may substitute 
back the corresponding solutions into the energy functional E{$} (1V.29I) in order to obtain 
its associated extremal value. An alternative "mixed" method would consist in prescribing 
the functional form (up to some variational parameters) of one or two of the three fields 
and determining the residual unprescribed field from its differential equation due to the 
principle of minimal energy. 

In any case, it is very instructive to first deduce the system of extremal equations 
due to that functional E{$} (IV. 29}) . The extremalization with respect to the spherically 
symmetric approximation t p L4o(r) yields just again the former Poisson equation flIII.36j) . 
Therefore the former link (IHI.37I) of the spherically symmetric approximation H4 (r) to 
the wave amplitude $(r) still persists which entails the continued validity of the Poisson 
constraint (1V.27I) . But the converse is not true, i.e. the Schrodinger-like wave equation 
(1III.34P receives an aditional coupling to the anisotropic part ^A u \r) of the gauge po- 
tential ^A (r, HV-3|) . Indeed, the extremalization of E{$} (IV.29I) with respect to the 
wave amplitude $(r) generates the following modification of the original equation (1III.34)) : 

(V.32) 

Thus it becomes evident that the gauge field influences the spherically symmetric am- 
plitude $(r) not only via its average potential H4 (r), i.e. the spherically symmetric 
approximation (IV. 8p . but also via its first anisotropy correction ^A u \r). However, this 
influence occurs in a very plausible way; namely the first anisotropy correction ^A u \r) 
is simply added to the spherically symmetric approximation l p l4o(r) in order to build up 
a slightly modified (but still spherically symmetric) potential, to be considered as some 
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kind of effective potential! But observe also that the coupling of the amplitude $(r) to 
the first anisotropic correction ^A u \r) is (roughly) ten times (~ ^) weaker than its 
coupling to the average potential ' p L4 (r)! On the other hand, it will readily turn out 
that the anisotropic correction ^ p ^'"(r) couples three times stronger to the amplitude 
$(r) than does the average potential f p L4o(r), see equation flV.34j) below. Therefore it is 



reasonable to take into account the anisotropy corrections for the gauge potential but not 
yet for the wave amplitude. 

Finally, the extremal equation of E($} (1V.29j) with respect to the anisotropic correction 



MA"'(r) must be determined. Since the latter field is contained only in two terms of 
namely the "anisotropic" energy Ei^ (IV.30P and the "anisotropic" mass equivalent m|^c 2 
(IV.3ip . the variational proces (5\\\, say) with respect to ^ p M'"(r) becomes very simple 

8 m E m = =► Sm (j5&> + Ag } • n£>) = , (V.33) 

n{e} 



with the "anisotropic" energy E^n given by equation (IV. 30}) and the mass equivalent 



Min^c 2 by (IV.31j) . The result of this extremalization process is the following Poisson-like 
equation 

£ + " •W. (r) = _ B(w «V>. (v .34) 
ar r ar r z J r 

Evidently, this equation owns some properties which are worthwhile to be mentioned 
briefly. Firstly, the coupling of the anisotropy potential * p }A (r) to the spherically sym- 
metric wave amplitude $(r) is encanced by a factor of ^ ~ 3, cf. the corresponding 
Poisson equation (1III.36I) for the spherically symmetric approximation (i. e. the average 
potential) t p l4o(r). This hints at the circumstance that the deviation from the spherical 
symmetry is relatively large so that one (erroneously) tends to expect a correspondingly 
large amount of anisotropy energy E^ (IV.30I) in comparison to the isotropic energy E^ 
flV.14j) . Secondly, the last term ( {P}A 2 ^ ) on the left-hand side of equation flV.341) is 



necessary in order to ensure the former boundary conditions flIV.16j) and (IIV.19j) for the 



anisotropy corrections of the gauge potential. Indeed, in the asymptotic region of three- 
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space (r — y oo, $(r) — > 0) the potential correction '(r) has to obey the equation 

which imposes upon the solution ^A u \r) the condition to vanish faster (~ -4,) than the 
Coulomb potential (~ -), i.e. more precisely (see App. D) 

{ P } A »i (r) * I {2v + 3)(2l/ + 2 ) . ; (V.36) 

of course, this hints at the quadrupole character of this field (see App. D). 

Finally, there is a further interesting point in connection with that "quadrupole equa- 
tion" (1V.34j) . This concerns the fact that the "anisotropic" Poisson constraint (1V.26I) 
vanishes whenever we can use an exact solution of the Poisson-like equation (1V.34P as our 
trial potential ^A m (r), quite independently of the selection of a trial amplitude $(r) (see 
App. D). This is the same effect as does occur in connection with the "isotropic" Poisson 
constraint Nq (IV.25I) . Thus, provided that we are able to use exact solutions of both 
gauge field equations (IIII.36P and (IV. 341) . then we can omit the corresponding constraint 
terms from the energy functional E[$] (1V.29I) . In this case, this functional becomes cut 
down to the following simple form: 

E m => 2E kiQ + 4 C] + ££} = 2E kiQ + 4 e} . (V.37) 

Here it is also assumed that the selected trial amplitude is normalized to unity so that 
the normalization condition (1111.460 is satisfied. 

4- Energy Spectrum of the Excited States 

After all these preparations, the calculation of the energy due to the excited states 
[y > 0) becomes now a simple matter: since the kinetic energy is proportional to the 
square of the variational parameter ft, cf. (1III.52h - f1III.53[) : and since furthermore the total 



74 



gauge field energy ( )V.21j) is a linear function of (3 (see equation (!E.9aP of App. E 
the extremalization process (lIII.40p . or more concretely 

5E< IV >(/3,v) 



d(3 

<9E< IV >(/3,z/) 
du 

with 



(V.38a) 
(V.38b) 



E {iv} (/3) „) = —{ (2a B (5) 2 ■ eU") ~ (2a B 0) • 45 W} , (V.39) 



yields again the desired energy values Ep n \ namely as the minimal values of the reduced 
energy function E? 

e ? } ^mL' ( v - 4 °) 

Here, the reduced energy function Ef , is again defined through 



e 2 



^ } = -^-4 a} M (V.41) 
with the "anisotropic" spectral function sl^(v) being itself defined through 

{e} / %2 
P £kin(^) 

The (dimensionless) total potential function e^(u) is given by equation (1E.9b|) of App. E 
and is obviously nothing else than the "anisotropic" generalization of the former e po t(' / ) 
(IA.28p . Thus the present energy function E^ Iv ^(/3, v) (1V.39|) is the "anisotropic" general- 
ization of the "isotropic" E[ Iv ^(/3, v) (1III.55H so that it merely remains to determine again 
the extremal values of the spectral function S^f 5 {y) (IV".42j) for any principal quantum 
number n p (= l p + 1); for the numerical details see App. F. 

The table 2 (below) displays the RST predictions of the positronium energy E^ } (TOO]) 
in an arrangement quite similar to the results of the spherically symmetric approximation 
in table 1 (p l44l) . Here it becomes obvious that the deviations from the conventional 
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results (last column in tables 1 and 2) are considerably reduced through the inclusion of 
the anisotropy corrections. Indeed, defining the "anisotropic " deviations by 

A T („J»>) = Ecom {n) K (V.43) 

-C'conv 

(last column of table 2) in a quite analogous manner to their "isotropic" counterparts 
A T (z/I n] ) (Q) (last column of table 1), one clearly realizes that the "anisotropic" devia- 
tions (1V.43P in table 2 amount to less than half of the "isotropic" deviations (table 1). 



Observe here that the average A{ T } of the deviations (IV.43[) . being defined by 



nE|Ar(^ W )l , (V-44) 



does amount now to A{t} = 3, 1% in comparison to the corresponding former devia- 



tion A[t] = 7,7% for the spherically symmetric approximation (table 1). The ground- 
state (n p = 1) is disregarded here because the presently applied electrostatic approxima- 
tion is supposed to fail for the lowest energy eigenstate (see the discussion of this in the 
precedent paper [TT]). 

Clearly such a result suggests to consider the higher-order approximations which even- 
tually do shift the RST predictions even closer to their conventional counterparts, see 
Fig.V.A below. 
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Table 2: RST predictions Ep ^V.4-0[) including the Anisotropic 

Corrections 



The "anisotropic" RST predictions Ep (1V.40P deviate only half as much from the 
conventional predictions -Econv as is the case for the spherically symmetric approximation 
of lowest order (table 1). Indeed, the average deviation A{ T } flV.44j) is now A{ T } = 3, 1% 
in place of Am = 7, 7% for the spherically symmetric approximation. Observe also that 
for quantum number n p ~ 15 the deviation is minimal as predicted by an inspection of 
the "isotropic" and "anisotropic" deviations A{„} and Ai-tVl" - ') in Fig. F. II, App.F 
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Fig. V.A Energy Spectrum Ep n ^ ftV.40\) inclusive Anisotropy 



Corrections 

The inclusion of the anisotropy corrections diminishes the deviations of the RST pre- 
dictions from their conventional counterparts by (roughly) 50%, cf. Fig. III. A: This im- 
provement suggests that the regard of the anisotropy effects of higher approximation 
order will shift the RST predictions (solid lines) even closer towards their conventional 
counterparts (broken lines). ^ 



VI Improving the Spherically Symmetric Approxi- 
mation 



Concerning the essential point with the anisotropy corrections, it is instructive to first 
recall those results being collected in table 2 on p. [77J These results let appear the 
magnitude of the anisotropy corrections in the range from 3% (for the first excited state, 
n p = 2) up to 0.3% (for the excited state with principal quantum number n p = 15). Evi- 
dently, the consideration of the anisotropy corrections in Sect. V improves the "isotropic" 
RST predictions (table 1) by (roughly) 50%. Surely, the higher-order corrections of the 
anisotropic type can be expected to shift the RST predictions further towards their con- 
ventional counterparts (II. 4ft ; but for the moment we are rather interested in testing the 
extent to which further improvements of the predictions may emerge within the framework 
of the spherically symmetric approximation itself isotropic corrections). 

In this sense, we study now the following three-parameter generalization $&(r) of the 
former two-parameter ansatz $(r) flIII.50p - flIII.5ip : 



$ 6 (r) = $^(6 + b ir )e-" r . (VIA) 

Here, the additional variational parameters are the (real) constants bo,bi whose presence 
entails a certain generalization of the normalization constant flIII.5ip so that the general 



normalization condition flIII.381) becomes specified down to the following form 

p 2 + 4(> + 1) • p oPl + (2u + 3) (2v + 2) • p\ = 1 . (VI.2) 

Here, the parameters po,pi are merely some convenient modification of the original ansatz 
parameters &o>&i flVI-ip . i.e. 



^-\/r r (^ + 2 )(^TT- 6 o (VI.3a) 



Vi^\l\ ■r(2z, + 2)^ 3 .& 1 . (VI.3b) 
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Clearly, for p\ — > 0,po — > 1 one is led back to the former two-parameter ansatz f lULSOj) - 
(jIII.511) . Observe also that the present normalization condition (1VI.2j) is a slight gen- 
eralization of the former condition (IV. 41) of ref. [H] which emerges from the present 
condition (IVI.2[) through putting v to zero {y — > 0). 



1. Improved Energy Function 

This new trial amplitude $?,(r) (IVI-lj) must now be used in order to calculate anew the 
kinetic energy E kin (lIII.42j) - (lIII.44j) and the gauge field energy E^ (1III.45I) with its corre- 
sponding mass equivalent M^c 2 (1III.48j) . Both objects will then obey the Poisson identity 
(1III.47P provided the electrostatic potential H4 (r) is taken again as the exact solution of 
the spherically symmetric Poisson equation flIII.36j) . Thus, both constraints (IIH.46|) and 
(IIII.47I) will be automatically satisfied and therefore the non-relativistic energy functional 
Er$i flIII.4ip will consist exclusively of its physical contributions (i. e. kinetic energy and 
gauge field energy): 

E m =► E[Jf ] = 2E kin + 4 6] • (VI.4) 

Clearly, this is again the non-relativistic version of the former relativistic case 

Turning here first to the kinetic energy -Ekin, one recalls that this energy contribution 
is the sum of the radial part ^E kin and the longitudinal part ^E kin , cf. equation flIII.42j) . 
Both parts must now be calculated by use of the generalized trial amplitude $&(r) fIVI.ip . 
This yields first for the radial part ^-E^m the usual form: 

(r) £kin = ^~ (2/3a B ) 2 • {r) e kin (^Po,Pi) • (VI.5) 

Here, the kinetic function of the radial type ^kin^, RbPi) is defined in terms of the 
dimensionless form <& up (y) of the original trial amplitude $b(r) (IVI.ip through 



(r) £kin(^,Po,Pi) = Wo 1 , ^ / dyy 



T(2z/ + 2) J ~"" [dy® vp ( y \ 
®u P {y) = y"(po + Pi2/)e~2 , see App. G ) . 



2 



(VI.6) 
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The explicit calculation lets emerge this radial kinetic function in the following shape 

<" £kl >,P„, Pl ) = + ^ ■ P0P1 + V -±± ■ Pi , (VI.7) 

i.e. a quadratic polynomial with respect to the variational parameters Po,Pi- Clearly, 
when the parameter p\ tends to zero (p\ — > 0) the present radial function ^kin^, po,Pi) 
( 1VI.7|) collapses to the corresponding radial part of the result ( 1III.53P which is due to the 
former two-parameter trial amplitude $(r) (IIII.50p . 

A similar line of arguments may also be applied to the longitudinal kinetic energy 
^-Ekin (IIH.44I) which then adopts an analogous form 

WE kin = 7^(2/3a B )% ■ W e kin (",Po,Pi) (VI. 8) 

with the longitudinal type of kinetic energy function ^'£un{v, Po,Pi) being found as 

£kin(z/ ' ^ Pl) = 2,(2, + 1) + 2^TT + Pl • (VL9) 
Clearly this quadratic form is again the longitudinal counterpart of the radial object 
£kin(^> Po>Pi) (IVI-Tf) : and if the variational parameter pi tends to zero po — > 1) one 
ends up here with the longitudinal part of the former £kin(^) (HII.53I) . Thus summarizing 
the result for the kinetic energy E k \ n flIII.42j) as the sum of the radial and longitudinal 



parts, one ultimately arrives at 

e 2 

£kin = 7— (2a B /3) 2 -e kin (v,Po,Pi) (VI. 10) 

with the total kinetic function £kin(^P0)Pi) emerging here as the sum of the radial and 
longitudinal parts: 

£km(v,P0,Pl) = (r) £kin(^,P0,Pl) + W £kin(^,P0,Pl) " ^ 

As a consistency check one easily verifies that in the limit p\ — > this three-parameter 
result (j VI. lip collapses again to the former result £kin( z/ ) (IIII 53j) for the two-parameter 
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variational ansatz. For v — and l p = 0, the present kinetic energy Ey, n (1VI.10I) agrees 
with the former kinetic energy E kin (1IV.22aj) - (1IV.22bj) of ref. [14] . 

The gauge field energy as the second constituent of the energy functional E^ 
(1VI.4P is considerably more complicated than the kinetic energy E kin . Here, the Poisson 
identity is of help in order to find the simplest form of this energy contribution (see 
App.G for the details). Summarizing the results, one finds the desired gauge field energy 
E [ l ] (I111.45P to look as follows: 

E § = ~— (2a B /3)-£pot(^Po,Pi) • (VI.12) 
a B 

Here, the total potential function £ po t(^, Po, Pi) is obviously the generalization of the former 
potential function e pot {v) (IHI.54I) which for the simple trial amplitude $(r) (IHI.50P is ex- 
plicitly given by equation (IA.28|) of App. A. But clearly, the present generalized trial am- 
plitude $b(r) ( IVI.lj) must now entail a more complicated potential function e pot (^, po,Pi); 
and indeed, the latter function is found to be of the following form: 

£ P ot(f,Po,Pi) = v\ ■ £ P ot(f) + • E\(v) + Apfyj ■ e n (v) + Ap pl ■ e m (v) 

+ p\ ■ Eiy (is) , 

cf. equation (IG.16P of App. G. 

This is obviously a fourth-order polynomial with respect to the variational parameters 
Po,Pi (1VI.3al) - (1VI.3bp where the four additional coefficients ei, ...,£rv are ordinary func- 
tions of the second variational parameter v. Concerning their specific shape, one finds 
(see App. G) that there are certain interrelationships so that it is sufficient to specify 
in detail only three of them (e.g. £ pot , ej, en) and the remaining two (i.e. em and £rv) 
are then related to the first three in a simple way. More concretely, the first auxiliary 
potential function e\(v) is given quite generally in terms of the (dimensionless) potential 
a u (y) through 

1 f°° 



T(2z/ + 2 



1 f°° 

J dyy 2 » +l e-y-a uH (y) . (VI. 14) 
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Here, the potential a u (y) is the solution of the Poisson-like equation 

y^-e-y (VI.15) 



d 2 a v (y) t 2 da J y) ,„_, 



dy 2 y dy 

which essentially is a dimensionless form of the original Poisson equation (lIII.36j) . There- 
fore the solution can easily be specified as 



oo 

n 



„ to) = -r(2, + 2).jl-e-».£ f p- T ^-«j , (VI.16) 

cf. equation (IG.13aj) of App. G; and if this is used in order to explicitly calculate the 
auxiliary function £\{v) (1VI.14I) one ends up with 

2v + 2 \ 1 ^ n r(4z/ + 3 + n) 



1 t 

2v + l ] T(2v + 3) • 2 4 ^+ 3 ^ 2" V(2u + 2 + n 



= 1 ! y»T^ + 3 + n) 

Y(2v + 2) -2 4y + 3 ^2 n r(2i/ + 3 + n) v ; 

\ / n=0 ' 

Next, the second coefficient en{u) of £ po t( z/ , Po^Pi) i n HVI. 13[) turns out to be the sum 
of two contributions, i. e. 

eniy) = e' n (v) + (VI.18) 
where the first part is related to the well-known e P ot(^) through 

e' ll (v) = {2v + 2) 2 -6 pot \ u+h (VI.19) 

and the second part is given by 

(i/ + 1)(2i/ + 3) f 1 A n r(4z/ + 4 + n) 

eil ^ V ' 2v + l 1 r(2z/ + 4) • 2 4 ^+ 4 ^ 2 n ' r(2z/ + 2 + n) 

s f 1 n r(4z/ + 4 + n)l 

= 2 '- + 1 )( 1 - r(2, + 2).2^ g^- rU + 4 + „) } • < VL20 > 

Finally, as mentioned above, the last two coefficients Em(v) and eiy(i>) in (I VI. 131) may be 
expressed in terms of the precedent potential functions which are to be shifted as follows: 

e m {y) = (2i/ + 2f ■ eA, (VI.21a) 

f-r 2 

e iy (u) = (2u + 3) 2 • (2u + 2) 2 ■ s pot \ u+1 . (VI.21b) 
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Thus the total potential function e pot (^,po,Pi) (1VI. 13|) ultimately appears as a well- 
defined function of the variational parameters v and po,Pi', and if this is substituted back 
into the gauge field energy E$ (1VI.12I) we obtain this quantity, as desired, as a function 
of all four variational parameters {(3; u, Po,p±} 

e 2 

£r = (2a B 0) ■ {Po ■ £ P otO) + ^P S oPi • e\{v) + 4pgp? • e-aiy) + Ap Q p\ ■ e m (v) 

a B (VI.22) 

+ pt-ejy(v) } . 

Here, one is easily convinced that for pi — y this four-parameter energy of the gauge 
field actually collapses to the two-parameter result flIII.54p being due to the simpler trial 
amplitude flHI.BOj) . 

But now that both the kinetic energy E kin and the gauge field energy are at hand 
in their improved versions, cf. f lVI.lOp -f lVI.lip and fTVI.22j) . one can add them together in 
order to find from the energy functional E^ (TVI.4|) the corresponding energy function 
E^ Iv ^(/3; v,po,pi) which then appears in the form 

e 2 

E^\f3-u, Po ,pi) = — {(2a B f3) 2 ■ e kin (u,p Q ,Pi) - (2a B (3) ■ e pot (u,p ,pi)} . (VI.23) 
a B 

Of course, this is again the p-generalized form of the simple predecessor E^ IV ^ v) flIII.55p . 



2. Principle of Minimal Energy 

With an improved energy function E^ Iv ^(/3; v,Po,pi) being at hand now, one can next 
apply the principle of minimal energy flIII.39P - flIII.40p in order to look for the extremal 



points of that function. First, turn here to the determination of the first variational 
parameter (3 in terms of v and p by exploiting the first of the equations (1IH.40p . More 
concretely, the requirement 

a E W(/V,p , Pl)=0 (vi.24) 
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yields the equilibrium value of as 

2 £kin(^,P0,Pl) 

and if this is substituted back into the original energy function E^ Iv ^(/3; is,p ,pi) (1VI.23j) 
one obtains the reduced energy function Et(V, Po,Pi) m the following form: 

e 2 

E T (i/,p ,Pi) = E^\f3;u, Po , Pl ) =- p (v,p ,Pi) , (VI.26) 

p=p t 4a B 

where the improved spectral function S p (u,p , P i) is given now by 

S P WPcPi t= / r • VI.27 

£kin(^, P0,Pl) 

Clearly, this is again the p-generalization of the simple spectral function f lIII.56bl) which is 
due to the simpler trial ansatz (lIII.50p . Indeed, for p\ — » the present three-parameter re- 
sult S p (u, po, pi) (1VI.27P becomes reduced to its one-parameter predecessor S p (u) (]III.56b|) . 



3. Simplified Groundstate Demonstration 

The improving effect due to the presence of the additional variational parameters 
Po,Pi may be roughly estimated by fixing the parameter v in equation (1VI.27|) and thus 
discussing the spectral function S p (u,po,pi) exclusively as a function of po,pi- If one 
knew the optimal value of u, one could find in this way the desired extremal value of the 
energy function Et(z^, Po,Pi) flV"I.26[) . But for the moment we will be satisfied with only a 
moderate improvement of the results of Sect. Ill which are due to the simpler variational 
ansatz (1111.501) . i. e. we will base our improvement upon those values v* which are 
displayed in table 1. Recall that those results yielded a deviation of (roughly) 10% from 
the conventional predictions within the framework of the two-parameter amplitude $(r) 
(HTL5QD. 
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For the sake of a transparent but sufficiently detailed demonstration of our improve- 
ment, we will restrict ourselves here to the groundstate situation (n p = 1) and, properly- 
speaking, this would imply that we have to choose = —0.2049... (see table 1, first line). 
However, in order to keep our groundstate demonstration as uncomplicated as possible, 
we are satisfied with choosing v = 0. This choice reduces our four-parameter amplitude 
^bi^) (IVI-ip to a three-parameter trial ansatz which (for p\ = 0) incidentally predicts 
exactly the same groundstate energy as the conventional theory (II. 4p . see the discussion 
of this in ref.s [HJ [18]. It is true, a more realistic choice for the groundstate would be 
z/* = —0,2049 (see the first line of table 1), but this would render the intended demon- 
stration unadequately complicated, inasmuch the groundstate {n p = 1 <^ l p = 0) is in 
any case badly predicted by the present electrostatic approximation. 

By these arrangements, the two-dimensional problem of looking for the extremal values 
of the energy function Et(za Po,Pi) flVI.26j) - flVI.27j) becomes reduced to a one-dimensional 
problem, namely to look for the extremal values of the function E^(j9 ,Pi) 

e¥ 1 ( Po ,Pi)=E t (z/,po,Pi) • (VI.28) 
Or equivalently, this means to look for the extremal values of the reduced spectral function 
Sp (po,Pi) = S p (0,p o ,pi) 

Sl% , Pl ) = £lot ^ P0,Pl \ ■ (VI.29) 

% £kin(0,p ,Pl) 

Here, the general kinetic function ^km^, Po, Pi) fIVI.llj) becomes reduced to (observe l p = 
0) 

£kin(0,R),Pl) = ^Po + 2 Pi ^ ^kin(Po,Pl) (VI.30) 

and finally the reduced potential function £ po t(Po,Pi) 4= £ po t(0,Po,Pi) is deduced from 
(^PcPi) QVI-13P as follows 

£pot(po,Pi) =p^£pot(0)+4p^ 1 -£ I (0)+4p^-£ii(0)+4p p?-£iii(0)+^-£iv(0) , (VI.31) 



i.e. by reference to the table in App.G 

,»..,(/',,./'!> = ^Po + 3PoPi + yPoPi + yPoP^ + ^Pi • <YU2} 
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Thus the reduced spectral function S p (p) f]VI.29|) is ultimately found to be of the following 
form 

q[0] ( gpot(RbPi) 2 _ (pj + 6 PoPi + 17 ^i + 25 ^i + f Pi)' n n 

S v {Po ' Pl) ~ sUpo,pi) ~ pTT^pJ • (VL33) 

For determining the local extrema of this spectral function flVI.33j) . one has to observe 
the constraint (IVI.2j) . For our present choice of v = 0, this constraint implies the restriction 
to a one-dimensional compact subspace of the two-dimensional IR 2 : 

Pi + 4p Pi + 6pl = 1 • (VI.34) 

Since such a closed subspace is topologically equivalent to a circle S 1 , one may introduce 
a circular coordinate (a, say) in order to parameterize both real numbers po,Pi'- 



Po = cos a 



sin a (VI. 35a) 

Pi = ^ ■ (VL35b) 

Indeed, one is easily convinced that the constraint (I VI. 34 p is automatically satisfied by this 

parametrization. Furthermore, the kinetic energy function £ r kin(po ) Pi) (IVI.3Q[) becomes 
now also a function of this angular parameter a, i.e. 

£km(po,Pi) =3- £km(a) = -g (l + 2\(a0) (VI. 36a) 

Ti(a) = 2 sin a ^sina — v^cosaj . (VI. 36b) 



The same arguments do also convert the potential function £ po t(Po?Pi) (IVI.31 [) into a 
function of a 

Spot (Po, Pi) Spot (a) = g ( X + ^i( a )) (VI.37a) 

(1 \[2 7 \ 
- sin cH cos a ■ sin 2 a — \f2 cos a sin 3 a) . (VI. 37b) 
2 4 8 J 
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Consequently, the angular dependence of the spectral function Sp\a) f]VI.33|) must look 
as follows 

(l + Vxia)) 2 

S® (Po, Pi) =► 4? («) = \ ^7 / • (VI-38) 

1 + li{a) 

This function has been discussed in some detail in ref. [2] (see equation (IIV.44I) of 



that reference); and it has been found that there occur two maxima and two min- 
ima which then generate four energy values E\ a . . . Eid in the following arrangement : 



e 2 



rel. maxima: „ 



E la = — 1 ~ -6, 8029 . . . [eV] (VI.39a) 



c 2 



E lh = — 0, 086717499 . . . ~ -0, 590 [eV] (VI.39b) 

4ae 



e 2 



i • • 4aR 
rel. minima: „ 



E lc = — 1, 033319474 . . . ~ -7, 030 [eV] (VI.40a) 



e 2 



E ld = — 1, 128194657 . . . ~ -7, 675 [eV] . (VI.40b) 

4a B 

Here, the first result E\ a flVI.39al) is clearly the most striking one, because it agrees 
exactly with both the zero-order approximation (a = 0) and the conventional ground- 
state prediction (11.41) . where the latter is commonly considered to be exact within the 
framework of the generally accepted non-relativistic quantum mechanics. However, the 
occurrence of two minima and two maxima in the first-order spectrum demands an expla- 
nation: From the physical point of view there should occur for N = 1 just one maximum 
and one minimum, corresponding to the groundstate and the first excited state. But evi- 
dently the present approximation procedure generates spectral functions with additional 
spurious extremal values which are not reproduced in the next higher approximation 
steps. Therefore it may appear desirable to get some confidence in the presently used 
approximation formalism. This is attained by reproducing the above results through a 
different parametrization of the trial amplitude (App.H). 

Summarizing the simplified groundstate treatment, one gets a first estimate of the 
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groundstate energy at E ld = -7, 675 [eV] flVI.40bj) which deviates in an unacceptable 



way from the conventional prediction E^L = — ^j— — —6,8029 [eV], (cf. equation ( II. 4j) ). 
However, this does not mean that RST is unable to reproduce the conventional positro- 
nium spectrum fll.4j) . As we will readily see, the spectrum is actually reproduced by the 
present approximation formalism with sufficient accuracy as far as the excited states are 
concerned. The groundstate seems to represent an exceptional situation which must be 
studied separately. 

4- Excited States 

For the treatment of the excited states one has of course to relax the condition v — 
which may be applicable exclusively for the (simplified) groundstate. Or otherwise, if one 
wishes to stick to v — 0, one has to further generalize the present trial ansatz (IVI.lj) in or- 
der to work with the more general hydrogen-like wave functions (the latter approximation 
method has been applied in ref. [14]). However, for the present elaborations, we rely upon 
the simpler variational ansatz (IVI.lj) with general v and integer quantum number l p of 
angular momentum (IIII.21j) . This approach must generate the same energy spectrum as 
does the method of the hydrogen-like wave functions in combination with vanishing angu- 
lar momentum l p = 0, provided the angular momentum degeneracy emerges in RST quite 
analogously to the situation in the conventional theory based upon the Hamiltonian H 

o. 

Thus our procedure will be a slight generalization of the precedent treatment of the 
simplified groundstate; i.e. for parametrizing the general spectral function S p (v,po,pi) 
(IVI.27j) (with non-zero variational parameter v) by the S* ^coordinate a we first have to 
look for both functions ekin^PcbPi) (jVI-lip and ^pot^, Po^Pi) HVI. 13[) as functions of that 
new coordinate a, and afterwards one determines the maximal value of the corresponding 
function S p (u, a). 

S P ^a) = S i^\. (VI.41) 
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For such a reparametrization by the angle a it is very convenient to first eliminate from 
both functions Evwiy-, Po,Pi) an d £ P ot(^,Po,Pi) the parameter p as far as possible since its 
reparametrization is somewhat more complicated than its associate p\. 



Po 
Pi 



Po{a) 
Pi (a) 



cos a — a/2(z/ + 1) • sin a 
sin a 



(VI.42a) 
(VI.42b) 



Clearly, this is the generalization of the predecessors ( 1VI.35ap - (1VI.35bj) for the simplified 
groundstate [v — 0); but indeed it guarantees the validity of the more general normaliza- 
tion condition flVI.2j) . Thus, the kinetic function £kin(^, Po,Pi) flVI-lip reappears now in 
terms of the S^-coordinate a as 

1 



£kinW P0,Pl, 



2v(2v 



+ sin a ■ \v + 




sm a cos a 



u + 1 



(VI.43) 



Evidently, this result becomes reduced to the former case £kin( zy ) (HII.53P for a — > 
(i.e. pi — > 0,po — > 1); and it becomes reduced to the simplified groundstate situa- 
tion e kin {a) flVI.36aj) - (IVI.36bjl for l p = 0, v 0. 

The next step must consist in looking for the reparametrized numerator e pot (i/, a), of 
the spectral function S p (u, a) (1VI.41|) ; i.e. we have to elaborate the transition £ P ot(^Po,Pi) — * 
£pot{v, a) where the original potential function £ P ot( z/ ) Po,Pi) is given by equation flVI.131) 
and the angular parametrization of po,p± by equations (IVI.42aP - flVI.42bp . The result then 
looks as follows: 

EAv) 



EpotW P0,Pl) 



-poU 



M + ejM + 




(VI.44) 



'2(i/ + l) 

Clearly, this is the generalization to arbitrary v of the former potential function e po t(a) 



flVI.37aP - flVI.37bp which in turn is recovered from the present more general result (1VI.44I) 
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by putting v to zero [v — > 0). This claim can easily be validated by observing the 
definition of the auxiliary functions (see also the table in App.G) 

ei(i/) = ei(i/) - 2(i/ + 1) • 
e 2 (v) = Ae 

=£mM-4(i/ + l 



£ P ot(^) 

,,(/') - 16(i/ + 1) ■ ei(i/) + [I6(i/ + l) 2 - 4(i/ + l)(2z/ + 3 



(VI.45a) 



■ e u {y) + 2(6z/ 2 + llu + 5) • £l (u) 



■e pot {v) (VI.45b) 
(VI.45c) 



£4(1^) = 



- 4(2z/ 3 + 5z/ 2 + 4i/ + 1) • e pot (i/) 
eiv(^) - 8(1/ + l)(2i/ + 3) ■ e n (i/) + 32(i/ + l) 2 (2z/ + 3) ■ £l (i/) 

- 4(1/ + l) 2 (2z/ + 3) (61/ + 5) • e pot (u) 



(VI.45d) 



Thus, putting v to zero yields 



ei(0) = 
£ 2 (0) = 

£ 3 (o) = ; 



1 

4 
3 

2 



£4(0) 



3 
1 



(VI.46a) 
(VI.46b) 
(VI.46c) 
(VI.46d) 



and this simplifies the general e pot (i/, a) flVI.44j) actually down to the former £ r po t(a) 
flVI.37aD - flVI.37bj) . 

With both the numerator e pot (z/, a) and denominator e^^v, a) of the spectral func- 
tion S p (u, a) (1VI.41P being known now, one can let search the maximum S}p\*) of S p (u, a) 
by means of an appropriate numerical program and finally can then tabulate the corre- 
sponding positronium spectrum Ep 

^2 



EN 



4a B 

4 n] (*) : = sj^ ] ,pt ] ) 



(VI.47a) 
(VI.47b) 



see table 3 below. 
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As a check for the correct working of the program, one can alternatively determine the 
energy spectrum also in the reparametrized form (1H.26aj) - (1H.26bj) which must identically 
reproduce the spectrum due to the present form (lVI.47aP - (1VI.47bj) . 



Tip 


J^conv 




An] 
Is* 


N 
p* 


t±p\V* ) 


dVT 47a 1 flVT 47h \ 
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5 
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n oonoc 
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p 
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-0.02001 
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Q 
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-U.lUboU 


n nn/i oq 
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-U.U1 / oo 


11 QtC 
11. OO 


9 


-0.08399 


-0.07470 


17.00395 


-0.01532 


11.06 


10 


-0.06803 


-0.06070 


19.59761 


-0.01367 


10.77 


11 


-0.05622 


-0.05033 


22.26932 


-0.01232 


10.49 


12 


-0.04724 


-0.04241 


25.01475 


-0.01120 


10.22 


13 


-0.04025 


-0.03624 


27.82998 


-0.01025 


9.97 


14 


-0.03471 


-0.03133 


30.71171 


-0.00944 


9.74 


15 


-0.03024 


-0.02736 


33.65704 


-0.00874 


9.51 



Table 3: Improved RST predictions (Spherical Symmetry) 

The improved RST predictions Ep (IVI.47aj) (spherically symmetric approximation) 
do differ essentially from their conventional counterparts Ecom f ll.4|) if compared to the 
corresponding "anisotropic" improvements of table 2 (p. [77|) . Defining here the devia- 
tion A[ p ] quite analogously to the former cases A[x] ( 1III.58j) and A{ T } (1V.44P as 



Ar 



1 U 

-V 

14 ^ 



(n) 



n=2 



(n) 



']4f A >i n, )> ( vi - 48 ) 

/ n=2 
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one finds for the present improvement of the spherically symmetric type 

A w = 10,7% (VI.49) 

This is an irrelevant improvement in comparison to the analogous average of 10,9% due 
to the simplest approximation (table 1, 2 ^ n p ^ 15). Observe also that for the ground- 
state (n p = 1) one finds ~ which is preliminarily assumed for the sake of simplicity in 
order to exemplify our general approximation scheme, see the precedent Section VI. 3. 
Through this fortuitous agreement, the former groundstate prediction Ei& flVI.40bj) is 
practically identical with the result Ep of the present more rigorous calculation (first 
line, n p — 1). 
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VII Conclusion 



The present results may give to us now a hint at the direction for searching further 
improvements. It is true, both higher approximations (i.e. the "anisotropic" one of table 2 
and the "spherically symmetric" one of table 3) yield certain improvements in the sense 
that the corresponding RST predictions do better approach the predictions E^olv (11-41) 
of the conventional theory. However, the improvement of the spherically symmetric kind 
(table 3) are inconsiderable because they amount to merely some 0,2 percent. Thus 
the conclusion is that the original results of the spherically symmetric type (table 1) do 
already exhaust the possibilities of this symmetry. But the "anisotropic" improvements 
(table 2) do reduce the RST deviations (from the conventional predictions) by some 50% 
in the range 2 < n p < 15. Therefore one expects that the next higher approximation 
step of the anisotropic type (i.e. the regard of the radial auxiliary potential ^ p ^ v (r) 
in the expansion (1V.3|) ) will show a further clear step of the RST predictions towards 
the conventional energy spectrum (11.41) . Here, the regard of the anisotropy can not be 
restricted to the interaction potential ^Aq alone but must refer also to the wave functions 
themselves. More concretely, this means that the product ansatz (IIII.17al) - (IIII.17bj) for 
the wave amplitudes must receive some generalization. Taking all these results together, it 
seems that the conventional levels (11.41) play the part of a lower bound for the successively 
higher RST approximations. In any case, further clarification of this question must be be 
left to a more extensive investigation. 
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Appendix A: 
Spherically Symmetric Potential 



The spherically symmetric approximation is defined by assuming the electric interac- 
tion potential ^A (r) ( III.30p to be SO(3) symmetric, cf. equation (IHI.ip . The principle 
of minimal energy (III.40P then yields the Poisson equation (IIII.36P for the symmetric 
potential ' p l4o(r) which may be recast to the following dimensionless form 

d?a v (y) 2 dd„(y) _ y y 2 ^ 1 
dy 2 y dy 6 F(2u + 2) 1 ' ' 

(y = 2/3r) , 

provided the original potential t p l4o(r) (IHI.ip due to the trial amplitude $(r) is rescaled 
as follows 

[ p U (r) = 2f3a s ■ ~a v (y) , (A.2) 
and similarly for the trial function $(r) flIII.50p 

<£(r) = J- ■ . 2/3 $„(y) (A.3a) 
®»(y) = y u ■ e~> . (A.3b) 



The solutions a u (y) of the dimensionless Poisson equation ( lA.ip can be represented 
for integer values of 2v in form of a recurrence formula: 

e~ y ^ y n 

< +l (y) = My) - (2z/ + 2)(2zy + 1) -]L- r (A- 4 ) 

with the lowest-order potential a (y) for v — being given by 

a (y) = - (1 - e~y) . (A.5) 

y 



95 



Thus, the lowest-order potentials (i.e. for v = 0, ^, 1, |) are found to be of the following 
shape 

01 (y) = ~a (v) ~ ^ = - { 1 - e- y ( 1 + |) | (A.6a) 



2 y I V 2 

o 1 (y) = 2j(y)-^(l + y) = i|l-e-*^l + ^ + ^)} , (A.6b) 

e " ^ ^ 1 r, .^^31,1, 



3^^ = 01^)- — { 1 +y+j) = ~{ 1 - e ~ v { 1 + zy + iy +2iy )] > ( A - 6c ) 

see fig. A.I below for a sketch of these lowest-order potentials. For general (but still 
integer) values of 2v the dimensionless potential a v {y) t ]A.4jl looks as follows: 

Here it must be stressed that such a potential a„(y) does actually exist for any real- valued 
> — =), where the general expression for this situation is then given by an infinite sum 



v \ 



oo 

n 



a„(y) = —!— [l-e-y-Y y 2v+n , (A.8) 

yy ' 2v+l\ f^T{2v + 2 + n) y I ' v ' 

see appendix D of ref. [TT] . 

The value a v (0) of the potentials a v (y) (IA.5j) - (IA.6cj) at the origin (y = 0) is easily 
found by expanding the exponential function e~ y in the usual way which then yields 

5o(0) = 1 (A.9a) 

3i(0) = i (A.9b) 

ai(0) = i (A.9c) 

5|(0) = i, (A.9d) 

see also the sketch A.I of potentials a v (y) below. Obviously the general case is given by 

4,(0) = ^ ; (A.10) 
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and this result does hold for any real- valued power v > — \. Indeed, this claim can easily 
be verified by exploiting the general shape of the solution of the Poisson equation (lA.lj) 

1 f p~y' ■ v ,2l, ~ l 

aJy) = / d 3 y' ^ ,, (A.ll) 

vKU> 47r.r(2i/ + 2)7 \\y-y'\\ 1 ; 

which satisfies also the usual boundary condition at infinity (y — > oo) 

lim a v {y) = - Wu . (A.12) 

y-yoo y 

In terms of the original potential t p U (r) (1A.2j) this boundary condition ( ]A.12j) yields just 
the well-known Coulomb form 

lim [p U (r) = — (A. 13) 

t— >oo T 

which thus is found to be the asymptotic limit form for all potentials H4 (r) being due 
to the class of trial functions $(r) (1111.501) . irrespective of the special value of v. 
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-1.4- 



groundstate (v = -0.2049...) 



Fig. A.I: Spherically Symmetric Potentials j\A.5\) - ljiA.6c\) 

The potentials a v (y) (IA.8[) are finite at the origin (y — 0), cf. flA.lOp . and approach the 
Coulomb potential (~ ^) for y — > oo, cf. ( 1A.12|) . For the excited states the parameter v 
is positive (see the table 1 on p. HH ) and therefore the potentials a u (y) have vanishing 
derivative at the origin (y = 0), but for the groundstate (z/* = —0, 2049 . . .) the derivative 
becomes infinite, see equation (IA.20j) below. 
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The fact that the spherically symmetric potentials H4o(r) are finite at the origin 
(r = 0) does of course not mean that the corresponding field strength E p felt by the 
electron 



EJr) = V [pl A 



o [r 



dr 



(A. 14) 



is in any case a regular vector field for r — > 0. Indeed, introducing the dimensionless 
form e u (y) of the radial component ^E r (r) 



through 



E p (f) = ®E r (r)e r 



(A.15) 



(A.16) 



one arrives at the following link of (dimensionless) potentials a v {y) and field strengths e v (y): 

da u {y) 



dy 



(A.17) 



But since the spherically symmetric potentials d v (y) are explicitly known both for 
integer 1v (I A. 71) and general v (IA.8j) . the field strengths e u (y) are also known, e.g. for 
integer 1v: 

(A.18) 



My) = -4fi-^-E ? r > ) > 



or similarly for general v 



e v (y) 



E 

n=0 



,2i/+n 



r(2z/ + 3 + n) ' 



(A.19) 



From this latter result it is now easily seen that the field strength does behave at the 
origin as follows: 



2,' 



hme^) = -hm f ^ y 



— oo , v < 
, i/ > . 



(A.20) 
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This says that the field strength E p (r) (1A.14I) is a non-singular unique vector field only 
for v > 0, see the figure A. II below. Thus it is only for the groundstate v ~ —0, 2 that 
the electric field strength turns out to be singular at the origin (y = 0). Observe also that 
(for all v) the field strength e u (y) adopts the Coulomb form (broken line) at infinity, i.e. 

lim e v {y) = -\ , \/u (A.21) 

quite similarly as for the potential a u (y) (IA.12I) . 
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1 2 3 4 5 6 



j i i i i i 




Fig. A.II: Electric Field Strength e„(y) QA.18l) - (A.l&) 

The electric field strength e v (y) is non-singular only for the excited states which have 
positive variational parameter v. However, for the groundstate [v ~ —0.2) the field 
strength becomes singular, albeit less singular than the Coulomb field (~ — p). This 
guarantees that the energy content fill. 361) of the gauge field remains finite. 
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In the present context of positronium binding energy, the essential point refers now 
to the value of the non-relativistic energy functional Km (1III.41I) upon the class of trial 
functions $(r) 50[) . Indeed, the desired approximate spectrum of binding energies 
emerges as the set of extremal values of Er$i on that chosen class of trial functions $(r) 
flIII.501) . More concretely, the value of E[$] upon that class yields a function EP v J(/3, v) of 
the two variational parameters f3 and v which looks as follows [TT] : 

E M(p, v) = — [{2a B I3f ■ e ki » - (2a B 0) • e pot (v)] (A.22) 

h 2 \ 



Me 2 , 

The extremalization procedure due to the principle of minimal energy (III.40I) requires now 
to look for the stationary points of this function E*- Iv ^(/3, v): 



d(3 



(A.23a) 
. (A.23b) 



dv 

It is easy to see that these two requirements can be reduced to the one requirement 

-it 1 = < A - 24 > 

with the reduced energy function Et(^) being defined through [11] 

E T (u) = ■ S p (u) . (A.25) 

The spectral function S p (v) emerging here is itself defined in terms of the variational 
parameter v through 

S p {u) = , (A.26) 

with the kinetic function e\a n {v) depending upon the variational parameter v and upon 
the angular momentum quantum number l p (1IH.33|) in the following way [11] 



e^{v) = ^—(- A + l ^\ ■ (A.27) 



2z/ + 1 \4 2v 
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Furthermore, the potential function e po t(u) is found to be somewhat more complicated 
but can be represented by virtue of the Poisson identity in two equivalent ways [11] 

L T(4v + 3 + m + n) 



-pot l 



2^+3 z-^i 2 m + n T(2v + 3 + m) ■ T(2u + 3 + n) 

(A.28) 

if 1 n r(4z/ + 2 + n) 

" 2z/ + 1 I ~ 2^+ 2 ' ^ 2" ' T(2z/ + 2) • T(2z/ + 2 + n 

It follows from the present results in combination with the principle of minimal energy 
(III40I) that the non-relativistic positronium spectrum may be obtained by looking for the 
stationary points of the spectral function S p (u) (1A.26I) for any quantum number l p . But 
since our trial function $(r) (1III.50j) is too simple in order to get the whole spectrum of 
excited states due to the one chosen l p , one will obtain in this way merely the "ground- 
state" for any l p , cf. figure 2 of ref. [11]. The result of this extremalization procedure is 
presented by table 1 on page HH 
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Appendix B: 

Auxiliary Potentials {p} A|, {p Um, {p} A v 

for general v 



For an estimate of the radial auxiliary potentials (lIV.llaP - (IIV.lle|) due to general v it 
is most convenient to first recast these potentials to a dimensionless form, quite similarly 
as it was done for the spherically symmetric approximation l p l4o(r), cf. equation (IA.2j) ; 
i. e. we put 



Wa(y) = (2/^r 1 - feU,(r)= / dy' - 

Jo [y 2 + y> 2 ] 



«%(y) = (2/3« s )- 1 - ^A w (r) = y 2 - / dy' - 

Jo [yi + y> 2 ] 



\{y) = (2/3a s )~ l • H4 v (r) = y* - / dy' V - • tftf 



o [y 2 + y' 2 ] 



(B.la) 
(B.lb) 
(B.lc) 



$(r) 



y = 2(5r) . 



The differential links flIV.20p - fllV.2ip read now in this dimensionless form 

{ ^{y) = -\y^(y- { Hy)) 

(Vfc( y) = -L ylA- 

yy) 35 y dy 



1A(V- (Ill)-/ } 



(B.2a) 
(B.2b) 



For a treatment of the excited states it is necessary to admit the variational parameter 
v in the trial amplitude <3>(r) flIII.50j) to adopt quite general (real) values, not only the 
special value v = which is adopted in Sect. IV. 2 for the simplified treatment of the 
groundstate. For such a more general situation, one may still take as the starting potential 
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ip}A\(r) again the spherically symmetric potential H4o(r) which obeys the spherically 
symmetric Poisson equation (1III.36j) but now with non-zero v: 

I f£ + 2 d \ M = _ 2/3a a (2/3r)2 ,-i e - 2/3 , ^ (B 3) 



(2/3) 2 V^ 2 rrfry uv 7 T(2z/ + 2) 
cf. also the dimensionless form (lA.lj) of this Poisson equation and its solution (1A.8j) . Thus 



our starting potential ^a u (y) is now just that spherically symmetric approximation ( ]A.8j) : 
i. e. we put 

(% v {y) = - (2/3a s )-i bUo(r) = ( 1 - e~» £ t^-^-t I ■ (B.4) 

V n=0 



vr y w tt2i/+1\ ^T(2z/ + 2 + n 

\ n=0 

For integer values of 2z/ (i.e. v — 0, |, 1, |,2, |, ...) this starting potential appears as a 
finite sum, namely [TTj 

l " T»l 2„+l „|J 



Clearly, for z/ = one recovers here the former starting potential ^A\(r) f)IV.25[) . 

But now that the more general starting potential ^a v (y) has been fixed, one can use 



it in order to calculate the associated third auxiliary potential ^ l % u (y) ( 1B.2aj) : 

run , n 1 d 2 / nv , ,\ 2 y^+V^ ,„ . 

As a consistency check, one puts here v = which yields 

(IM) ao(l/) = |- ye"* = A (2/3r) e" 2 ^ , (B.7) 

and thus the third radial potential ^A\\\(r) for v = becomes 

H4„,(r) = 2/3a s • V% (y) = A a s /3 2 re" 2 ^ (B.8) 

which is nothing else than the former result (!IV.26aP being used for the simplified treat- 
ment of the groundstate. Once the third auxiliary potential ^ w 'a v {y) is known, cf. flB.6j) . 
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one can proceed to calculate by means of it the corresponding fifth auxiliary potential 
{y) ~aM dE2b]): 



1 3 d 
35 ^ dy 



' 1 d 

_y A dy 
1 



V 



a>v{y) 



(B.9) 



105vrr(2z/ + 2, 
The consistency check for v = yields here 

16 



2u+l 



Mr. 



ao(y) 



2T 



2/ 



105tt I 8 2 
i.e. for the potential ^Ay{r) due to v — one obtains 

H4 v (r) = 2/3a s ^(y) = (2/3a s )e~^ - i (2/3r) 2 - \ (2/3r) 



(B.10) 



(B.ll] 



and this again is just the former result flIV.26bj) for the simplified groundstate. 

Furthermore it is also rather obvious that it is the auxiliary potential of lowest order, 



i. e. 



] &v{y) (IB.4p . which obeys the Coulomb-like boundary condition (IIV.18P at infinity 
{^a v {y) — > ^ • ^) whereas all auxiliary potentials of higher order (e.g. ^% u (y) (1B.6j) and 
i/(y) (1B.9P ) tend to zero more rapidly than the Coulomb potential (~ -). Similarly, 
the boundary conditions at the origin (y = <^ r = 0) are also satisfied: the higher- 
order potentials such as ("Ofi 

u{y) flB.6j) and ( v u u (y) flB.9j) do actually vanish at the origin 
(y = 0) and thus validate the boundary conditions (HV.16j) . On the other hand, it is 
only the auxiliary potential of lowest order (i.e. Cfe 

v{y) (IB.4p ) which adopts the required 

non-zero values (IIV.13P at the origin (y = 0). 

Summarizing, the expansion ( 1IV.9|) generates a hierarchical set of radial auxiliary 
potentials (1IV.1 lap - (1IV.1 left etc. of ever decreasing magnitude; for the three lowest- 
order potentials ®a v {y) flES), m a u (y) (ESJ) and W a u (y) (P , see the figure B.I below. 
With increasing order of this "anisotropy expansion" the deviation of the corresponding 
electrostatic potential ^A (r, ■&) (IIV.9P from the exact potential (IIV.5P becomes reduced 
step by step. In the present context, the importance of such an expansion refers to the 
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gauge field energy (1II.36P which undergoes a corresponding expansion with respect 
to the magnitude of the anisotropy and thus provides us with an estimate of the relative 
weight of the anisotropic terms (compare, e.g., the "anisotropic" groundstate energy 
flIV.45j) to its "isotropic" counterpart (IIV.46j) . where a deviation of only 0.4% arises). 
Though the treatment of the groundstate [y = 0) on the basis of the present anisotropic 
correction ^A\\\(r) (IB. 81) appears to be meaningful, the generalization to the excited 
states v > turns out to be unacceptable (see App. C). 
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2 

3k 



1 

371 



(I) 



a. 




Fig. B.I Radial Auxiliary Potentials ^ai, ^'"^ai, { v }<ii 

(EJfy, (EJfy and (E7W 



The hierarchy of magnitudes of the potentials ^a„(y), ^a v (y), ^a v (y), . . . persists 
also for v > (here v = 1) similarly as shown for v = by Fig. IV. A. Nevertheless the 
associated series expansion (IIV.9I) of the anisotropic potential { p bA (r, fails to describe 
the excited states correctly, see App. C. 



108 



Appendix C: 



Failure of the Expansion (ilV. 9\ ) for the Excited 

States (v > 0) 



As plausible as the result of the simplified groundstate {y = 0) may appear (see the 
discussion below (IIV.46|) ). this way of dealing with the anisotropic corrections cannot be 
transferred to the excited states! Indeed, we will readily see that the energy predictions 
of the excited states would come out in an unacceptable magnitude. Since the simplified 
groundstate discussion is based upon the "anisotropy expansion" (IIV.9I) . this expansion 
must be afflicted with a serious deficiency, at least as far as the excited states are con- 
cerned. Therefore, an improved expansion is presented in Sect. V, but for the sake of 
completeness we will subsequently work out the reduced energy function E*(z/) due to the 
(rejectable) series expansion (1IV.9j) . in order to see more clearly what is going wrong with 
this expansion. 

To this end, we compute the corresponding energy function E$(/3,i/) which is to be 
understood as the generalization of the groundstate function E$(/3) f)IV.43|) . such that 



Here it is convenient to scale off the atomic energy unit (a. u.) by introducing the spectral 
function Sp through 



so that we are left with the problem to elaborate this spectral function Sp (u) for general 
v. Clearly, its extremal values will then yield the wanted energy spectrum of positronium 
(though to be rejected). 



E*O3)=E*G0,i/) 



(C.l) 



E*(i/) 




(C.2) 



6.8029... [eV]) 
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Fortunately, it is possible to explicitly calculate that spectral function Sp\v) (10.2)) 
for integer 2u. More concretely, the general gauge field energy appears now m 

terms of the dimensionless auxiliary potentials ^a u (y) (IB.4j) and ( n % u (y) (IB. 61) as a sum 
of three terms, cf. equation flIV.32j) . with each of the three terms being found to appear 
in the following form: 



E^ = -teWa.) (|) 2 [dyy* (^j = -2& ■ 



4; } = -^c(2/3a s 



3tt 



'o 

2 /-oo 



160 



d yy s « 



dy 

7 fdV% v (y) 



dy 



(C.3a) 
(C.3b) 



+ 



-2/3e 2 ■ ("0 £pot(r/) . 

(C.3c) 



But since the radial auxiliary potentials ('^(y) and ^'"^(y) are explicitly known, cf. 
(|B.4p - (lB.6p . one can calculate here also explicitly the three auxiliary potential functions 
Spotty), ^"fepot(^) and C'fepotty) for general v. The result is 



0) 



eprt(f) = epotM , cf. (PQgp 



OO 



Spot I f J 



1 T(4z/ + 3) 



4!2 4 ^ T(2z/ + 2) 2 
9 + 7^ r(4z/ + 3) 



poti ; 192Q 2 4 - • T(2z/ + 2) 2 



(C.4a) 
(C.4b) 

(C.4c) 



These three terms add then up to the total "anisotropic" gauge field energy 

4 e} = £, {e} + ^ e > + J3fr> = -fic (2/3a s ) • 4 e >(^) 
where the "anisotropic"' potential function (v) is evidently given by 

M 



M 



(C.5) 



(C.6) 



Obviously, the present gauge field energy 

#r dHH) with (p4a|) - (p4cj) is nothing else 



than the generalization of the anisotropic groundstate result E^ pV.32[) - (lIV.33cj) to the 
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excited states. Observe also that the Poisson identity = (cf. (1IV.38P ) is obeyed 
strictly for the spherically symmetric approximation, but this identity is of course not 
obeyed by the present anisotropic approximation procedure, neither for the groundstate 
nor for the excited states! 

Therefore we have to calculate now the non-relativistic version M^c 2 of the (rela- 

(e) 
G 



tivistic) mass equivalent M^c 2 (III.37I) in order to build up again the constraint term N, 



(1IV.38j) - (1IV.39l) . but now for general values of the variational parameter v. This means 
that for the trial function <f>(r) in both parts Mf e} and Mjf fllV.40p - fllV.41l) of the mass 
equivalent one has now to make use of the general form flIH.50p - flIH.5ip . i. e. 

M[ c) c 2 Mf c} c 2 = - (|) 2 he J drr H4,(r) ■ $ 2 r 2 ^- 2/3r (C.7a) 

M^c 2 M^c 2 = - (j?J 2 he J drr H4,„(r) ■ ^r^e" 2 ^ , (C.7b) 

where the radial auxiliary potentials ^A\(r) and ^Um(r) are specified for general v (in 
dimensionless form ) by equations (IB . 5 1) - f lBT6|) of App. B. Thus one finds by explicitly 
carrying out the required integrations in f lC.7aP -f rC.7bp 



fc 2 = -he (2« s 0) ■ e pot {v) = £, {c} (C.8a) 



c 



hc(2a s f3)^(v) (C.8b) 
3 1 T{AU + 3) fC8c^ 



Clearly for v = 0, this result leads us back to the previous groundstate results (1IV.40P 



(lIV.4ip . But even for v ^ 0, the present results display a structure very similar to the 
former groundstate situation (II V. 42 p . Namely, the striking feature concerns here the 
Poisson constraint term Nq^ flIV.381) which for the present situation of the excited states 
adopts the following form: 

N§ } = Ef, c} + E^ 1 - M^c 2 = -he (2a s f3) ■ 4° } (C.9) 

with 

4V) ^ ( "^Pot(^) + (lll) ^pot(^) - /4°V) • (CIO) 
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This similarity comes about through the present identity relation (lC.4al) for the potential 
functions e pa t(u) and ^pot^) which entails for v ^ the same equality of the energy 
and mass equivalent M^c 2 , cf. the remark above (1IV.42[) . as is the case for the previous 
groundstate situation (1IV.40P with v = 0. Consequently, both electrostatic contributions 
(IC.5j) and (IC.9j) to the energy functional flIV.29j) add up in order to yield now 



Eg } + \)V ■ N£ } = -he (2a s (3) e£ } (v) - 2e G (u) = -he {2a s p) ■ e {c \v) (C.lla) 

e W(„) = e « - 2e G (jy) = % ot (u) - «% ot (v) - ( "% ot (v) + 2^\u) . 

(C.llb) 

The crucial point with this total electrostatic contribution (IC.llaP to the energy func- 
tion E$ (/?,!/), to be conceived as the value of the energy functional E[$] (1IV.29h on the 
chosen trial configuration, is now that it is of the same structure as in the case for the 
spherically symmetric approximation (1A.22|) and for the anisotropic groundstate flIV.43j) ; 



namely, both the gauge field energy and its mass equivalent M^c 2 do contribute a 
linear function of the variational parameter j3. On the other hand, the kinetic energy E^m 
has been found to contribute a quadratic function of /3, cf. (1III.52I) - (1III.53I) : and therefore 
the present energy function E$(/3, v) for the anisotropic excited states displays again the 
formerly encountered shape for the isotropic situation (cf. (1A.22|) ): 

E*(/3, v) = — [(2a B (5) 2 ■ e khl (u) - (2a B (3) ■ ^ c \u)] , (C.12) 
as 

where merely the simple potential function e pot (u) ()A.28j) is to be replaced by the present 
more complicated electrostatic function e^(y) (IC.llbj) . Clearly, as a consistency check 



one puts here v = and is led back to the former groundstate result E$(/3) (1IV.43j) . 

This similarity to the former configurations of spherical symmetry (Sect. Ill) and 
of the anisotropic groundstate (Subsect. IV. 3) enables us to write down immediately 
the result for the present anisotropic situation. The extremalization process (1111.40 j) . 
being exemplified for the spherically symmetric approximation through equations flA.23aj) - 



(lA.23b|) in App. A, reduces now the two-parameter energy function E$(/3, v) flC.12p again 
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to the wanted one-parameter energy function E*(i/) (IC.2I) . Here, the spectral function 
Sp iy) is formally the same as previously for the spherically symmetric approximation, 
cf. (lIII.56ap - (IIII.56bl) . with merely the former potential function e pot (u) being replaced 
now by its present "anisotropic" counterpart e^{v) (IC.llbj) . i.e. 



, 2 



( £ { e }(»)Y 

S^(u) = \ ) > . (C.13) 

Consequently, the last step consists again in looking for the maximal value of Sp\v) (y* 
minimal value of E*(i/) (j0.2j) ) by means of an appropriate computer program. 

For our present purposes, it may be sufficient to display the result for some representa- 
tive values of the variational parameter u, see the table below. For the sake of comparison 
with the spherically symmetric approximation, one can introduce here the relative shift 
(A) due to the present anisotropic corrections 

where the energy function Et(z^) due to the spherically symmetric approximation is de- 
fined through equations (1III.56aj) - (lIII.56b[) . From the results of this table it becomes 
obvious that the anisotropy effect lifts up the energy levels of the spherically symmetric 
approximation, where the shift increases with increasing variational parameter v. For 
instance, for v = 0, l p = one finds a shift of —0.41% (first line of the table); whereas for 
= 1.8 (i. e. the first excited state) the energy shift amounts to roughly 8% which is the 
same deviation as is found for the difference between the spherically symmetric approxi- 
mation and the conventional result (see the second line of table 1 on p. l4"l"|) . However, the 
unacceptable feature of this result is the fact that for the excited states [y > 0) the energy 
levels of the spherically symmetric approximation become lifted instead of lowered, as re- 
quired by the principle of minimal energy! Therefore the present anisotropy expansion 
(IIV.9P must be replaced by a more systematic procedure (~» Sect. V). 
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V 


7 


o v \v\ (1111. dodi 


\v\ (IL.13I) 


A(zy) (10. 141) [7o J 














1.0000 


1.0041. . . 


-0.41 


1 








i 

2 





0.7812. . . 


0.7715. . . 


1.24 


1 


0.1562. . . 


0.1543. . . 


1.22 


1 





0.6302. . . 


0.6089. . . 


3.38 


1 


0.2100. . . 


0.2029. . . 


3.38 


3 
2 





0.5278... 


0.4969. . . 


5.85 


1 


0.2262... 


0.2129... 


5.88 


2 





0.4546. . . 


0.4153... 


8.64 


1 


0.2273. . . 


0.2076. . . 


8.67 
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Appendix D: 

Exact Solutions of the Quadrupole Equation ( | V. 34\) 



For the construction of solutions ^A u \r) to the equation f lV.34j) it is very convenient 
to pass over again to dimensionless objects, i. e. we first define the dimensionless version 
of H4'"(r) through 

X\y) = (2(3a s )- 1 ^A u \r) , (D.l) 



so that the quadrupole equation flV.34j) reappears in the following form: 

6\ 10 y 2 "- 1 ■ e~ y 



tt T(2z/ + 2) 



(D.2) 



(y = 2/3r) 



A first integral of this equation can easily be found here by observing that the (ordi- 
nary!) differential equation (1D.2j) may be rewritten in the following form 



1 d 

y 4 dy 



10 y 



2u-l 



7T (2I/+1)! 



(D.3) 



dy v r 

where, for the present purposes, we restrict ourselves to integer values of 2u (i. e. 2v = 
1, 2, 3, 4, ...). By a first integration step in f]D.3|) one easily arrives at the intermediate 
result (for v > 1) 



10 . w . 1 
— (2z/ + 3)(2z/ + 2) - 

7T 



(D.4) 



dy v y* J n y 

which can be directly checked by means of a simple differentiation process. For the next 
integration step it is very advantageous to first consider the special value v = 1 for the 
variational parameter. In this case, one obtains from the intermediate result ( 1D.4|) 



7T 



T 4 ■ 



(D.5) 
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This result for v = 1 may now be used in order to further integrate in equation (ID. 41) 
which yields for integer 2v 

2v+3 



,„„ N (2z/ + 3)(2z/ + 2) mi. . 10, w , „ 1 /"°° , , 

ZU 71 n=6 n ' A 

(Z/>1) (D.6) 

(for v — 1, the second term is to be omitted). As a check, one lets act the differential 
operator (A y — p-) upon the result in order to find that the differential equation (1D.2|) is 



actually satisfied. In the last step, one explicitly calculates the residual integral on the 
right-hand side of equation flD.6j) and thus obtains the final result (for u > |) 

= {2 " + ^ + 2> (2, + 3)(2, + 2) ^e-» 2 f 1 g ^ (»-) (D.7) 

n=6 ' m=0 y 

(for z/ = 1, the second term is to be omitted). 

From this result it is clearly seen that the behaviour of the first anisotropic correction 
A^(y) is determined by A[ (y) flD.5j) in both the asymptotic region (y — > oo) and around 
the origin (y — > 0). Since Ai\y) (ID.5[) behaves as follows 



lim Al\y) ~ - • 1 , (D.8) 



one obtains for the general A„{y) (ID. 71) 



lim ^U'(y) ~ -(2i/ + 3)(2v + 2) • (D.9) 

2/->oo 7T 



which is nothing else than the claim (1V.36j) at the end of Sect. V. 3, rewritten in dimen- 
sionless form. 

In order to gain some insight into the properties of the anisotropy corrections A„ (y) 
in the vicinity of the origin (y = 0), one writes down explicitly some of the lowest-order 
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cases: 



^"(j,) = 40 y ^ 1 / 2 _ 5 3 5 4 _ 5 5 \ 
1 vy; 7T ^ n! 3tt I 6 y 14 y 48 y J v 7 

n=5 ' 

A "> = \ A "'^ ~ m yh ' y ^w^~h l + ■■■) < D10b > 

Observe here that with increasing v (— > 1, |, 2, ...) the anisotropic correction A^(y) 
becomes smaller and smaller in the vicinity of the origin (y — > 0), whereas at infinity 
(y — > oo) the correction A"\y) becomes larger and larger, cf. ( ID. 9}) . For a sketch of the 
three lowest-order corrections (lD.10a|) -f lD.10cj) see Fig. D.I below. 
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Fig. D.I Anisotropy Corrections A)J)(y) {ID. 6\) The exact solu- 
tion (y) ( 1D.6j) of the quadrupole equation (1D.2|) behaves like y~ 3 in the asymptotic 
region (y — > oo), cf. (1V.36j) ; whereas in the vicinity of the origin (y — > 0), the behaviour is 
for all v like y 2 , cf. (ID.10ap - (1D.10bj) . The magnitude of this quadrupole correction *4.i! (y) 
(< 1/10) is relatively small in comparison to the magnitude of the spherically symmet- 
ric approximation a u (y), cf. flA.101) . so that the anisotropy approximation can indeed be 
treated as a pertubation of the spherical symmetry, see the table of relative deviations A^j 
in App. F. 
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Appendix E: 

Electrostatic Energy E^ of the Anisotropic 

Configurations 



Since exact solutions of the "quadrupole" equation (IV. 34}) are at hand (see App. D), 
one can use them in order to calculate both the anisotropy energy Et } fTOOj) and its mass 
equivalent M an } c 2 fTOT]) which of course must be numerically identical. This identity is 
itself a helpful mutual check for the numerical correctness of the electrostatic field energy 
and its mass equivalent. Subsequently, we will present the results for some (half-)integer 
values of the variational parameter v which in Sect. V are then used in order to build up 
the "anisotropic" energy function E^ Iv ^(/3, u) for the excited states, i.e. the anisotropic 
generalization of the "isotropic" EP v l(/3, v) flIII.55j) . 



First, it is very convenient to introduce here again the dimensionless objects, i. e. the 
anisotropy energy E^ flV.301) reads in terms of the dimensionless potential correction 

A'M {ED 

Similarly, the corresponding mass equivalent M an } c 2 flV^T]) also reappears in terms of that 
dimensionless A^}(y) as follows: 

^ {c} c 2 = r dyy^e-y-Al\y) . (E.2) 



16(2z/ + l)!7 

The vanishing of the "anisotropic" Poisson identity N an ^ (1V.26|) says now that both right- 
hand sides of equations (lE.ip and (IE. 2D must be numerically identical! The subsequent 
table at the end of this appendix displays the results in units of (e 2 (3) for some values of 
the variational parameter v. 

Now in order to deduce a practicable recipe for calculating the anisotropy energy 
(lE.ip and its mass equivalent Mll^c 2 flE.2j) it is very convenient to refer the general case 
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(for integer values of 2u) to the special case for v — 1, quite similarly as it has been done 
for the calculation of the anisotropy corrections A l l l (y), see equation (1D.7j) of App. D. 



First, one recasts that equation (ID.7P to a somewhat more concise form, i.e. we put 

4»(y) = (2i/ + 3)(2i/ + 2). {^M-^B^y)} , (E.3) 



20 7i 

with the short-range constituent £>j,"(y) being defined through 



2u+3 n-6 , m 

B^)*^-.^!^!-^), (E.4) 



TV ^ 6??/^' 
n=6 ' m=0 y 



and obeying the following differential equation 



a »-fJ b "<^^-^ (e - 5) 

(w #!,"(?/) = for z/ = l) . 

Next, one could substitute this new form A„ (y) (1E.3|) into the equation (1E.1|) and 
would then find the desired anisotropy correction E^n being expressed in terms of Ai\y) 
and B"\y). We do not present here that equation, but rather prefer to calculate the 
corresponding mass equivalent Mli^c 2 (1E.2j) . namely on behalf of its more convenient 
handling. Indeed, if we introduve two auxiliary objects Kj, and ~h u through 

K v = / dyy 2 » +1 e-y -B™(y) (E.6a) 
Jo 

L ^40(2^i (E-6b) 
then the mass equivalent Mll^c 2 (IE.2[) reads in terms of these newly introduced objects 
IV = _^.(^ + 3)(2, + 2) f (2^2)1. u _ 5.^1 



(2i/ + 1)! 
-e 2 /3-/iIn } M- 
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The straightforward integration yields 

K - V* nl 1 (2z/ + 3 + n-m)! 

^ ( n + 6 )! 2 2-+4+n— ' ( n -m)\ [ ' 

n=0 v y m=0 K ' 

K = l-t , * . ( ^- 2 + "". (E.8b) 
Z^ n i 2 2^i+n 2z/-2! v ; 

n=0 v ' 

The subsequent table presents the value of the "anisotropic" mass equivalent Mli^c 2 
( 1E.7P for some representative values of the variational parameter v. Since the "anisotropic" 
Poisson identity vanishes for our exact solution A„\y) flE.3|) - flE.4|) of the (dimen- 
sionless) quadrupole equation flD.3j) . this table simultaneously displays the "anisotropic" 



gauge field energy Ein flE.ip for those specified values of v. This then admits to deter- 
mine also the total gauge field energy (IV.2ip as the sum of its "isotropic" part E^ 
( TV. 15ft and its "anisotropic" counterpart flE.ip : 



4 e} = 4 cl + Eg = -(e 2 /3) {2e pot (") + 

= -e 2 /3 • 2el:}(u) (E.9a) 

4SM epotM + ^SV) • (E.9b) 

In order to better estimate the magnitude of the anisotropy corrections, it is useful to 
introduce the relative magnitude A&£(y) through 

= ^j = = t • (E.10) 
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It is instructive to see from this table that the anisotropy correction in the vicinity 
of the first excited state n p = 2 (which has ~ 1.8, see table 1 on p. 144)) lowers the 
electrostatic interaction energy of the spherically symmetric approximation by some 
3%. The smallness of this anisotropy effect is the reason why the spherically symmetric 
approximation produces preliminary predictions with acceptable exactness, i. e. roughly 
8% deviation for the first excited state n p = 2 (second line of table 1 on p. |4"4"|) . Ob- 
serve, however, that the present anisotropy corrections (~ 3%) must improve the RST 
predictions for the binding energy by a factor of two (i. e. roughly 6%) because the total 
binding energy Et amounts to only half the value of the interaction energy E$ . For this 
reason, the deviation of the RST prediction for the first excited state (n p = 2) becomes 
reduced from 8.8% (table 1) to 3% (see App. F). 
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Appendix F: 

Extremalization Process ( 1 V. 38al) - (V. 38b\ ) for the 
Anisotropic Configurations 

If all constraints are duly regarded, the original energy functional Ep] ultimately con- 
sists of two terms only, cf. (1V.37[) . i.e. the kinetic energy E kin (1III.52p - (1III.53l) and the 
gauge field energy (1V.21|) which is the sum of the "isotropic" part (IV. 14}) and the 



"anisotropic" part Ein flV.30j) . From the principal point of view, both kinds of energy 
-Ekin an E^ are of equivalent conceptual importance; but from a more technical viewpoint 
the gauge field part E^ (and especially its anisotropic constituent E$ (1V.30I) ) requires 
greater effort for its exact calculation (for general v). Therefore we are satisfied for the 
moment with specifying E^} for integer 2u only, see App. E. But this then implies that 
also the total energy function E^ Iv ^(/3, v) (1V.39|) and its reduced form E^} (1V.40I) can 
be specified only for integer values of 2u. On the other hand, for the extremalization 
procedure (IIII40I) . as being specified by equations flV.38al) - flV.38bl) . we need a continu- 
ously differentiable energy function E{^} whose stationary points ( } v} = 0) determine 
the wanted energy spectrum. Therefore, as an auxiliary construction, we are satisfied for 
the present purposes with a supporting polynomial (of appropriate rank) which has in 
common with the spectral function Sp^(u) ( 1V.42[) its values on some suitable integers 1v. 
The extremalization process ( dSp , — = 0) is then carried through with the help of this 



supporting polynomial. 

For given quantum number l p (= 0, 1, 2, 3, ...) of angular momentum it will not be 
necessary to take into account the values of ' (y) on all integers 2v but one merely 
considers a limited number of supporting points (6, say) around the expected stationary 
point of Sp a \i>). For an estimate of the latter stationary point one may resort to table 1 
(p. HI]) because the anisotropy correction will be sufficiently small so that the stationary 
points due to the spherically symmetric approximation will not be shifted too far away. 
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Finally, it may also be interesting to see the extent to which the "isotropic" energy curve 
Et(^) (1III.56aj) - (1III.56bp becomes lowered to the "anisotropic" Ef\ (1V.41|) as consequence 
of the anisotropy effect. To this end, one defines the corresponding relative energy shift 
A{^} through 



A 



E 



M 



M 



E T (z/) 



"pot 



|E T (z/) 
» 



4ot0) 



SJv) 



1 - 



-{e}, 



^pot(^) 



(F.l) 



see the table below. 

Perhaps it is instructive to exemplify the proposed extremalization process by consid- 



ering in some detail the first excited state 



2 f~* L 



1). In order to write down 



the values of the spectral function S p (v) (IV.42|) on the integers 2z/, one takes over the 



potential function e^(v) from equations (lE.9al) - (lE.9bj) of App. E in combination with 
(1E.7|) . whereas the kinetic function e^iy) has already been specified by equation (lIII.53j) . 



In this way, one sets up the subsequent table for the values of the spectral function 



{a}, 



V) 



(on the integers of 2z/) in the vicinity of its extremal value 



1.7942 for the first 



excited state {n p = 2) in the spherically symmetric approximation, see table 1 on p. 14*4"]) . 



For an evaluation of the present RST results recall also the conventional prediction ( II. 4p . 
which for the value of the spectral function due to the first excited state (n c = 2) would 
yield (|) 2 = 0.25. (Compare this to the subsequent values of Sp^(v) for v — | and v — 2) 



(TOT]) 



(ED) 



( 181 y 
V384/ 

: 0.22217 



-5.76% 



12 



( 1533 \ z 
7 V 4096 / 

~ 0.24013 



-6.13% 



!()•(- 



10240 7 

0.24200 



-6.44% 



10 ( 65533 V 
3 ' V 245760 7 

~ 0.23701 



-6.70% 



84 



133933 V 



0.22911 



-6.93% 



It is clearly seen from this table that the "isotropic" energy curve (IIII.56aj) - 
(IIII.56bl) becomes lowered in the vicinity of the first excited state (which has vf* ~ 1.79) 
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by roughly 6%; and this evidently yields an important improvement of the 8.8% deviation 
of the spherically symmetric approximation from the conventional results (-w 2nd line of 
the table 1 on p. |4*4"|) . Recall here also the failure of the former series expansion (1IV.9P of 
the gauge potential { p ^A ( r ; i?) as it is expressed by the table in App. C. Indeed, in place 
of the unrealistic raising of the energy curve by 8% we observe now a lowering by some 
6%! 

For exemplifying our approximative method of determining the binding energy for the 
selected quantum number n p = 2, we construct the supporting polynomial S^f) as 

S 2 (v) = - 0.000000111342528740459 ■ z/ 10 + 0.00000441158546612208 • v 9 

- 0.0000783955928235642 • z/ 8 + 0.000825779502976810 ■ z/ 7 

- 0.00574182989858913 • z/ 6 + 0.0277779729669865 • z/ 5 - 0.0959826304956328 • z/ 4 
+ 0.238726544780364 • v z - 0.419428837290107 • v 2 + 0.460369057103840 ■ v 

+ 0.015702789114075 

(F.2) 

which can easily be checked to agree with {v) (IV.42I) on the values v — 1, |,2, |,3 
(see the above table for S^{v) with the deviation A{z/}). Next, the maximal value of 
S 2 (v) O is found as 

S 2 \ = 0.2424... (F.3) 

z I max v ' 

at the value v\ of the variational parameter 

z/| 2} = 1.83299... (F.4) 

(see the figure below). Thus the RST prediction for the energy of the first excited 
state (Ep 2 ^, say) is obtained as 

4 2} = ■ 5 2| max = -6, 8029 • 0.2424 . . . [eV] = -1, 649 . . . [eV] . (F.5) 
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Fig. F.I Positronium Energy and Maximal Value of 
the Spectral Function S^{y) (jV.4S\ ) 



From the smallest (but integer) values of 2v (see the table on p J 1 1 41) one constructs the 
interpolating polynomial S^v) (IF. 21) as an approximation to the spectral function S\°^{v) 
(IV.420 . Its maximal value (1F.3P determines the energy e£ 2} (Q) of the first excited 
state (rip = 2 -H- l p = 1). This "anisotropic" result deviates now from the corresponding 
conventional prediction (11.41) by only 3%. Thus the present anisotropic correction improves 
the "isotropic" RST prediction (solid line) from 8,8% deviation (cf. table 1 on p. HI]) to 3% 
deviation (with respect to the conventional prediction, broken line) which is a reduction of 
the deviation of roughly 66%, see also table 2 on p. [77] for the case of higher excited states. 
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The 66% reduction of the deviation (from the conventional results) naturally leads us 
to the question of how large this decrease may be for the higher excited states (l p > 1)? 
In order to get some feeling for this, one may sketch the "correction function" A{z/} (F.l) 
over the variational parameter is, see Fig. F. II below. Here it is satisfying to observe a 
decrease of deviation which ranges from 5,76% for v — 1 up to 11% for v = 100. The 
latter case (is ~ 100) is due to a quantum number of n p ~ 35 (see table 1). It becomes 
obvious from this diagram that the "anisotropic" correction A{z/} is countercurrent to 
the "isotropic" deviation At[^] due to the spherically symmetric approximation Sect. Ill, 
being defined through 



Indeed, this object quantifies the extent of deviation of our "isotropic" RST predictions 
from their conventional counterpart ( II. 4 p and is displayed in the last column of table 1 
on p. HU and when this deviation function Ax(z4 n ') (1F.6P is also represented in Fig. F. II, 
it becomes evident that for increasing is it runs countercurrently to the present correction 
function A{z/} (IF. 1|) . But this then entails the consequence that both curves must intersect 
at a certain value of is (is x , say); and such an intersection means that the conventional 
and RST predictions will agree in the vicinity of that v x (see Fig. F. II below which 
yields v x ~ 30 corresponding to quantum number n p ~ 15). See also the discussion of 
this effect below table 2 on p. [77J 





£tK J ) 



(F.6) 




'conv 
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Fig. F.II Relative Anisotropy Correction A{z/} ( IF. 1\ ) 
and "Isotropic" Deviations At(^) ( E 61 ) 

The consideration of the anisotropy effect in the first approximation order (~ ^ p ^y4'"(r)) 
yields a relative lowering of the "isotropic" RST predictions Et(V) (1III.56ap by A{z/} 
percent, cf. (IF.ip . Since these "isotropic" predictions yield always higher energy than their 
conventional counterparts (11.41) . their lowering by the anisotropy effect yields considerable 
improvement of the RST predictions for the low excited states (table 2 on p. ITT)) . But 
since the "anisotropic" deviation A{i/} increases up to 11% (for v — > oo), the decreasing 
"isotropic" deviation At(4"') (last column of table 1 on p. |4"4"1) must become compensated 
by A{z/} for some value v x of v (here v x ~ 30, i.e. n p ~ 15); and in the vicinity of this v x 
the RST predictions will then agree with the conventional results, see the discussion below 
equ. (1F.6P and table 2 on p. [77J 
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Appendix G: 

Poisson Identity and Generalized Trial Amplitude 



For the numerical exploitation of the "isotropic" Poisson identity ( 1III47I) it is very 
convenient to rescale again both the gauge potential H4 (r) and the generalized trial 



amplitude $(r) (IVI.ip - flVI.3bj) . For the latter object, one puts 



^^l^mrT)-^ (ai) 

(y ■= 2/3r) 

so that the dimensionless trial amplitude <& vp (y) appears in the following form 

®*p(y) = y v (po+piy)e~* ■ (G.2) 

Next, it is necessary to look for the (exact) solution of the Poisson equation f lIH.361) 
when the trial amplitude $(r) in that equation is given by the three-parameter ansatz 
f lVI.lj) - (1VI.3bp . i.e. in dimensionless form ( 1G.1|) - (1G.2[) . The corresponding Poisson equa- 



tion reads explicitly 

d 2 ~a vp {y) 2 d~a vp {y) = (gg ■ y 2 ^ 1 + 2 PoPl ■ y 2 » + p\ ■ gV 
dy 2 y dy ' T(2z/ + 2) 

([ p U (r) = 2/3« s KM) ■ 



(G.3) 



Obviously, the generalized potential d up {y) collapses here to the former potential a u (y) 
(1A.8P if the parameter pi tends to zero, because the present equation (1G.3I) becomes 
in this limit (pi — > 0) identical to the former Poisson equation f lA.ip for the original 
ansatz (lIII.50p - flIII.5ip . Indeed, the present generalized Poisson equation (1G.3[) is closely 
related to its precursor flA.ip . This becomes readily more obvious by exploiting the fact 



that flG.3p is a linear inhomogeneous equation where the inhomogeneity is essentially a 
quadratic polynomial with respect to the variational parameters po,pi- For this reason, 
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My) = Wo „T^ , (G.4) 



the solution a vp (y) may be written as a superposition of three solutions a v (y) of the simple 
precursor (1A.1I) : 

pI ■ a»{y) + ZpoPi ■ a u+ i(y)+p\ ■ a u+1 (y) 
T{2v + 2) 

where the functions a v (y) are (exact) solutions of the (slightly modified) precursor equa- 
tion ([Qjl 

d\J^ + 2 daAy± = y ^ . e-v. (G.5) 

Quite analogously to the former situation ( 1A.17j) one introduces here the dimensionless 
field strength e vp (y) through 

Mv) = (G.6) 
dy 

which then is subjected to the source equation 

deupjy) , 2 (go • y 2 "" 1 + 2poPi • Z/ 2 " + gj • y^ +1 )e^ 

% e ^- r(2z/ + 2) • 

This is the same situation as with the Poisson equation flG.3j) so that the field strength 
e up (y) can also be written as a sum of three terms, i. e. we may put 

p 2 ■ e»(y) + ZpoPi ■ + pi ■ ev+i(y) 

eM = , (G.8) 

so that the simple field strength e u (y) obeys the source equation 

d ^ + ~e v {y)=y 2v - 1 -e- V ■ (G.9) 
dy y 

Comparing now this source equation for the field strength e u (y) to the Poisson equation 
0G.5p for the potentials a u (y) one recovers the field strengths e u (y) as the simple derivatives 



of the corresponding potentials a u (y), i.e. 

da v (y) 

^ ' = dy ' (G.10) 
Consequently, since the solutions of the Poisson equation (1G.5P can easily be obtained 
for integer 2u, cf. (1A.7|) 

(21/ + 1)! f _„ ^ 2u+l-n y"1 
a ^ ) = — I 1 " 6 ^ 21/ + 1 VJ' (G ' n) 



131 



and the solutions of the source equation f lG.9j) for the field strengths e u (y) are obtained 



from this result by the simple differentiation process (IG.lOj) as 



^(^J^/fV), (G , 2) 

cf. (IA.18[) . Clearly, the latter two results for integer 2v may be generalized to arbitrary 
real values of the variational parameter v and do then appear in the following form: 

aM = -r(2„ + 1) • (l - e-y ■ J2 v{2u l 2 + n) y 2 » +n ) (G.13a) 



n=0 
x ,2v+n 



e„(y) = r(2» + 2)e-».j: T -^-— , (G.lSb) 



n=0 



T(2v + 3 + n 



cf. ( [Qjl and 



But now that both the generalized potentials a up (y) (1G.4j) and e up (y) ( 1G.8j) are ex- 



plicitly known, one can realize the Poisson identity (1III.47P within the present framework 
of the generalized trial ansatz (1VI. lj) -( TVI.2[) . First, consider the gauge field energy 
( 1III.45P which may be rewritten in terms of the dimensionless objects as follows 



4 cl = -hc(2a s f3)j^dyy 2 (^^j 2 = -hc(2a s (3) J™ dy y 2 el p (y) . (G.U) 

Inserting here the dimensionless field strengths e up (y) (1G.8I) recasts the gauge field energy 
Eft into the following shape 

E§ = -—(2/3a B )-s pot (u,p , Pl ) . (G.15) 

This is indeed a result which looks very similar to the former case 54[) for the simpler 
trial amplitude ( lIH.50p - flIH.5iP ; but the present potential function e pot {y,pQ,px) emerges 
now as a polynomial of fourth order with respect to the additional parameters po and pi, 
i. e. 

£ P ot(v,Po,Pi) = epot(^)•Po + 4 PoPl• £ I(^)+ 4 PoP^ £: II( zy )+ 4 PoP^^III(^)+P^^Iv(^ / ) , (G.16) 



132 



with the potential functions being given by 

eu(v) 



-pot 



V 



dy y 2 
dyy 2 
dyy 
dyy 
dy y 2 



T(2u + 2) 
e v {y) ■ e v+ i 



(y) 



2 



r(2z/ + 2) 2 

(y) ■ e u+ i (y) + \e u (y) ■ e v+1 (y) 



2 e ^ 



T(2z/ + 2) 2 

r (j/) -eH-i(y) 



T(2z/ + 2) 2 



(G.17a) 
(G.17b) 
(G.17c) 
(G.17d) 
(G.17e) 



J(2i/ + 2) 

Since the field strengths e^y) are explicitly known, all the potential functions (IG.17ap 
( !G.17eP can now be explicitly calculated. For the first one, i. e. (lG.17aj) . one finds 

L T(Au + 3 + m + n) 



-pot 



E 



2^+3 Z_-c 2 m+n r(2i/ + 3 + m) -r(2z/ + 3 + n) 

m,n=0 ' 



(G18) 



2f int. 



2u+l 

E 

m,n— 
(m + n>2) 



(m + n — 2)! 



2m+n- 



m\ nl 



2u+l 

*E 



; n - 2)! 



n=2 



n! 



Here, the result in the first line refers to general (real) values of the variational parameter 
v and the second line presents the specialization to integer values of 2v. Clearly, this first 
potential function e pot {v) (IG.18P is nothing else than the corresponding result ()A.28j) for 
the simple trial function ( 1III.50P which has p± = 0; and therefore the present generalized 
potential function £ po t(^,P0)Pi) ( 1G.16[) must of course collapse to e po t(^) ( ]A.28jl for p\ — > 0. 

But for no n- vanishing variational parameter p ly the other potential functions ej, £pv 
flG.17bp - flG.17eP yield important contributions and must be explicitly known for the ex- 
tremalization process flIII.40p due to the principle of minimal energy (the role played by 
the third variational parameter v 2 is played here by the parameter p l5 cf. flVI.ip - flVI.3bp ). 
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Thus one finds by explicit integration for the first additional potential function £j{v) (1G.17bD : 

( ) 2u + 2 V — T(4u + 4 + m + n) 

£l{U) " 2 4 -+ 4 ^ 2 m + n ' T(2u + 3 + to) ■ T(2z/ + 4 + n) 1 j 

m,n=0 ' 

2 " int - ^ \ , N 1 / 1 2 ^ 1 (2i/ + n)f 

Next the second potential function £11 (f) flG.17c|) obviously splits up into two contri- 
butions where the first one is closely related to the original function e pot (u) flG.17aj) with 
a shift v — > v + I so that this second function appears as 

e u {y) =£n(^) + 4(^) (G.20a) 

/• e,. , 1 (w) ■ e„, 1 (v) 1 
£» = y <W U+ f" 2 l + ; + } r J = (2„ + 2) 2 ■ e pot (^ + i) (G.20b) 



£ h» "\ J d yy 



1 /\,„,2 e ^(^) • e ^+l(^) 

T(2z/ + 2) 2 



(i/ + l)(2i/ + 3) ^ _J_ r(4^ + 5 + TO + n) 



1v int. 



2^+5 2 m+n r(2z/ + 5 + m)r(2z/ + 3 + n) 

m,n=0 ' 

n=2v + 3 

i x - 1 (m + n — 2)1 x ^ 1 ^ 

+ l)(2z/ + 3U 2^ M. l^ n ( n -l)~ 1^ 



n(n — 1) ^— ' nfn — T 

m,n=0 n =2 v y n=2 v ' 

(m+n>2) 



A similar situation does occur also for the third potential S\\\{y) ( IG.17dj) . Indeed, 
comparing this to the first potential E\{y) flG.17b[) one easily realizes that both potentials 



differ essentially by the shift v — y v + ^, i. e. one finds 

e m (u) = (2z/ + 2) 2 -ei{u + i) . (G.21) 

And finally, such a shift z/ — > z/ + 1 of the variational parameter p is also observed in 
connection with the original and the fourth potentials e po t(v) (IG.17aj) and erv^) flG.17el) 
so that one arrives at the following result: 

evfiy) = {2p + 3) 2 (2z/ + 2) 2 • e pot {u + 1) . (G.22) 
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The generalized potential function £ po t{v,p) (1G.16j) is an important object for the 
calculation of the positronium spectrum since it is the essential constituent of the gauge 
field energy Eft (10.15D . The point here is that Eft must obey the Poisson identity 
(1111.47ft ; otherwise the energy functional E[$] (1III.41 [) could not be reduced to its physical 
part (1111.55ft alone which consists solely of the kinetic and potential energy (recall 
here that for the approximative treatment of the anisotropic configurations in Sect. IV 
the Poisson identity is not satisfied so that the energy functional must embrace also the 
Poisson constraint term in addition to the kinetic and potential energy, cf. ( 1IV.29j) ). 
Therefore it is a satisfying reinsurance to see the Poisson identity being satisfied here 
even for the approximative (!) configurations & vp (y) (1Gr.2[) and its associated potential 

aup{y) 4S3D- 

In its rescaled form, the mass equivalent M^c 2 (1111.48ft reads 



, el 2 hc(2a s /3) 



c v(2u + 2 ) J d yy~ a ^(yy^ P (y) (G.23) 



with the trial amplitude <& vp (y) being given by equation (1G.2|) and the corresponding 
gauge potential a vp (y) by equation (1G.4|) . Thus, inserting the latter two objects into the 
present mass equivalent M^c 2 (IG.23P yields obviously a polynomial of fourth order with 
respect to the variational parameters po,p±, just as is the case with the gauge field energy 
(lG.15ft - (1G.16j) . The Poisson identity (1111.47ft demands now that both fourth-order 
polynomials must be identical which implies the identity of the coefficients in front of the 
various products of the parameters Po; Pi- 
First, the lowest-order identification (~ pfy lets emerge the zero-order function £ po t(^) 
(IG.17aj) in terms of the potential a v (y) in the following alternative way 

e p<*(") = -FTTT^^ / dyy 2u+l e-» ■ a v {y) , (G.24) 



T(2z/ + 2) 2 

i. e. for integer 2v 

1 f 1 ^2v + l-n {2u + n)\ 



-pot 



V) 



2u+l\ (2v+ 1)!2^+ 1 ^ 2 n n\ 
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whose generalization to arbitrary v is then given by equation (1A.28jh 

Next, the identification of the linear terms (~ v\vi) lets emerge the first-order function 
e\[y) ( lG.17bj) alternatively in terms of the potential a u (y) as follows 

eiiy) = — F7o~T~ov? \ dy y 2u+2 e~ v ■ &v{y) (G.26) 



T(2z/ + 2) 2 

1 



T(2z/ + 2) 2 
i. e. for integer 2v: 

< \ i 1 1 ^ 2u+l-n (2u+l + n)\ 



(2i/ + 1) ■ (2i/ + 1)! 2 2 ^+ 2 ^ 2™ n! 
2z/ + 2 f 2i/ + l 2 ^2z/ + 2-n (2^ + n)! 



~~ ToTj J_ 92^+1 



2z/ + l (2z/ + 2)!2 2 ^+ 1 ^ 2™ n 

I 71=0 



or more generally for arbitrary z/: 



r s _ 2z/ + 2 1 \- n r(4z/ + 3 + n) y 

Sl{V} ~ 2u + l\ r(2z/ + 3)-2^+3^2"'r(2z/ + 2 + nW 1 ' 



± x it 

r(0u Z 91 . 94^+3 9" ' 



n r(4v + 3 + n) 



r(2f + 2) ■ 2 4l/ + 3 ^ 2 n T(2u + 3 + n) ' 

\ / n=0 ' 

Furthermore, the identification of the quadratic terms (~ "p\ ■ p\ ) confirms the former 
splitting flG.20aj) into two parts s'n{v) and £n(^) with the first part e'n{v) being given 
by equation flG.20bj) . The second part s'diy) is found in terms of the potential a u (y) as 
follows 

e'M = - 2 . T{2 \ + 2? j d yy 2v+ ^~ y - a ^y) ( G - 29 ) 
i 

- 2-T(2z/ + 2) 2 

i.e. explicitly for integer 2v 



dyy 2 » +1 e-y-a u+1 (y) 



v + Dll 1 V 1 (2u + 2 + n), ■ 



(z/ + l)(2z/ + 3) J 1 y^ 2 12z/ + 3-n (2i/ + n)! 

N I 92^+1 / j On 



2v+l (2v)\2 2 "+ l ^ 2 n 2z/ + 3 n\ 

I v 7 n=0 
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For arbitrary v, this result is generalized to 



(i/ + l)(2i/ + 3) 
2i/+ 1 



1 



oo 



n T(4z/ + 4 + n) 
2™ ' T(2u + 2 + n) 




1 - 



r(2z/ + 4) ■ 2 4 »+ 4 



E 



(G.31) 



n=0 



1 



oo 



n T(4^ + 4 + n) 
2» ' r(2z/ + 4 + n) 



(^ + 1)1 



T(2u + 2) ■ 2 4iy + 4 



E 



n=0 



Finally, the last two auxiliary functions em(u) (lG.17djl and £iv(^) (lG.17eh are recov- 
ered by the identification process to undergo the former relationships flG.2ip - fjG.22j) and 
therefore need not be reproduced explicitly once more. Thus the overall result is that 
the total potential function e po t(v,p) is built up by only three "independent" auxiliary 
potential functions s pot (i/) (1A.28[) . s\{y) fjG.28j) and s"i(u) (10.311) . The subsequent table 
displays the values of these auxiliary functions on the relevant values of the variational 
parameter v. 



137 



CO 

GO 



V 


epot(f) 


ei(u) 




em(f) 






(IA.28ll.llG.18ll.llG.25ll 


AG. 191. (IG.21II, (IG.28II 


1G.20al-dG.20cl. (IG.29ll-IIG.31ll 


(|G.2l| 


JG.22} 





1/2 


3/4 


17/8 


25/8 


33/4 


1/2 


5/16 


25/32 


27/8 


231/32 


837/32 


1 


11/48 


77/96 


895/192 


837/64 


965/16 


3/2 


93/512 


837/1024 


6115/1024 


10615/512 


59475/512 


1.7942 


11410011/70398355 


18307008/22197935 


30171115/4469191 


67286322/2579651 


127537747/790294 


2 


193/1280 


2123/2560 


37371/5120 


30927/1024 


101997/512 


5/2 


793/6144 


10309/12288 


53081/6144 


169995/4096 


1290317/4096 


3 


1619/14336 


24285/28672 


572891/57344 


447661/8192 


480429/1024 


7/2 


26333/262144 


447661/524288 


5950665/524288 


9128151/131072 


87461775/131072 


3.758 


0.095054241197409 


0.857009038635841 


12.054442498170000 


78.107662353796800 


788.549642958726000 


5.8740 


0.066169154715619 


0.876608961672518 


17.889125526774400 


166.345557346382000 


2381.422474398040000 


8.1307 


0.050122450469590 


0.890244891770568 


24.190554653983300 


297.684754381852000 


5602.879738648750000 


10.5044 


0.040003065289686 


0.900420995992488 


30.876199420891700 


477.632653223513000 


11254.456630209200000 


20.9538 


0.021304380884411 


0.924772043678155 


60.694395455665500 


1784.286219086780000 


79307.664262054600000 


35.7017 


0.012897315290706 


0.940258135564422 


103.307086974352000 


5068.138932269630000 


374712.658866740000000 


51.9196 


0.009013462373609 


0.949470915666053 


150.512127748210000 


10638.361918298200000 


1131556.682871470000000 


69.3583 


0.006816023482166 


0.955719638189442 


201.500567239779000 


18927.246550543200000 


2673278.362712070000000 


87.8627 


0.005418519284556 


0.960299247611835 


255.774202228416000 


30335.614989561300000 


5407202.467179070000000 


107.3506 


0.004458053255881 


0.963836463562330 


313.064859152140000 


45264.878166797100000 


9832433.121974210000000 


127.8215 


0.003759305734937 


0.966679154371855 


373.353836428638000 


64172.139638711000000 


16566692.535880100000000 


149.3660 


0.003227635139863 


0.969039351893033 


436.897193397184000 


87644.083923089700000 


26402600.563888400000000 


172.1617 


0.002807948458097 


0.971054276427006 


504.212697547382000 


116472.927171035000000 


40397270.423721200000000 


222.5187 


0.002182090441513 


0.974384926745674 


653.142502031965000 


194728.948425671000000 


87151050.944645500000000 


280.9497 


0.001734286200261 


0.977095438492802 


826.236162746333000 


310705.515336457000000 


175365719.699678000000000 


348.2858 


0.001402947397976 


0.979356411386679 


1025.982804702670000 


477935.192145890000000 


334116453.208982000000000 


422.8567 


0.001158190262783 


0.981230916467000 


1247.450242037150000 


705138.107234676000000 


598116579.333488000000000 


501.4654 


0.000978443179650 


0.982771788905577 


1481.148639701370000 


992496.232749923000000 


997902240.046021000000000 



Appendix H: 
Alternative Parametrization of Trial Amplitude 



The normalization condition ( 1VI.2|) for the trial amplitude $&(r) (IVI.lj) reduces the pa- 
rameter space M 2 of the two real- valued ansatz parameters bo, b\ to some one-dimensional 
compact subspace which is topologically equivalent to the circle S 1 . This entails ulti- 
mately that the energy function Et(V, Po,Pi) flVI.26j) must be extremalized on the half- 
cylinder given by the constraint (IVI.2j) together with the condition — | < v < oo. How- 
ever, it may be instructive to use as the underlying parameter space also a half-plane 
{ — | < z/ < oo, — oo < p < +00}. In this sense, one parametrizes the generalized trial 
amplitude $;,(r) (IVI.lj) in the following alternative way 



$ 6 (r) = $,r y (1 + 6r)e" /3r (H.la) 



« (r) = £^. 1 (H ,b) 



vr T(2z/ + 2) w 2 (z/,p) 
where the denominator -u;(z/, p) in (IH.lbj) is given by 



w(i/,p) = v 7 ! + 4(i/ + 1) • p + 2(i/ + l)(2z/ + 3) • p 2 (H.2) 
and the variational parameter p is related to the original ansatz parameter b through 

P*^. (H.3) 



Substituting this back into the original ansatz (1H. la|) for <&&(r), one arrives at the following 
form of the trial amplitude 

«,(r) = Jl , V . %M (H.4a) 
V ; V vr ^(21/ + 2) ™(f,p) 

$^(y)=^(l+ P2/ )e-i . (H.4b) 
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Here, the original normalization condition ( 1III.38j) on $&(r) transcribes to the present 
reduced amplitude $ vw (y) (lH.4bj) as 



1 



T(2z/ + 2) -w(y,p) 



(H.5) 



Kinetic Energy 



After the reparametrization of the trial amplitude is fixed now, one can next turn to 
the kinetic energy ( 1III.42j) -( 1III.44j) in terms of the new parameters p, w. Here, the 
radial kinetic energy ^E^ (1111.430 emerges in terms of the new parametrization as 



(r) E . 



KIN 



w 2 (u,p) 



(H.6) 



with the radial kinetic function ^Ekin(v, p) being defined by 



(r) £KiNO,p) 



1 



T(2i/ + 2) 



dyy 



d$ vw { y) 
dy 



-1 2 



(H.7) 



The explicit calculation of this integral by use of the dimensionless trial amplitude <& vw (y) 
( IFTibj) lets then emerge this radial kinetic function in the following form 



(r) £KIN(>, p) 



+ 



V 



2v+l y 2 V 



1)-P 2 



(H.8) 



4(2i/ + l) 

i.e. a quadratic function with respect to the additional parameter p (1H.3j) . 

Of course, the kinetic energy is an observable quantity and therefore cannot depend 
upon the special parametrization of the wave function $b(r). This means that both 
parametrizations (IVI.lj) -f TVl.3bj) and flH.laj) -f lH.5j) must generate the same kinetic energy, 
i. e. both results flVI.5j) and flH.6j) for the radial type of energy must be identical 



( r) £ kin (/3; v,Po,Pi) = {r) E K m(fc u,p) 



(H.9) 
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But this requirement says now that there must exist a diffeomorphism which relates the 
half- cylinder, being defined by the constraint ( 1VI.2j) together with — | < v < oo, to the 
half-plane {— oo < p < oo; — | < z/ < oo}; and this then must guarantee the identification 
(1H.9|) . The desired map is specified by the following simple transformation 

w = — (H.lOa) 
Po 

p= — , (H.lOb) 

Po 

since by this relationship both kinetic functions ^^(is, p , pi) (1VI.6h (1VI.7[) and ^£kin(^, p) 
(lH.7[) - (lH.8h become linked through 

lr) ( \ ^KIN^p) mi1 x 
1 «kin(f,P0,Pl) => T( ( H - U ) 

w l [y, p) 

so that the required numerical identity flH.9|) can actually be true. 

Similar arguments may also be used in order to set up the longitudinal kinetic energy 
W£kin(/3;M rfTTL44l 

00 

W E^{^p) = ^ll-\j^^ b {r) (H.12) 



e 



2 



-(2a B p) 2 l 2 



W £kin(^, p) 



' 2a B y °' ' v w(v,p) 2 ' 
with the longitudinal kinetic function being found as 

Of course, this kinetic function arises again also through applying the parameter trans- 
formation (lH.10aj) - (lH.10b|) to the former kinetic function ^£kin(^,Po,Pi) (1VI-9j) 

1 fekin(^,PO,Pl) T7 \ • ( H - 14 ) 

w z (v, p) 

Such a transformation must then also hold for the total kinetic function £kin(^, p) 

ejsm(y,p) = {r) e K w(v,p) + W £kin(^,p) • l 2 p (H.15) 



4i/(2i/ + 1) (2i/ + l 



+ 777— rr • P + + C • P' 



141 



i. e. the action of the parameter transformation flH.10al) - flH.10b|) on the total kinetic 
function ekinC^Po; Pi) flYLH]) reads 

£kin(^, p) 



w 2 (u, p) 



(H.16) 



Interaction Energy 

Clearly, the intrinsic consistency of the reparametrization demands that an analogous 
relationship must apply also to the potential energy y— M' e 'c 2 J (1VI. 12[) . The coun- 
terpart of this object (E$, say) is built up by the interaction potential ^A w (r) as the 
counterpart of H4o(r) in the same way as for the first parametrization, cf. flIII.45j) . i. e. 



he f . , fd^A, 







Here, the reparametrized potential ^A w (r) is the counterpart of the former H4o(r) and is 
also the solution of the Poisson equation (IIII.36I) . albeit with regard to the reparametrized 
source-term flH.4ap -f fH.4bp 



or rewritten in dimensionless form (^A w (r) = 2(3a s ■ a uw (y)) 

1 l> 2 (v) 

A v ~ a »M = - TV, ^ 57 S • (H.19) 

where the dimensionless form 3> vw (y) is already displayed by equation (1H.4b[) . This special 
form of the Poisson equation suggests again to compose its solution a vw (y) from the more 
elementary potentials a u (y) (G.13a) as the solutions of the simplified Poisson equation 
(G.5), i.e. we put again 

_ , , au(y) + 2p ■ a v+ i(y) + p 2 ■ a v+1 (y) 

a uw {y) = rfo — j: Y( S H - 20 

T(2u + 2) • w z {u,p) 
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where the elementary potentials a v {y) obey the simplified Poisson equation (G.5) and 
therefore are given by equation (G.ll) for integer 2z/, and by equation (G.13a) for arbitrary 
v (but v > — ~). Substituting this solution a uw {y) (IH.20|) back into the electrostatic field 
energy E$ (1H. 17j) (or equivalently into its mass equivalent) lets then appear the latter 
objects in the following form 

E M = = -£ (2/3ob ) . ^ , (H,!) 

with the alternative potential function £pot(j / ,p) being found as 

£pot(v,p) = £ P ot(f) + 4p • ei{v) + V • e u (i/) + 4p 3 • e m (^ +p 4 • ejy(v) . (H.22) 

The crucial point is here that this alternative form of the potential function just 
guarantees the numerical identity of both electrostatic field energies (G.15) and E$ 
f lfOT]) : 



4 Cl = Et ] (H.23a) 
!N C 2 _ m^c 2 ) . (H.23b) 



Indeed, applying the parameter transformation (lH.10ap -( 1H.10bj) to the original potential 
function e po t(i / , Po,Pi) (IVI-13j) yields the transition 

epot^PcPi) => (H.24) 



and this is just what validates the claimed numerical identity (IH.23a|) . 



Finally, the value of the total energy functional Et |V l on the selected trial configurations 
appears now again in the form of the following energy function 

E T ((3, u,p) = 2E Km ((3, u,p) + E$(P, u,p) (H.25) 
62 t(2Pa B y. £K Z iU ^ -Wa B ). £po ^ p) 



clb { w 2 {u,p) w 4 (z/, p) 

with the kinetic function £kin(^,p) being specified by equation (IH.15P and the potential 
function £pot(^, p) by f)H.22[) . 
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Extremal Configurations 

In order to determine the extremal values of the energy functional E' IV 1 on the selected 
trial functions, one has again to determine the local minima of the total energy E^(f3, u,p) 
( 1H.25I) for any l p . This is equivalent to looking for the minima of the reduced energy 
function Et(z^, p) 

E T {v,p) = -^--S v {v,p) (H.26a) 
Aa B 

3p(u,p) = f OT ^ P l r - (H.26b) 

The local maxima of the spectral function S-p (1H.26bl) determine now the binding energies 
as the local minima of the total energy E<r(V, p) (lH.26ap . 

For a brief demonstration of the present approximation method, one may consider the 
groundstate (l p — =>■ nj> — 1). Obviously, the roughest approximation is specified 
by putting both variational parameters v and p to zero. For this situation, the kinetic 
function e K iN (IH. 15[) becomes simplified to 

£kin(0,0) = ^. (H.27) 

Furthermore, the potential function £pot (1H.22H becomes reduced to 

£ POT (0,0)=£ pot (0) = i, (H.28) 

and the function w(u,p) flH.2j) collapses to unity 

w(0,0) = l. (H.29) 

Consequently, the spectral function S-p{v,p) flH.26bj) adopts the unity 

Sp(0,0) = l, (H.30) 

and therefore the positronium groundstate energy Kr(u,p) (1H.26aj) is found in this rough- 



est approximation as 



e 2 

E T (0, 0) = ~ -6,8029 . . . [eV] . (H.31) 

Aa B 
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This is exactly the value predicted by the conventional theory, cf. (II. 4j) ; but within 
the present RST framework it appears as an approximation] Therefore one expects that 
in the next higher approximation the corresponding RST prediction fore the groundstate 
must yield some lower energy value; and indeed, putting p to zero and thus admitting only 
two variational parameters (i. e. (3, v) yields a groundstate prediction of —7,23055 . . . [eV] 
in place of the conventional prediction flH.31j) . see first line of table 1 on p. HH 



Alternatively, one could have taken as the next higher approximation also v = and 
leaving p as the second variational parameter. In this case, the kinetic function ekin 
( 1H.15|) is simplified for the groundstate (l p = 0) to 

SKm(0,p) = \ + \p 2 (H.32) 



and, furthermore, the potential function £pot (1H.22j) becomes 



£pot(0,p) = epot(O) + 4p • ei (0) + Ap 2 • e n (0) + 4p 3 • e m (0) + p A ■ e w {0) , (H.33) 
i. e. by use of the table in App. G 

b P ot(0,p) = \ + 3p+yP 2 + + f ^ ■ (H ' 34) 



This together with w(0,p) flH. 2j) builds up the spectral function S-p (IH.26bj) to the follow- 
ing form 



Sv(0,p) = S [ ° l (p) = 1 - ^ 9 / M ; ^ J" J . (H.35) 



{1 + Qp + 17p 2 + 25p 3 + f p 4 } 2 
(1 + 2p 2 ) ■ [1 + Ap + 6p 2 
The present alternative parametrization serves as a technical consistency check for 

our numerical program of approximately determining the positronium binding energies 

= E T (u^\p n ); and it is very satisfying to see the positronium energies E Xa up to E xd 

(lVI.40aP - (lVI.40bj) of the compact parametrization (lVI.36al) - (1VI.36bh emerging also in the 

alternative open parametrization (IH.lOajl - flH.lObj) . see the Fig. H.I below. 

Concerning now the excited states, one has to relax the adhoc postulate v = of 

the simplified groundstate treatment and thus one has to admit any real value of v (> 
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— |). This means that one has to look for the maximum of the spectral function S p (u,p) 
(IH.26bl) as a function of the two variables v and p. This may be done by means of 
an appropriate computer program, for the first excited state (n p = 2 <^ l p = 1) one 
finds the maximum S [ p ] (*) = 0, 229683 ... of the spectral function S p (u,p) flH.26bl) at the 
equilibrium point 

vf ] = 1,258290... (H.36a) 
= 0, 677982 . . . (H.36b) 

Thus, the energy E$ of the first excited state is obtained in the present (spherically 
symmetric) approximation step as 

e 2 

M 2] = Sf(x) ~ -6, 8029 • 0, 2296 [eV] 

4 «b p (H.37) 

~ -1,56194 [eV] 

Comparing this to the corresponding first-order approximation (of the spherically sym- 
metric kind) Et(z4 2 ^) (see second line of table 1 on p.(?)) 

E T (z4 2] ) = -1, 55087 [eV] (H.38) 

one observes an improvement (of the deviation from the conventional result -Econv (II. 4p ) 
from 8,8% (last column of table 1) to presently 8,2%: 



EloL -E [ p ] 1 , 70058 - 1 , 56194 



% = 8, 2% (H.39) 



E [2] 1, 70058 

This is a relatively small improvement in view of those improvements due to the regard 
of anisotropy which yielded a corresponding improvement from 8,8% to 3,1% (see first 
line of table 2 on p.(x)). Similar magnitudes of improvement are also found for the higher 
excited states (see table 3 on p.(x)). Thus the conclusion is that the first-order anisotropy 
corrections (Sect.V) are much more important than the second-order corrections of the 
spherically symmetric kind. 
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Fig. H.I Spectral Function S^(p) 



The spectral function Sp (p) (1H.35I) possesses two relative maxima at p\ c = — 0, 161651 . 
and pid = 0,565338... with the corresponding maximal values Sp'lic = 1,033319... 
and iSj?|id = 1, 128194 . . . which by means of (lH.26al) yields the two energy minima E\ c 
(jVI.40al) and En flVI.40b|) . The two relative minima la and lb of Sp\p) are found 
at pi a = and px^ = —0, 634469 . . . which yield the two energy maxima Ei a flVI.39al) 
and £^ib flVI.39bl) of the angular parameterization in Sect. VI. 3. Both asymptotic values 
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of Sp\p) for infinity (p — >■ ±00) do agree 

121 

lim S®(p) = — - 0,6302... (H.40) 

so that the real line (—00 < p < +00) of the present parameterization may be compactified 
to the circle S 1 (0 < a < 2n) of the first parameterization in Sect. (VI. 3). 
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